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ABSTRACT 


This report describes the development of two radiation monitors (RADMONs) for use 
in detecting total radiation dose and high-energy particles. These radiation detectors 
are chip-size devices fabricated in 1.2-jxm CMOS and have flown in space on both 
experimental and commercial spacecraft. They have been used to characterize protons 
and electrons in the Earth’s radiation belts, particles from the Sun, and protons used 
for medical therapy. Having proven useful in a variety of applications, the detector is 
now being readied for commercialization. 




PREFACE 


The dosimeters reported on here operate on two principles. The total dose monitors 
use p-FETs (Field Effect Transistors) whose threshold voltage shifts as radiation- 
induced charge is trapped in the gate oxide. The p-FETs have the unique feature that 
they only need to be powered at the time of measurement and not during exposure to 
radiation. Since the dosimeter is generally measured infrequently, these dosimeters 
require low average power. The particle detectors use a special SRAM (Static Random 
Access Memory) that flips state in response to a high energy particle that strikes a 
powered memory cell. Because this detector is a modified SRAM, the power 
requirement is minimal. 

The radiation detector effort was started in the Product Assurance Technology Program 
at JPL in the 80’s. During the past seven years, these papers and reports were 
generated by workers at JPL and associated institutions. The first radiation detector 
was developed in 1984 and delivered to the Combined Release and Radiation Effects 
Satellite (CRRES) in 1985. The launch of the CRRES was delayed until 1990 due to 
the Shuttle Challenger accident. Since then, dosimeters have been delivered for use on 
MSX, the Clementine spacecraft, the Clementine booster, STRV-lb, Telstar, 
SAMMES-1, Intelsat, and STRV-2. These spacecraft are both experimental and 
communication satellites. All orbit the Earth except for the Clementine spacecraft 
which circled the Moon and now is in heliocentric orbit. In addition to space 
applications, dosimeters have been used to characterize the proton beam at the Loma 
Linda University Proton Therapy Facility. 

Currently, the dosimeter development is aimed at commercialization. The surrounding 
electrons are being placed in the same package as the dosimeters in order to provide 
a digital interface between the sensor and the micro-controller. This will greatly 
facilitate their use by eliminating the need to design the analog circuitry on a separate 
board. In addition, it will allow multiple dosimeters to be placed throughout the 
spacecraft to map the radiation field around the Earth and within the spacecraft. 


v 
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Abstract 

A bench-level test is being developed to evaluate memory- 
cell upsets in a test SRAM designed with a cell offset volt- 
age. This offset voltage controls the critical charge needed 
to upset the cell. The effect was demonstrated using a spe- 
cially designed 2-/im n-well CMOS 4k-bit test SRAM and 
a Po-208 5.1-MeV 0.61-LET alpha-particle source. 


I. Introduction 

The characterization of SRAMs due to ion induced upsets 
is essential to evaluating the susceptibility of SRAMs and 
ASIC registers to Cosmic Ray upsets [1]. This character* 
ization usually requires evaluations at a heavy ion source 
such as a cyclotron. Such evaluations are time consum- 
ing and expensive. In recent years laser pulses have been 
proposed as substitutes for the heavy ion sources. How- 
ever the laser simulations are limited by metal layers that 
frequently block the laser pulses from SEU-sensitive nod* i s 
and by the complexity of the ion-photon calibrations [2j 
[3j. 

In order to overcome these difficulties, six-transistor cell 
SRAMs were designed to be easily upset by a laboratory 
ion source. This simulation closely emulates the Cosmic 
Ray upset of SRAMs. The test SRAM (SEU chip) as 
shown in Figure 1 is a 4k-bit device. In the experiments, a 
Polonium 208 alpha-particle source was used to upset the 
SEU chip. The cells of the SEU chip were designed with 
an offset voltage which controls the critical charge needed 
to upset the cells. The purpose of this paper is to explore 
the potential for using test structures to quickly and inex 
pensively characterize the SEU sensitivity of an integrated 
circuit family. In the next phase of this development, the 
SEU chip's upset rate will be characterized by heavier ions. 
The purpose is to explore the possibility of extrapolating 
the alpha particle results to results obtained by heavier 
ions. Next, this approach will be extended to study the 
upset of latches designed from standard cells. Finally, the 
methodology will be included as a part of an integrated 
circuit qualification procedure for custom ICs. 
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Figure 1: The 4k-bit SEU chip designed to be up- 

set by alpha particles. This chip, designed for a 2-fim 
CMOS n-well process, is 3.84 mm by 4.05 mm and fits 
in a 28-pin dual-in-line package. The substrate doping is 
7.0 ± 1.5 x 10 M cm~ : *. 

II. SRAM Design 

A schematic diagram of the SRAM cell used in these ex- 
periments is shown in Figure 2. The pulsed current 
source is used to model the alpha particle strike. The 
offset volt-age, Voff is varied to alter the critical 
charge required to upset the cell. The cells are always 
biased so that VI -node is low and the V2-node is 
high. This means that the drain diodes Dpi and Dp2 
are reverse biased and susceptible to upsets by heavy 
ions. Since Mpw is “on” the V2-node, voltage is equal 
to Voff, which was varied from V DD = 5 V to 1.9 V. 

In order to emphasize ion hits on the Dn2-diode, its 
area was enlarged to 117 /xm 2 ; whereas, the Dpl-diode was 
designed with a minimum area of 22 /xm 2 . The layout of 
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the cell is shown in Figure 3 and its critical dimensions are 
listed in Table 1 . The dimensions of the MOSFETs were 
chosen so that the cell can be operated as a “standard” 
SRAM when Voff — Vnn- 

In order to minimize power line noise, the V^o, GND, 
and V, >| lines were made as wide as possible. This re- 
duces the possibility that one cell will trip another cell 
from the current transient demanded by the flipping cell. 


Table 1 : Dimensions of the transistors and drains in the 

memory cell. 


DEVICE 1 L(jxm) [ WQxm) | Ad (Min 




III. SPICE Analysis 

The principle of operation of the cell is illustrated by the 
transfer curves shown in Figure 4. When the cell is biased 
with V ( >,,,, above 1.9 V, the cell exhibits two stable states. 
For V o,,,, less than 1.9 V, the cell has only one stable 
state. As the critical V ( jj F is approached, cells begin to 
flip according to the distribution in cell parameters. The 
voltage span for this spontaneous flipping action is about 
30 mV for the cells in the SEU chips studied here. 

The critical charge, Q r , was determined for this cir- 
cuit using SPICE and MOSIS supplied parameters. The 
Vl-node and V 2 -node capacitances were modeled by fixed 
metal and polysilicon interconnect capacitances and by the 
drain depletion capacitances given by their areas and pe- 
ripheries and by the SPICE LEVEL 2 junction capacitance 
parameters. A current generator was placed on the V 2 - 
node as shown in Figure 2 and a triangle pulse with a 119 
riserfall shape was used to upset the cell. For a given pulse 
height, transient simulations were evaluated up to 100 ns 
where the response was compared to Vdx >/2 to determine 
if the cell had flipped. The pulse height was adjusted using 
the binary search algorithm until the difference in charge 
(area under the pulse) between successive simulations dif- 
fered by less than 0.1 fC. The resulting critical charge Q c 
versus Voff is shown in Figure 5. The results were invari- 
ant with current pulse widths between 20 ps and 500 ps. 
For pulse widths (PWs) in excess of 500 ps, critical charge 
curves deviated from the 500 ps curve. The deviation of 
the PW 2 ns curve from the PW < 500 ps curve is 
shown in Figure 5. Thus, the cell time constant exceeds 
500 ps. Since the alpha particle pulse width is about 200 ps 
[4], the response of this circuit exceeds the alpha particle 
pulse. As determined from this analysis, the details of the 
pulse shape is unimportant which simplifies the analysis. 


Figure 2 : Memory cell schematic diagram showing the 

placement of y 0FF and the bloated n-drain, Dn2. The di- 
mensions of the transistors and drains are listed in Table 1 . 
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Figure 3: Memory cell layout, designed for 2-fim CMOS 

n-well process, is 42.0 /xm by 45.0 /xm. The alpha sensitive 
region is the bloated n-drain in the p-substrate. The offset 
voltage, V off , is used to enhance its upsetability. 
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Figure 4: Memory cell transfer curves showing the effect 
of Vi >|.-|- on the stable points. 

IV. Alpha Particle Behavior in 
Silicon 

The energy loss of an alpha particle passing through Sil- 
icon is shown in Figure 6 [5], This figure shows that a 

5 MeV alpha particle has a range of 23.7 /im. The differ- 
ential of the energy curve leads to the Linear Energy, LET, 
curve shown in Figure 7. This shows that as the particle 
penetrates the Silicon its LET increases up to the Bragg 
peak and than falls to zero. The integral of the LET curve 
leads to the silicon deposited charge, Q s , curve shown in 
Figure 8. This curve is based on the electron-hole pair 
production rate of 3.6 eV/e-h pair. The charge deposition 
curve shown in Figure 8 indicates that a 5-MeV alpha par- 
ticle can deposit about 225 fC in Silicon assuming no over 
layers. This amount of charge is about 10 times larger than 
the charge needed to upset the cell as shown by the SPICE 
analysis presented in Figure 5. Thus it should be antici- 
pated that the enlarged diode, Dn2, will be upset by alpha 
hits that strike a considerable distance from the diode. In 
a later section the lateral hit distance was determined to 
be 6.2 fim. 

V. Experimental Results 

In the following experiments a 1 jxCi Po-208 alpha particle 
emitter was used. This source was placed 32.5 mm above 
the SRAM cell. This distance is very close to the range 
of alphas in air of 35.2 mm. The cell and the source were 
placed in a vacuum chamber. The flux at the SRAM cell, 
measured with a pin diode, was 214 alphas/cm 2 • sec. The 
partial vacuum pressures were measured using a Wallace 

6 Tiernan FA 112 pressure gauge. The offset voltage was 
applied for a duration that allowed 1 to 30 upsets to occur, 
and this measurement was repeated 30 times and averaged 
to give results with deviations less than 18 percent for data 
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Figure 5: Memory cell SPICE analysis for the critical 

charge. 
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Figure 6: Energy profile of a 5-MeV alpha particle in 

silicon with a range of 23.7 p,m. 






larger than 0.01 counts/sec and 6 percent for data larger 
than 0.2 counts/sec. For this experiment the “on” time for 
the offset voltage varied from 1 to 60 seconds. 

The response of the cell follows the detector equation 
where it is assumed that all the memory cells are initially 
placed in the same state: 

^ = 4>a(N, - AT) (1) 

where N is the number of flipped cells at time t, N, = 
4096 is the total number of cells in the SRAM, <f> is the 
alpha particle flux, and cr is the cell cross section. The 
cross section is determined from the initial slope for N 
approaching zero: 



1 dN 
<f>N t dt n 


( 2 ) 



Figure 7: Linear energy profile of a 5-MeV alpha particle 

in silicon. 



Figure 8: Deposited charge profile of a 5-MeV alpha 

particle in silicon. 


Figure 9: SEU chip upset rate characteristics induced 

by Po-208 alpha particles 32.5 mm above the SRAM in a 
vacuum. The cross section observed by extrapolating the 
peripheral hit and tail regions agrees with the designed 
cross section of 117 pm 2 . Repeatability of the measure- 
ments is demonstrated by the overlap of the data taken 
from 1.9 to 2.7 V (squares) and data taken from 2.7 to 
1.9 V (circles). 


VI. Vacuum Test Results 

The response of the SRAM is shown in Figure 9 for alpha 
particle irradiation in a vacuum. Three regions are ap- 
parent in the figure. Region I, the spontaneous region, is 
governed by cell parameter distributions and indicates that 
the transfer curves, as seen in Figure 4, no longer overlap 
(Voff < 1.922 V). Region II is termed the peripheral hit 
region because the cross section exceeds the design cross 
section of 117 pm 2 . For this cell the extent of the periph- 
eral hit area is 540 pm 2 . Assuming the bloated drain is a 
square 117 /im 2 and the peripheral hit region is a square 
540 /im 2 , then the lateral hit distance is 6.2 /im. Region III 
is the tail region; in this region the upset rate falls with a 
slope that is identical to Region I which indicates that this 
region is also governed by the spread in cell parameters. 

As seen in Figure 9, a change in Voff » AV<>| i of 0.5 V 
is required to bring the cells from the spontaneous flip re- 
gion where the critical charge, Q, is zero to a point where 
the alpha particles are just able to trip the cells (tail re- 
gion in Figure 9). By aligning the V r) , value in Fig- 
ure 9 where the spontaneous-flip line intercepts the phys- 
ical cross-section line (Voff — 2 V) with the origin of 
Figure 5 (Voff — 1-9 V), then AV u r f = 0.5 corresponds 
to Voff = 2.4 V in Figure 5 and to a critical charge, 
AQ C = 58fC. This leads to a collection depth, Dc, of 8 pm, 
as determined from Figure 8. The analysis for D c depicted 
in Figure 8 does not include the effect of over layers. The 
inclusion of over layer in this analysis, requires a detailed 
knowledge of the thickness and composition of the over 
layers. Such information is often unavailable. However, an 
estimate of the effect of over layers on the calculation of 
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D c C an be obtained by assuming the over layers have the 
same density as silicon. Using this assumption, the deter- 
mination uses the curve in Figure 8. The over layers add 
a few micrometers to the point at which AQ C is applied to 
the curve. Since the curvature of the curve of Figure 8 is 
slight, the error in D t is similarly small. 



Figure 10: SEU chip upset rate characteristics induced by 

Po-208 alpha particles 32.5 mm above the SRAM in a par- 
tial vacuum. The distortions that appear in the curves for 
pressures above 0.75 atm are due to the straggle induced 
by the passage of the alpha particles through air. 


VII. Partial Vacuum Test Results 

The response of the SRAM, as a function of partial vacuum, 
is shown in Figure 10. The region II behavior is affected by 
alpha particle straggle. This behavior is noticeable for 
pressures in excess of 0.75 atm. This behavior can be 
explained qualitatively as shown in Figure 11. For P = 0, 
The alpha particles penetrate 23.7 urn but as determined 
from the vacuum results, only 8 pm are effective in 
supplying charge to the junction. As the air is introduced 
into the chamber, the energy of the alpha particles as they 
enter the silicon is reduced which reduces the penetration of 
the alpha particles into the Silicon. For this experiment 
where the alpha particle source is placed approximately at 
the range of the alphas in air, the upset rate goes to zero at 
atmospheric pressure. 

The data shown in Figure 10 is analyzed in Figure 12. 
The V () j y values were determined from the data shown in 
Figure 10 at an upset rate of 1 count/second. The Voff 
values were converted to critical charge using the SPICE 
results shown in Figure 5. The experimental results are 
plotted in Figure 12. 

A simple model was fitted to the data shown in Fig- 
ure 12. The model parameters assumed an over layer cf 
2 /xm and a collection depth of 7 /xm. For pressures be- 


low 0.65 atm the range of the alphas exceeds the collection 
depth. In this region, the range of the alphas in silicon 
decreases as the pressure increases because the incident 
energy of the alphas decreases with increasing pressure. 
Thus, the alphas deposit more charge in the collection 
depth as the pressure increases and the critical charge in- 
creases with increasing pressure. At the peak in the model 
curve, the collection depth equals the alpha particle range. 
For pressures above 0.65 atm, the alpha particle range is 
less than the collection depth. Thus the deposited charge 
falls rapidly to zero. 

The model used in this analysis does not consider alpha 
particle straggling. At the end of its range, alpha particle 
straggling is about 1 /xm longitudinally and 1 /xm in the 
transverse direction [6]. A more complete model must take 
particle straggling into account. However, it is encouraging 
to see that a simple model provides a qualitative explana- 
tion of the data. 

A simple theoretical model was fitted to the data as- 
suming an over layer of 2 /xm and a collection depth of 
7 /xm. This model assumes that the collection depth 
“slides” along the critical charge curve (Figure 5) as the 
pressure changes from vacuum to atmospheric. As the col- 
lection depth slides along the curve the LET increases and 
so the AVoff increases. Once the end of range is reached, 
the charge is reduced and thus AVqff decreases. 


P = 0 0.2 0.4 0.6 0.8 1.0 atm 
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Figure IT. Alpha particle range variations due to energy 
loss in air before particles reach the silicon. 


VIII. Conclusion 

In conclusion, this test SRAM has been made sensitive to 
alpha particles through the use of a cell offset voltage and 
this has allowed a bench-level characterization in a labo- 
ratory setting. The experimental data was linked to alpha 
particle interaction physics and to SPICE circuit simula- 
tions through the alpha particle collection depth. The col- 
lection depth was determined by two methods and found 
to be about 7 /xm. In addition, alpha particles that struck 
outside the bloated drain were able to flip the SRAM cells. 
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Figure 12: Variation of the critical charge with pressure. 
For a collection depth of 7 pm, Q c increases due to an 
increase in LET with increasing pressure followed by a de- 
crease in Q c due to the collection of alphas at the end of 
their range. 


in integrated circuits, Proc. 10th Annual Int. Relia- 
bility Phys. Symp., Orlando, FL, Apr. 7, 1981, pp. 
38-42. 

[5] J. F. Ziegler, Handbook of Stopping Cross-Section 
For Energetic Ions In All Elements (Pergamon Press, 
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Wagner, N. Bordes, and J. M. Bradley, CMOS VLSI 
Single Event Transient Characterization , IEEE Trans. 
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This lateral charge collection was observed to be more thrn 
6 pm. This effort will be extended to the evaluation of 
more complex SRAMs and latches. 
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ABSTRACT 

This paper describes a methodology for predicting the 
SEU susceptibility of a standard-cell D-latch using an alpha- 
particle sensitive SRAM, SPICE critical charge simulation 
results, and alpha-particle interaction physics. Measurements 
were made on a 1.6-pm n-well CMOS 4k-bit test SRAM 
irradiated with an Am-24 1 alpha-particle source. A collection 
depth of 6.09 pm was determined using these results and 
TRIM computer code. Using this collection depth and SPICE 
derived critical charge results on the latch design, an LET 
threshold of 34 Mev cm^/mg was predicted. Heavy ion tests 
were then performed on the latch and an LET threshold of 41 
MeV cm 2 /mg was determined. 


I. INTRODUCTION 

The characterization of digital ICs to heavy ion induced 
upsets is essential for qualifying their use in critical systems to 
be used on spacecraft. These characterizations are usually 
carried out at a heavy ion source such as a cyclotron. Such 
sources are expensive to use and there is usually a large lead 
time in scheduling time on the source. 

The purpose of this paper is to introduce a technique 
utilizing test structures to quickly and inexpensively 
characterize the SEU sensitivity of standard cell latches 
intended for use in a space environment. This bench-level 
approach utilizes alpha particles to induce upsets in a low LET 
sensitive 4k-bit test SRAM. This SRAM consists of cells that 
employ an offset voltage to adjust their upset sensitivity and an 
enlarged sensitive drain junction to enhance the cell's upset 
rate [1,2]. 

We demonstrate that SRAM alpha particle data can be used 
to predict the latch LET threshold. The latch is modified by 
equipping it with an offset voltage so that its LET threshold 
can be lowered into the alpha particle range for direct 
comparison to the SRAM. The normal latch, with its offset 
voltage set to Vdd = 5 volts, is then tested with heavy ions to 
verify the SRAM LET threshold prediction. 

© IEEE 1991 


II. TEST SRAM AND D-LATCH DESIGNS 

The test SRAM and latch designs studied in this paper 
were submitted to the MOS Implementation System (MOSIS) 
and fabricated at a 1.6-p.m n-well double-metal CMOS/bulk 
foundry. A schematic diagram of the test SRAM cell is 
shown in Figure 1. This figure does not include the read/write 
transistors. The pulsed current source is used to model an 
alpha particle strike on drain Dn2 when calculating the critical 
charge of the cell with SPICE. This cell differs from that of a 
standard six-transistor SRAM cell in three ways: (1) the source 
of the p-MOSFET, Mp2, is connected to an adjustable offset 
voltage, V 0 , instead of Vdd 10 provide a control of the cells 
critical charge; (2) the drain area of n-MOSFET Mn2, Dn2, 
has been enlarged by a factor of four over minimum to 
enhance upset rates, thus reducing measurement time; and (3) 
the cell is imbalanced by widening Mn2 over minimum to 
enhance its SEU sensitivity versus V 0 response. The 
dimensions of the MOSFETs and their drains are given in 
Table 1. The bloated drain, Dn2, is fabricated in a p-type 
substrate doped to 5xl0 14 cm' 3 . This substrate doping value 
was determined through a dopant profiling measurement 
utilizing a large edgeless n-MOSFET. The Dn2 junction 


V DD = 5V 



Figure 1. Schematic diagram of test SRAM cell showing the 
placement of V 0 and the bloated n-drain, Dn2. The 
dimensions of the transistors and drains are listed in Table 1 . 
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Table 1. SRAM Memory Cell As-Drawn Dimensions. 


DEVICE 

L(nm) 

W (pm) 

AD (ttm2) 

Mnl 

1.6 

2.4 

17.92 

Mn2 

1.6 

3.2 

74.88 

Mpl 

3.2 

2.4 

12.16 

Mp2 

3.2 

2.4 

12.16 


has enhanced charge collection due to funneling compared to 
charge collected by the p-MOSFET drain formed in the n-well, 
where truncation of the ion-induced plasma track occurs. 

In operation all the memory cells are written into a 
"sensitive" state where Mn2 is turned OFF and Mp2 is turned 
ON, connecting V 0 to the bloated drain, Dn2. V Q is then 
lowered from Vqd = 5 V for a period called the stare time. 
Thereafter V 0 is returned to \foD and the cells read to 
determine the number of upsets. This cycle is repeated at 
different values of V 0 and can be repeated a number of times 
at a given V 0 to improve the resolution of the measurement. 

The latch results presented in this paper were obtained 
from a 64-cell modified transparent D-latch array fabricated 
through the same foundry as the test SRAM. The schematic 
diagram of this latch in the zero state showing the reverse 
biased SEU-sensitive drain diodes is shown in Figure 2. The 
offset voltage, V 0 was added to sensitize the latch to alpha 
particles. The source of Mp2 is tied to Vdd in the unmodified 
latch. The diodes designated with solid squares collect the 
most charge when struck by an ionized particle because they 
sit in the p-substrate. The layout of the latch cell is shown in 
Figure 3. This layout includes an inverter that generates an 
enable-bar, and inverters that buffer the VI -node and V2-node. 



Figure 2. D-latch in the zero state. The reverse-biased diodes 
are SEU sensitive: p-substrate diodes (solid squares) and n- 
well diode (solid diamond). The zero-biased diodes are not 
SEU sensitive: p-substrate diode (open square) and n-well 
diodes (open diamonds). 


The sensitive diode areas at the V2-node and V3-node that are 
designated with solid squares in Figure 2, are outlined with a 
bold line in Figure 3. These diodes each have an as-drawn 
area of 69 |im2. 

III. SPICE ANALYSIS 

The critical charge, Q c of the test SRAM cell and of the D- 
latch cell were determined as a function of V 0 using MOSIS 
supplied parameters and SPICE. The parasitic nodal 
capacitances were modeled by fixed metal and polysilicon 
interconnect capacitances and by the drain depletion 
capacitances given by their areas and peripheries and by the 
SPICE LEVEL 2 junction capacitance parameters. A triangle 
current pulse with a 1:19 rise:fall shape was used to upset 
these cells. For a given pulse height, the transient simulation 
was evaluated to 100 ns where the response was compared to 
Vdd/ 2 to determine if the cell had flipped. The pulse height 
was adjusted using the binary search algorithm until the 
difference in charge (area under current pulse) between 
successive simulation runs differed by less than 1 fC. The 
resulting critical charge versus V 0 response for the V2-node 
of the test SRAM cell and for the average of the V2-node and 
V3-node for the D-latch is shown in Figure 4. These results 
were found to be invariant with current pulse widths up to 500 
ps. Since the alpha particle width is about 200 ps [3,4], the 
response of these circuits exceeds that of the alpha particle 


Vo 



Figure 3. Layout of standard cell latch showing sensitive 
drains (solid squares in Figure 2) and offset voltage V 0 
connection. In the unmodified latch the metal 1 V 0 trace is 
connected through a via to the metal 2 Vdd bus. 
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Figure 4. SRAM and latch critical charge calculated using 
SPICE. 


current pulse. The critical charge of the latch V2-node and 
V3-node, which are connected through the ’’on" MOSFET 
Mn4, were averaged because their values differed little and it 
simplifies the LET-threshold analysis. 

As V 0 decreases the critical charge of the SRAM cell and 
D-latch decreases approximately linearly and goes to zero at 
Vo = V om . For the SRAM cell V om S * 1.8 V and for the 
latch V om L = 2.4 V. DC SPICE simulations of these bistable 
circuits show that the sensitive state and metastable state 
approach each other as V 0 decreases until they are equal when 
Vo = V 0 m* For Vo ^ V om the circuit spontaneously goes to 
the flipped state and is no longer bistable. V 0 m is a function of 
the threshold voltages and geometries of the n- and p- 
MOSFETs. Process induced dispersion in the MOSFET 
parameters give rise to a dispersion in V 0 m for the array of 
SRAM cells. The slopes of the critical charge curves in 
Figure 4 are defined as their upset capacitances, Cu = 
SQc/8Vo* The upset capacitance of the SRAM cell Cus = 56 
fC/V and of the latch cell Qjl = 850 fC/V. 

IV. TEST SRAM ALPHA -PARTICLE RESPONSE 

The response of the SRAM to alpha particles from a 4.6 
jiCi Am-241 source placed 0.508 cm above the SRAM in air is 
shown in Figure 5. This source-to-SRAM distance was 
chosen to minimize test time while maintaining essentially 
normally incident alphas [2], Also shown, is the spontaneous 
response which was measured without the alpha-particle 
source present. As V 0 is lowered from Vdd» the SRAM 
begins to upset at V 0 < 3.2 V in a region termed the tail 
region. As V 0 is further lowered, the upset rate increases more 



Figure 5. SRAM upset rate versus V 0 response to a 4.6 
pCi Am-241 source. The data used in this analysis were taken 
on 5-mV centers and are an average of 200 samples with a 2- 
second stare time. 

slowly; this region is termed the peripheral hit region. Finally, 
as V 0 reaches V om = 1.8 V, the upset rate increases very 
rapidly in what is termed the spontaneous flip region. The 
spontaneous response remained constant (negligible total dose 
shift) throughout the 33.3 hours of alpha particle exposure. 

In the peripheral hit region, cells upset when the alpha 
particles hit the Dn2 drain or come close to, but miss, the Dn2 
drain. This effect can be explained by a delayed field- 
funneling effect [4] and to diffused charge that is collected by 
Dn2. Since the cell critical charge decreases as V 0 decreases, 
the alpha particles can strike further from the Dn2 depletion 
edge and still upset the cell. 

The upset rate R, displayed in Figure 5, was obtained by 
dividing the average measured number of upsets N a that 
occurred in a stare time t by t, or R = Na/t, where 



and the number of samples, k = 200. The stare time t, in this 
case, was fixed at 2 seconds. The spontaneous upset data was 
also divided by t so that both responses saturate at Nt/t= 2048, 
where Nt is the number of cells in the memory array (4096 in 
this case). In the analysis to follow, the spontaneous response 
is approximated by the step function N t [l-u(V 0 -V 0 m)] = 
4096[l-u(V o -1.8)], where u(V) is the unit step function. This 
is justified for the SRAM data shown in Figure 5 because the 
slope of the spontaneous response is 2048 cells/34 mV- sec 
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which is much steeper than the SRAM alpha-particle response 
in the tail region. The response of the SRAM cells is given by 
the detector equation: 

^ = <*t>(Nt-N) (2) 

where a is the device cross section (Dn2 drain area), and <t> is 
the total integrated flux whose particles deposit energies 
exceeding the critical energy, E c required to trip a cell. The 
homogeneous solution of (2) with initial condition at t = 0 of 
no tripped cells (N = 0) is: 


N = N t ( 1 - e‘ a<I>t ) 


(3) 


The upset capacitance Cu. the spontaneous upset voltage 
V 0 m. and the hole-electron pair charge-energy factor for 
silicon K = 44.2 fC/MeV enable one to calibrate the offset 
voltage, V 0 axis in terms of the critical energy by: 



Ec = 


Cy 

K 


(Vo-Vom) 


or for this SRAM, E c = 1.267(V 0 - 1.8) MeV. 


(4) Figure 6. Relationship between measured upset rate, the 
spontaneous density function, the source spectra, and the 
SRAM critical charge response. 


The effective flux, <t> is a function of V 0 and is given by: 


oo 

<t>(V 0 ) = / <KV) dV (5) 

Vo 

where the source spectral flux 4>(V) = «J»(E) and the energy 
E = (Cu/K)(V - Vom)- Applying (5) to (3) and solving for 
dN/dV 0 one gets: 


dN 

dV 0 


dR 

dV 0 


= -at 


( 

N t j 


<KV 0 ) 


= -at<KV 0 ) ;N«N t (6) 

This is the condition under which the data in Figure 5 was 
taken, therefore dR/dV 0 is the average particle spectrum in the 
collection region of the sensitive drain Dn2. 


Figure 6 illustrates the relationship between upset rate R, 
the SRAM spontaneous density n, the source spectra <j>, and 
the SPICE derived critical charge curve. The spontaneous 
density, n, is the derivative of the spontaneous response data 
shown in Figure 5, and was found to be gaussian with a mean 
and standard deviation of V om = 1.8 ± 0.01 IV. The area 
under the n(Vo) curve is equal to Nt. The critical charge 
curve of Figure 6 , characterized by (4) when there is no 
dispersion in n, functions as a calibration for mapping the 
source spectra onto the V 0 axis. In this case, V 0 = V om = 
1.8 V represents zero energy. As V 0 is lowered from Vdd = 
5 V, no cell upsets occur until <I> > 0 which occurs when 


V 0 <3 V. As V 0 is further reduced below 3 V, the upset rate 
increases giving the tail-region response in Figures 5 and 6. 
When V 0 reaches 2.6 V, all of the source spectrum is effective in 
upsetting SRAM cells and the upset response saturates at a 
maximum. In the absence of the peripheral hit effect, the 
response would remain at this maximum value as V 0 was 
decreased below 2.6 V, as shown in Figure 6. 

The tail-region response shown in Figure 5 was smoothed 
and differentiated using a twenty-five point least squares 
difference algorithm [5,6]. This differential response, shown 
in Figure 7, is composed of a number of peaks plus the onset 
of the peripheral hit response. These peaks are due to to the 
selective degradation of the alpha-particle energies as the 
particles pass through the various regions found in the bloated 
drain, Dn2. The major peak, shown on Figure 7, is due to 
the alpha particles that pass through the thickest regions, those 
covered with metal 2 (M2). This M2 region is 38.2 pm 2 [2]. 
In the region covered with M2, the alpha particles lose more 
energy than in other regions so that they have the highest 
linear energy transfer (LET). The cross section for the M2 
portion of the bloated drain was calculated using half the total 
alpha-particle flux , d>(V om )/2 and found to agree closely with 
the as-drawn area [2]. With V 0 set at the peak, V 0 p, only half 
of the alpha particles have sufficient enerev (or LET) to flip 
the memory cells, thus O(V 0p ) = 0(V om )/2. 

A schematic view of the path of the alpha particle through 
the silicon is shown in Figure 8. Relating the full-width-half- 
maximum of the response shown in Figure 7, V 0 f to the 
energy straggle of the alpha particles in the collection region 
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Am-241 Source 


X 0 

Xl 

X 2 

X 3 

x 4 

X 5 


Au Am-241 Source Coaling 


Air 
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Si Collection Layer 


End of Range Layer 


E 0 = 5.476 ± 0 MeV 

5X 1 = 1.59 pm 

Ei = 4.75 ± 0.057 MeV 

6X2 = 508 pm 

E 2 = 4.70 ± 0.057 MeV 

6 X 3 = 7.59 pm ( 4.3 pm ) 

E 3 = 3.46 ±0.076 MeV 

5 X 4 = 6.09 pin 

E 4 = 2.23 ± 0.087 MeV 

5X5 = 0.01 pm 

e 5 = 0 


Figure 8 . Schematic view of the alpha-panic le path from the 
source to the SRAM. The over-layer thickness 5X3 = 7.59 pm 
corresponds to the Dn2 region covered by metal 2. An 
average over-layer thickness for the whole Dn2 area of 4.3 pm 
was obtained. 


Figure 7. Test SRAM differential upset rate response showing 
peak at V 0 p = 2.78 V with full- width-half-maximum of V 0 f = 
0.170 V. 

of the silicon, the Si thickness X 4 - X 2 = 13.58 pm was 
obtained [2], The full-width-half-maximum, V 0 f in Figure 7 
is related to the one standard deviation of a gaussian, V oa in 
Figure 6 by V 0 f = 2.35Voa- Next the Si collection layer, 8 X 4 , 
was determined by matching the charge deposited by the alpha 
particles in 8 X 4 to the shift in the offset voltage, 8 V 0 = (V 0 p - 
V om ) ± V oa = 0.98V ± 0.07V [2]. Finally the over-layer 
thickness, 8 X 3 was calculated by 8 X 3 = (X 4 - X 2 ) - 8 X 4 = 
7.59 pm. This over-layer thickness is applicable to 50% of the 
Dn2 drain area (that covered by metal 2). An effective Dn2 
over-layer thickness of 4.3 pm was obtained using an average 
V 0 p of 2.66 V. This value of V 0 p was determined by first 
obtaining a reduced differential response where the peripheral 
hit response was removed (fitting an exponential to the 
dR/dVo response between V 0 = 2 V and 2.2 V and then 
subtracting this exponential function from the dR/dV 0 data 
shown in Figure 7). Then the average V G p is the weighted 
average of the reduced dR/dV 0 response. 

Using the offset voltage shift and upset capacitance, the 
SRAM-cell collected charge C fcs = QjS - 54.9 ± 3.9 fC. 

Also, the energy deposited in the silicon collection layer 8 X 4 , 
and its straggling ^4 shown in Figure 8 is: 5E = (E 4 - E 3 ) ± 
E a 4 = Q c s/K = 1 .24 ± 0.088 MeV. 


V. D-LATCH ALPHA PARTICLE RESPONSE 

The alpha-particle response of the latch array was obtained 
in the same manner as for the SRAM except that the Am-241 
source was placed 0.558 cm above the latch chip in air. The 
results, shown in Figure 9, include the spontaneous data 
obtained before and after exposing the part to alpha particles. 
The shift observed in these responses is due to total dose 
effects such as the charging of gate oxides. The average of the 
pre- and post-radiation spontaneous data was used in the 
analysis of this part This small shift, due to total dose effects, 
was neglected for the SRAM because the 5V 0 shift in the 
SRAM alpha-particle response was much larger than that of 
the latch. 

The dispersion in the latch alpha-particle response, shown 
in Figure 9, is primarily that of the spontaneous response 
because the offset voltage dispersion due to the alpha-particle 
energy spectrum is much smaller than the dispersion in the 
spontaneous response. This is a consequence of the reduced 
energy detection sensitivity of the latch (K/QjL = 0 052 
V/MeV for the latch versus K/QjS - 0.789 V/MeV for the 
SRAM). The dispersion in the latch alpha-particle response, 
V oa = 34 mV, was found by Fitting the data to a gaussian 
function. As shown in Figure 9 the shift in the offset voltage, 
5V 0 , was obtained between the points on the average 
spontaneous curve and the alpha-particle curve where, on the 
average, half of the latch cells are flipping (Nt/2 = 32). Using 
the SPICE derived upset capacitance of the latch, Cjjl = 850 
fC/V, the collected charge required to upset a modified latch 
cell with its LET threshold lowered into the alpha-particle 
range is 8Q c l = 50.2 ± 28.9 fC which agrees with that 
obtained with the SRAM, Q c s. 
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Figure 9. D-latch alpha-particle response. The data used in 
this analysis were taken on 5-mV centers and are an average of 
10 samples with a 10-second stare time. 


When the offset voltage for the SRAM and the latch is 
swept through the alpha particle 8 E peak, the charge required 
to upset their cells is equivalent. Cfc = K 8 E where 5E is the 
energy deposited in the silicon collection layer 8X4 which is 
the same for the SRAM and latch because they are both 
fabricated through the same process (p-type substrate doped to 
5X10 1 * cm-3). For the same alpha-particle beam the charge 
loss mechanisms (X4 is less than the alpha particle range) are 
equivalent in the SRAM and the latch. Since Q;L = Q cS- * s 
valid to apply the collection depth 8 X 4 obtained with the 
SRAM to the analysis of the latch. 

For the latch, the LET threshold for Vdd = 5 V operation 
can be predicted from the SPICE derived upset capacitance, 
CuL. the measured mean spontaneous offset voltage, V om , 
and the collection depth obtained from the SRAM alpha- 
particle response, 8 X 4 : 

. „„ C(jl (Vdd ■ v om) MeV cm2 

LET = Tr — c-J = j4 (/) 

K p 8 X 4 mg 


where, Cul = 850 fC/V, V om = 2.49 V, K = 44.2 fC/MeV, p 
= 2320 mg/cm3 for Si, and 8 X 4 = 6.09 x 10-4 cm. 


VI. D-LATCH HEAVY ION RESPONSE 

Heavy ion tests were performed with 317 MeV 1-127, 
279 MeV Br-79, and 255 MeV Fe-56 at the Brookhaven 
National Laboratory on the latch array with Vdd = V 0 = 
5 V. These measurements were made by gating the ion beam 
on and then off when a fluence, F = 1.13 x lO^/cosO was 
reached, where 0 is the angle if incidence. This fluence was 
chosen to maintain an experimental error of less than 3.2% 


and caused V om to shift by 0.4 mV for each heavy ion beam. 
The stare time was set so that an average of two upsets would 
occur each test cycle. This assured that essentially all N( = 
64 latch cells in the array were susceptible to upset and 
allowed the effective cross section to be calculated by: 


CT N t F cos0 (8) 

where N is the total number of upsets that occur while the 
ion bean is gated on. 

Using the average over layer thickness, 8 X 3 = 4.3 pm, and 
the collection depth, 8 X 4 = 6.1 |im obtained from the SRAM 
alpha particle measurements, the energies of the ions used in 
the test were computed from TRIM [7] at the entry point of the 
Si collection layer, X3 and at the exit point of the Si collection 
layer, X4. The average LET of the ion in the collection layer 
was then calculated by: 


E (X3 )-E (X4) 

p8X4 


(9) 


These results are shown in Table 2 and in Figure 10. 


Table 2. Heavy Ion Test Results. 


ION 

PLOT 

SYMBOL 

ANGLE, 0 
(degrees) 

a 

(cm2) 

LET 

(MeVcm2/mg) 

1-127 

D 

0 

232 

58.8 

Br-79 

■ 

45 

209 

57.3 

Br-79 

O 

35 

144 

49.5 

Br-79 

• 

0 

63.8 

40.3 

Fe-56 

X 

45 

33.0 

37.5 


CM 

E 

D 

Z 

o 

h 

O 

Lii 

CO 

(O 

CO 

o 

oc 

o 



LET, (MeV cm 2 /mg) 


Figure 10. Heavy ion test results performed on latch with V 0 
- V D D = 5 V. 
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ABSTRACT 

Because of the recently discovered importance of protons 
to the upset of spacebome electronics, a custom 4k-bit SRAM 
chip was tested with protons. The SRAM was developed to 
determine the Single Event Upset hardness of CMOS latches 
using alpha particle measurements, by adjusting an offset 
voltage that reduces the charge required to upset a cell. The 
proton experiments were designed to observe both proton and 
silicon recoil produced ionization. The silicon recoils were 
generated by protons undergoing nuclear coulomb scattering. 
It was discovered that silicon recoil produced charge can be 
collected from very deep in the silicon substrate. This paper 
describes a calibration procedure for the SRAM detector. 
Source spectra were acquired with this chip by measuring 
the number of upset cells versus offset voltage. 

I. INTRODUCTION 

A custom 4k-bit SRAM chip, developed to measure the 
Single Event Upset hardness of CMOS latches through alpha 
particle measurements [1], can also function as a proton 
detector. The SRAM detector has an offset voltage that adjusts 
the cell critical charge to upset, allowing a threshold to be set, 
below which its cells will upset. This threshold can be set 
below the charge deposited by protons passing through, or 
stopping in, a charge collection region. Because of the recently 
discovered importance of protons to the upset of spacebome 
electronics, as observed on the CRRES (Combined Release 
Radiation Effects Satellite) [2], the SRAM detector was tested 
with protons. 

The proton experiments were designed to observe both 
proton and silicon recoil produced ionization. The silicon 
recoils were generated by protons undergoing nuclear coulomb 
scattering. A charge collection region thickness of 6.64 Jim 
was measured for proton produced ionization, which compares 
favorably with the charge collection region thickness of 6.33 
|im for alpha particle produced ionization [1]. The charge 
collection region thickness for silicon recoil produced 
ionization was discovered to be nearly equal to the protons 
range in the silicon substrate. 


This paper describes a calibration procedure for the SRAM 
detector, allowing spectrometers to be designed for measuring 
proton, helium, and heavier ion environments inside spacecraft 
computers. The detector was calibrated for protons using the 
Caltech Tandem Van de Graaff. The SPICE circuit simulation 
program is used to compute an effective calibration curve and 
this curve is then used, with the proton data, to compute an 
effective charge collection depth, allowing calibration [1,3,4]. 


n. SRAM DETECTOR DESIGN 

The SRAM design was submitted to the MOS 
Implementation System (MOS IS) and fabricated at a 1.6-p.m 
n-well double-metal CMOS/bulk foundry. A schematic 
diagram of the SRAM cell is shown in Figure 1. The pulsed 
current source is used to model a particle strike on drain 
Dn2 when calculating the critical charge of the cell with 
SPICE. This cell differs from that of a standard six- transistor 
SRAM cell in three wavs: (1) the source of the p-MOSFET, 
Mp2, is connected to an adjustable offset voltage, V Q , instead 
of V DD to provide a control of the cell’s critical charge; (2) the 
drain area of n-MOSFET Mn2, Dn2, has been enlarged by a 
factor of four over minimum to enhance upset rates, thus 
reducing measurement time; and (3) the cell is imbalanced by 
widening Mn2 over minimum to enhance its SEU sensitivity 
versus V 0 . The bloated drain, Dn2, is fabricated in a p-type 
substrate doped to 5xl0 14 cm“3. This substrate doping value 
was determined through a dopant profiling measurement 
utilizing a large edgeless n-MOSFET. 

In operation all the memory cells are written into a 
"sensitive" state where Mn2 is turned OFF and Mp2 is turned 
ON, connecting Vq to the bloated drain, Dn2. V 0 is then 
lowered from VdD = 5 V allowing the SRAM to accumulate 
upsets at a given V 0 value. Thereafter V 0 is returned to V dd 
and the cells read to determine the number of upsets. This 
cycle is repeated at different values of V 0 . 


© IEEE 1992 
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Figure 1. Schematic diagram of the SRAM cell showing the 
placement of V 0 and the bloated n-drain, Dn2 

m. SPICE' ANALYSIS 

The critical charge, Cfc of the SRAM cell was determined 
as a function of V 0 using MOSIS supplied SPICE parameters. 
The parasitic nodal capacitances were modeled by fixed metal 
and polysilicon interconnect capacitances and by the drain 
depletion capacitances given by their areas and peripheries and 
by the SPICE level-2 junction capacitance parameters. A 
triangle current pulse with a 1:19 riseifall shape was used to 
upset these cells. For a given pulse height, the transient 
simulation was evaluated to 100 ns where the response was 
compared to Vdd/ 2 to determine if the cell had flipped. The 
current pulse height was adjusted using the binary search 
algorithm until the difference in charge (area under current 
pulse) between successive simulation runs differed by less 
than 1 fC. These results were found to be invariant with 
current pulse widths up to 500 ps. Since the proton width is 
about 200 ps [1], the response of these circuits exceeds that 
of the proton current pulse. The resulting critical charge 
versus V 0 for the V2-node of the SRAM cell is shown in 
Figure 2. This curve is the effective calibration curve used 
with the proton data to calibrate the SRAM detector. 

As V 0 decreases the critical charge of the SRAM cell 
decreases approximately linearly and goes to zero at V 0 = V om 
= 1.8 V. DC SPICE simulations of these bistable circuits 
show that the sensitive state and metastable state approach 
each other as Vo decreases until they are equal when 
V 0 = V om . For V 0 < V om the circuit spontaneously goes to 
the flipped state and is no longer bistable. V 0 m is a function of 
the threshold voltages and geometries of the n- and p- 
MOSFETs. Process induced dispersion in the MOSFET 
parameters gives rise to a dispersion, V omo =0.0l 1 V, in V om for 
the array of SRAM cells. The slope of the critical charge 
curve in Figure 2 is defined as the upset capacitance, C u = 
8Qc/5V 0 = 56 fC/V. 



OFFSET VOLTAGE V Q ( V ) 


Figure 2. SRAM critical charge calculated using SPICE 
showing error bars associated with spontaneous upset data [1], 
The SRAM calibration curve. 

IV. SRAM PROTON RESPONSE 

The response of the SRAM to protons from the Caltech 
tandem Van de Graaff is shown in Figures 3 and 4. Each data 
point was obtained by setting the offset voltage, V G , and 
allowing the SRAM to accumulate upsets until an 
independently measured fluence, F, was reached. Upon 
reaching the preset fluence the SRAM memory was cleared 
and the process repeated for k runs. The sum of the number of 
upsets for k runs was recorded by a NIM counter. 

The normalized cross section, A = A m /A a , where A m is 
the measured cross section and, A a = 68.8 jim^, is the as- 
drawn Dn2 sensitive drain area. The data points above the A 
= 1 line in Figure 3 contain a contribution from protons hitting 
outside the sensitive cell areas. Some of the charge deposited 
by these peripheral hits is collected by the Dn2 drain [1]. If 
this effect were not present, all data points would have an A ^ 
1 . The alpha particle data, shown in Figures 3 and 4, was taken 
from Figure 5 in Reference 1 and its A values were derived 
with Equation 6 in Reference 4 using the sensitive drain area, 
Aa- 

For the proton data, the average number of upset cells per 
run, <n>, in k runs is <n> = N/k, where N is the accumulated 
number of upsets recorded by a NIM counter, k = 50 for all 
data points except for the 5.0 MeV data at two points, k = 10 at 
V 0 = 2.1 and 2.2 V. <n> is related to the normalized cross 
section, A, by the equation: 



where F is the independently measured fluence and Nj = 
4096, is the total number of SRAM cells. Equation 1 is the 
homogeneous solution of Equation 2 [also Equation 2 in 
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reference 1] with the initial condition at t = 0 of no upset cells, 
ie. <n> = 0: 


= A m <KNT - <n>) (2) 

where <f> is the flux and F = $ t is the fluence. 



OFFSET VOLTAGE V (V) 
o 

Figure 3. The SRAM response to direct ionization from 0.55 
and 1.0 MeV protons and 4.7 MeV alpha particles. The mean 
offset voltage values are shown at one half the cross section. 
The measured cross section is normalized to unity at the Dn2 
drain area. 

The dispersion in F is, F a = V F/A c , where A c = 0.02 + 
0.0 - 0.0004 cm^ j s p ro ton counter cross section. The 
counting system was calibrated with an alpha source calibrated 
silicon surface barrier detector. This calibration compared 
favorably with Moliere scattering in the gold foil used to 
produce a uniform proton beam over the SRAM and the 
position of the counter relative to the SRAM. The number of 
counts equals the fluence divided by 50 + 0.0 - 1.0. Each k run 
was for 1,000 counts, except for 5.0 MeV data at offset 
voltages above 1.94 V. At 1.95, 1.96, and 1.98 V the number 
of counts was 1E4, and at 2.1 and 2.2 V the number was 1E5. 
The dispersion in <n> is given by: 


<n a > = 


km 


( m M/2 
' iNi 

i=l J 


(3) 


where m is the number of experiments each consisting of k 
runs. The number m = 1 for all data points except for the 0.55 
MeV data at two points, m = 6 at V 0 = 1.95 V and m = 2 at 
V 0 = 2.0 V. The dispersion in A is given by: 


A^ 2 _fF ^ 2 / <n CT > \ 

A) I^fJ + ^FA m (N T -<n>)J 


2 


(4) 


The data points in Figure 3, with 1 > A > 0.01, are shown in 
the Figure 4 as the percent probability that particle i will 
deposit a charge equivalent voltage V 0 i that is greater than the 
energy threshold equivalent offset voltage V 0 . The A values > 

1 cannot be plotted on probability paper and the large errors 
on the A values < 0.0 1 place little weight on these points. The 
percent probability A values are fit to a line to measure the 
mean, Vop^, which is equivalent to a probability density 
gaussian peak centroid, and the standard deviation, V 0 p<j, of 
these integral gaussian distributions. 

The 0.55 MeV mean is a measure of the proton energy 
minus the energy lost in the device overlayers. The 0.55 MeV 
proton is stopped in the charge collection region of the Dn2 
drain. This is true because the standard deviation is the sum, in 
quadrature, of the proton beam energy dispersion, the drain 
overlayer material straggling, and the SRAM spontaneous 
upset voltage, V om > standard deviation, V oma . The proton 
beam dispersion was generated by a 1.2 pm thick gold foil 
placed up stream to produce a uniform beam over the SRAM. 
The 0.55 MeV proton beam energy was measured with an 
alpha source calibrated silicon surface barrier detector at 557 ± 
11.5 keV. No contribution to the SRAM energy peak 
dispersion from straggling in silicon below the overlayer 
material is measured for 0.55 MeV protons. The 1.0 MeV 
proton data shows a significant contribution to its measured 
standard deviation from straggling in silicon below the 
overlayer material and is therefore a delta energy peak. 



Figure 4. The proton and alpha particle normal voltage 
distribution means and standard deviations. 

The upset capacitance C u , the spontaneous upset voltage 
V om p ± V 0 ma, and the hole-electron pair charge-energy 
factor for silicon K = 44.2 fC/MeV give the mean and 
standard deviation voltages, shown in Figure 4, in terms of 
energy by the following equations: 
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Eopp- £ ( v opfi “ v omjj) 


C u 

Eopa = ^opcr 


(5) 

( 6 ) 


where (Qj/K) = 1.267 MeV/V, V 0 mp = 1.80 V, and V oma = 
0.011 V. 

Applying Equations 5 and 6 to the data mean and standard 
deviation values gives, Eq^ ± Eq^ = 1 .140 ± 0. 139 MeV for 
4.7 MeV alpha particles, Egp^ ± Eq^ = 0.237 ± 0.020 MeV 
for 0.55 MeV protons, and E 0 pjx ± = 0.366 ± 0.032 MeV 

for 1.0 MeV protons. 

The Dn2 drain overlayer and charge collection region 
thicknesses, shown in Table 1, are computed from Equation 5 
and the alpha particle and proton range tables for silicon [5], 
Silicon is used for both the overlayer and charge collection 
region materials. 

Table 1. Proton and alpha particle delta range and delta energy 
in the SRAM silicon equivalent overlayer and silicon charge 
collection regions. 


Particle 
E (MeV) 

Alpha 

4.7 

Proton 

0.55 

Proton 

1.0 

Overlayer 
AE (MeV) 

0.68 

0.313 

0.191 

AR (pm) 

4.32 

4.32 

4.32 

Collection 
AE (MeV) 

1.14 

0.237 

0.366 

AR (pm) 

6.33 

— 

6.64 


Because the 0.55 MeV proton stops in the charge collection 
region, the difference between the incoming proton energy 
and the energy given by Equation 5, 0.237 MeV, equals the 
energy lost in the silicon equivalent overlayer, 0.313 MeV and 
the overlayer thickness is the difference in their respective 
ranges, 4.32 pm. After passing through 4.32 pm, the 1 .0 MeV 
protons deposit 0.366 MeV, from Equation 5, in the charge 
collection region. The difference in the ranges of the proton 
entering the collection region and exiting the region after 
losing 0.366 MeV is 6.64 pm, the effective charge collection 
region depth for protons. After passing through a 4.32-pm 
silicon equivalent overlayer, the alpha particles have an 
effective charge collection region depth of 6.33 pm. A 5.1-pm 
effective charge collection region depth for this CMOS/bulk 
process was measured with 255 MeV iron. This measurement 
required using a D-latch circuit having an offset capacitance of 
850 fC/V [1], 

The SRAM has a 4.32-pm silicon equivalent overlayer 
thickness and a 6.64-pm effective charge collection depth for 
normal incident protons. These depths are a function of the 
SPICE computed upset capacitance, C u , and would change for 
different computed values. But for the purpose of predicting 
the SRAM response to low energy protons and alpha particles 


this calibration procedure is functional. Analysis of proton 
straggling in these silicon regions demonstrates the 
functionality of the SPICE assisted calibration [4]. The energy 
straggling is given by the following equation: 

E<j = ( 47tz2e 4 nZAR ) 1/2 (7) 

where z and Z are the proton and silicon atomic numbers, n is 
the number density of silicon, and AR is the silicon thickness. 
The standard deviations given by Equation 6 is equivalent to: 

Eopa “ t (Eoma)^ + (EAua)^ + (Esio)^ 1 ^ (8) 

The spontaneous upset standard deviation is 14 keV, the 1.2- 
pm Au foil produces a 12 keV standard deviation in the 
proton beam energy, and 4.32 pm of Si produces a 9 keV 
standard deviation. It follows from Equation 8 that Eop a = 21 
keV for the 0.55 MeV data which compares favorably with the 
value of 20 keV given by Equation 6. For the 1.0 MeV proton 
beam, 11 pm of Si produces a 14 keV standard deviation. 
Applying Equation 8 gives, Eop a = 23 keV, which is 22 keV 
less, in quadrature, than the 32 keV given by Equation 6. This 
difference is equivalent to a 0.31 -pm standard deviation in the 
effective charge collection depth. 

The large tail on the 5.0 MeV data shown in Figure 5 is 
generated by Rutherford Scattering (RS). Because 5.0 MeV 
protons lose energy as they penetrate the substrate to a total 
range of 213 pm, the RS cross section must be summed over 
the effective range. The normalized cross section. A, for RS is 
given by the equation: 

R/ AR 

A = n £ (c c (Tl))i(AR)i (9) 

i=l 

where a c is the RS cross section, n is the silicon number 
density, and R is the effective RS path length of the protons in 
the SRAM substrate. a c is given by: 

BTi 

<r c (Ti) = (4n/ BTi ) J (doc/ d£2 ) dT 4 (10) 

M 

where Tj is the proton energy before scattering and (AoJ dft ) 
is the Rutherford differential scattering cross section. The 
maximum proton energy transferable to a silicon nucleus, in a 
head-on collision, is BTi and T 4 is the minimum silicon recoil 
energy required to upset an SRAM cell. B = (4Mi M2 )/(Mi + 
M 2)2 , where Mi and M2 are the proton and silicon mass 
numbers. T 4 is given by the equation: 

T4=y(V 0 -V 0p|I ) UD 
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where V 0 pp. = 1.852 volts, given by the V 0 value at A = 0.5 
for the 5.0 MeV data in Figure 5. 

The proton effective path length, R, used in Equation 9 is 
given by: 

R = R(Tmax)-R(Tmin) (12) 

where R(T max ) is the incoming proton path length in silicon 
and R(T m ui ) >s the proton path length when its energy has just 
dropped below that required to produce an upset, in a head-on 
collision. 



Figure 5. Rutherford scattering compared to the 5.0 MeV 
proton data assuming that all the silicon recoil charge is 
collected. 

The Rutherford scattering curve falls below the 5.0 MeV data, 
indicating that the sensitive node cross section is larger for RS 
than for direct ionization. The protons undergoing RS do not 
need to pass directly through the sensitive node to cause a 
flip. 



Figure 6. The maximum depth from which silicon recoil 
produced charge can upset the SRAM. 


The data indicates that silicon recoil produced charge is 
collected from very deep in the substrate. The effective 
Rutherford scattering path length, R, in Equation 9, is plotted 
as a function of offset voltage in Figure 6, and is the 
maximum depth from which silicon recoil charge can upset 
the SRAM for 5.0 MeV protons. 

V. CONCLUSIONS 

The SPICE assisted calibration utilizing 56 fC/V and the 
proton data identified a 4.32-pm silicon equivalent overlayer 
and a 6.64 ± 0.31 ftm effective charge collection depth for 
protons and a 6.33-pm collection depth for alpha particles. 
These collection depths can be used to predict the SRAM 
detector response to proton produced ionization in space. The 
SRAM cross section is larger for Rutherford scattering than for 
direct ionization. Rutherford scattering charge is collected 
from very deep in the SRAM substrate. 
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ABSTRACT 

Data from a MOSFET matrix on two JPL CRRES chips, 
each behind different amounts of shielding, is presented. Space 
damage factors are nearly identical to ground test values for 
pMOSFETs. The results from neighboring rows of MOSFETs 
show similar radiation degradation. The SRD (Space Radiation 
Dosimeter) is used to measure the total dose accumulated by 
the JPL chips. A parameter extraction algorithm that does not 
underestimate threshold voltage shifts is used. Temperature 
effects are removed from the MOSFET data. 

INTRODUCTION 

The MEP (microelectronics package) of the CRRES 
(Combined Release and Radiation Effects Satellite) included 
12 custom JPL test chips [1]. The chip layout, identifying the 
MOSFET matrix, is shown in Figure 1. This paper presents 
MOSFET matrix data on chip A3, located on the MEP outer 
board, and chip B4, located on the MEP middle board. The 
middle board is located behind the outer board which shields 
chip B4. 

The MOSFET matrix layout is shown in Figure 2. Each 
row is identical, except for the two field oxide devices in row 
one, and each pair of columns contains a different size 
transistor. The drawn transistor width and length, W/L, in 
micrometers and the channel type, n or p, are shown in Figure 
2. Neighboring rows of transistors, of the same size and in the 
same column, show nearly identical radiation degradation. 
Data is presented from orbits 500 to 700, with a major solar 
flare occurring at orbit 596 on March 23, 1991. Space 
pMOSFET damage factors nearly equal ground test values and 
their threshold voltages track SRD (Space Radiation 
Dosimeter) total radiation dose through the flare. 

The on-chip MOSFET matrix cell addressing circuitry is 
shown in Figure 3 [1]. This addressing circuitry is controlled 


by digital signals coming from the row and column decoders 
and works without effecting MOSFET measurements, even 
when considerable radiation degradation occurs [2]. In this 
experiment the pMOSFET matrix failed at 60 krad. 

These, and earlier results [3], based on radiation-induced 
threshold voltage shifts, have laid the ground work for space 
experiments utilizing pMOSFETs as flight dosimeters. On- 
chip pMOSFET dosimeters will be flown on the JPL neural 
net and SRAM chips on the United Kingdom STRV (Space 
Technology Research Vehicle). 

CRRES ENVIRONMENT 

The CRRES was launched on July 25, 1990 into a 9 hour 
50 minute 18.2° orbit with a 350 km perigee and 33,580 km 
apogee. The satellite experienced a major solar flare resulting 
in an increased dose rate to the CRRES chips after orbit 596. 



Figure 1. The JPL CRRES chip containing a MOSFET matrix, 
timing sampler, and lk-bit SRAM. 
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The CRRES SRD (Space Radiation Dosimeter) [4] data 
was used to measure the total dose seen by the CRRES chips. 
The SRD consists of four silicon surface-barrier diode charge 
particle detectors, each shielded by an aluminium dome of a 
given thickness. In this paper SRD data from dome 2 was 
taken as representative of the dose for the middle board, and 
data from dome 1, for the outer board. The data is shown in 
Figures 4, and 5, respectively. The step increase in the dose 
shown in Figure 4 at orbit 590 identifies the solar proton flare 
and the increase beginning at orbit 596 is due to electron belt 
formation induced by the flare. The electron belt is 
responsible for the large increase in total dose observed by the 
SRD and the MOSFET matrix. The total dose curve slope in 
Figure 5 identifies four dose rate regions. Only two dose rate 
regions are identified in Figure 4. 
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Figure 2. The MOSFET matrix layout showing the row and 
column placement, channel type (n or p), and transistor size 
(W/L) in jim/|im. 



Figure 3. The MOSFET matrix cell layout showing digital 
switch addressing circuitry and the MOSFETs under test, 
N-XT and P-XT. 


JPL CRRES CHIP 

The CRRES chip was fabricated using a non-radiation hard 
3-pm CMOS/bulk p-well process with a gate oxide thickness 
of 500 A (50 nm). As shown in Figure 1, each CRRES chip in 
the MEP contains: (a) a MOSFET Matrix with 32 addressable 
MOSFETs for characterizing the effects of total dose on 
MOSFET parameters [1], (b) a Timing Sampler for measuring 
radiation-induced degradation of signal propagation delay, and 


(c) a lk-bit SRAM for evaluating SEUs. This paper presents 
results from only the MOSFET matrix. 



ORBIT NUMBER 


Figure 4. CRRES MEP middle board dose as measured by 
the SRD 232.5 mil (5.906 mm) A1 dome number two. 



Figure 5. CRRES MEP outer board dose as measured by the 
SRD 82.5 mil (2.096 mm) A1 dome number one. Four dose 
rate regions are identified by the total dose curve slope. 


A total of 120 IV data points where measured per 
transistor. This extensive collection of MOSFET data allows 
the radiation effects evaluation of critical MOSFET 
parameters and represents a first attempt at acquiring 
MOSFET curve-tracer data from space. The MOSFETs were 
biased in the OFF state when not being measured, which is the 
worst case bias condition [5]. ID was measured for 10 settings 
of VD and 12 settings of VG # The measurement resolution for 
VDD is 5 V ±2%, VDS, and VGS ranges from 0 to 5 V ± 
lmV, and IDS from 2 nA to 1mA ±0.5%. If the measured 
current is less or equal to 1 nA then it is set to 1 nA [6]. In 
order to reveal spurious data, MOSFET I-V curves were 
plotted every 100 orbits in the linear, saturation, and 
subthreshold regions. This paper presents results taken from 
the saturation region for orbits 500 to 700. 
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EXPERIMENT RESULTS 


Threshold voltage damage factors, threshold voltage shifts, 
and channel mobility degradations are evaluated over orbits 
500 to 700 on MOSFET matrix transistors located on two 
CRRES chips. Chip A3(00) is located on the outer board and 
chip B4(03) is located on the middle board. Figure 6 shows a 
damage factor measurement on device 00 VP 12, which is on 
chip (00), is a threshold voltage damage factor on a pMOSFET 
(VP), and is located in row 1 at column 2 (12) in the 
MOSFET matrix. This matrix location has a W/L of 6/3 
pm/jim as shown in Figure 2. Orbits 596 to 608, shown in 
Figure 6, are defined by the low data density caused by the 
high dose rate region shown in Figure 5. The absolute value of 
threshold voltage in Figure 6 allows pMOSFET values to be 
plotted as positive numbers. 



TOTAL DOSE, D (krad) 


Figure 6. Outer board damage factor measured over orbits 500 
to 660 on MOSFET 00VP12. The data density identifies the 
four dose rate regions shown in Figure 5. The space damage 
factor, IVTpDl, is very close to the ground test value of 21.7 
mV /krad. 


The pMOSFET damage factor, IVTpDl = 21.4 mV/krad, is 
very close to the Co-60 ground test value of 21.7 mV/krad [1]. 
The zero dose threshold voltage, IVTpOl, is very close to the 
value, 0.56 V, measured during the first few orbits. A damage 
factor of 22 mV/krad was previously measured over the first 
200 orbits using the SRD to measure the dose [6]. The space 
and ground measured pMOSFET damage factors are nearly 
equivalent because the threshold voltages track the total dose 
as shown in Figure 7 for device 00VP12. Figure 8 compares 
the SRD measured dose to the 00 VP 12 MOSFET measured 
dose, using the ground measured damage factor of 21.7 
mV/krad. The MOSFET curves in Figures 7 and 8 are nearly 
identical. The value of presenting pMOSFET data, as 
dosimeter data, as shown in Figure 8, is recognized by its 
ability to explain the anonymously high damage factor shown 
in Figure 9. 
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Figure 7. The pMOSFET 00VP12 threshold voltage tracks the 
SRD dome 1 total dose through the solar flare. The pMOSFET 
threshold voltage was scaled to dose using IVTpDl =21.4 
mV/krad and IVTpOl = 0.58V shown in Figure 6. 



ORBIT NUMBER 

Figure 8. The pMOSFET 00VP12 dose and SRD dome 1 dose 
versus orbit number. The pMOSFET dose was derived from its 
threshold voltage scaled to dose via the ground test damage 
factor of 21.7 mV/krad and IVTpOl = 0.56 V. 

The pMOSFET IVTpl values plotted in Figures 7 and 8 are 
temperature corrected [6], and scaled to dose using the 
following equation : 

iVT p l = IVTpOl + IVT p dI • D (1) 

where the damage factor is given by: 

vt p d=™ (2) 

In Figure 7, D is the SRD dose. The pMOSFET threshold 
voltages were converted to dose using the SRD dose and 
IVTpDl, and IVTpOl from Figure 6. In Figure 8, D is computed 
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from Equation 1 using (a) the pMOSFET threshold voltages, 
(b) IVTpDl = 21.7 mV/krad derived from Co-60 ground test 
[1], and (c) IVTpOl = 0.56 V measured during the first few 
CRRES orbits. 

This data indicates that pMOSFETs used as on-chip 
dosimeters are capable of tracking the high dose rates induced 
by large solar flares. The outer board data is truncated at orbit 
660 due to MOSFET failure at 62 krad. The pMOSFET is a 
good dosimeter right up to its catastrophic failure dose. 

Figure 9 shows a threshold voltage damage factor 
measurement on device 03 VP22. This pMOSFET is located on 
the middle board in row 2 and column 2. The space measured 
damage factor of 25.6 mV/krad is too large because of a dose 
rate region mismatch between the SRD and the MOSFET 
shown in Figure 10. The ground test measured damage factor 
of IVTpol = 21.7 mV/krad, IVTpOl = 0.55 V measured during 
the first few CRRES orbits, and Equation 1 are used to 
compute the MOSFET dose shown in Figure 10. The high 
dose rate region between orbits 596 and 618, shown in Figure 
10, increases the damage factor fit line slope shown in Figure 
9. The mismatch is due to the low quality of the shielding 
over the middle board. The shield quality is a function of the 
uniformity of its thickness and density. The middle board is 
shielded by the outer board; whereas, the SRD is shielded by 
the high quality aluminum dome. Gaps in the outer board 
shield allowed more electrons to penetrate creating the high 
dose rate region between orbits 596 and 618 at the middle 
board, as shown in Figure 10. 



Figure 9. Middle board device 03VP22 measured over orbits 
500 to 700. The SRD dome 2, shown in Figure 4, is used for 
total dose. The damage factor is high because of low quality 
shielding over the middle board. 

MOSFET data on the middle and outer boards is used to 
compute damage factors for W/L = 6/3 jim/jim for both n and 
pMOSFETs in matrix rows 2 and 3. The damage factors are 
shown in Figure 11. The pMOSFET values were very close to 
the ground test measured value of 21.7 mV/krad. The ground 


test measured value of - 37.9 mV/krad for the nMOSFETs is 
greater than the space measured value, - 9.4 mV/krad, because 
of nMOSFET annealing. 



ORBIT NUMBER 


Figure 10. Dose measured by MOSFET 03VP22 using the 
ground test damage factor of IVTp^i =21.7 mV/krad and 
IVTpOl = 0.55 V. The MOSFET dose identifies a third dose 
rate region between orbits 596 and 618 that is not shown by 
the SRD dose. 


> 



SRD TOTAL DOSE, D (krad) 

Figure 11. Linear fits used to compute damage factors for 
MOSFETs on the middle and outer boards. The pMOSFET 
damage factors are very close to the ground test value of 21.7 
mV/krad. The MOSFET W/L = 6/3. 

The MOSFETs located in neighboring matrix rows show 
nearly identical radiation degradation. The pMOSFETs fail at 
62 krad, as shown in Figure 11, and the nMOSFETs fail at 
39.5 krad at orbit 597. The pMOSFET failure curve shown in 
Figure 1 1 above 62 krad is not included in the damage factor 
linear fit 

All of the data presented up to this point represents 
equivalent geometry MOSFETs with a W/L of 6/3. Other W/L 






values, for the middle board 03 chip, are shown in Figure 12. 
This data shows that equivalent geometry pMOSFETs 
experience similar radiation degradation by comparing 
devices from neighboring matrix rows. 
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Figure 12. Middle board data showing that equivalent 
geometry pMOSFETs experience similar radiation 
degradation. 


MOSFET ANALYSIS 


This equation is used to fit the IV data as shown in Figure 13. 



GATE VOLTAGE, VG (V) 


Figure 13. Values for 00P12 (chip A3(00), p-MOSFET, row 
1, column 2) at orbit 500. The absolute value symbols are not 
shown. 

The temperature coefficients for the threshold voltage and 
channel mobility are given by, 


A simple model was used to extract the threshold voltage, 
VT, and mobility, p,, from the saturation region of the 
MOSFET using the largest VD for both n- and p-channels. In 
the saturation region the drain current is given by, 


vtt-^ (7) 

|l(D=|l(To)(^) n (8) 


and 


3 (VG - VT)2 
2 (1+ 0 (VG - VT)) 


(3) 

where the values, VTp=-2 mV/°C for n-channels, -1 mV /°C 
for p-channels and the n-factor in Equation 8 is equal to -1.2. 

(4) These equations are used to remove the 35 orbit temperature 
periodicity, shown in Figure 14, from the MOSFET data [6], 


where 6 is the channel conduction factor, KP=|iC , ox is the 
intrinsic conduction factor, and p is the channel mobility. C' 0 x 
is the gate oxide per unit area, 0 is the mobility degradation 
factor and W and L are the as-drawn channel width and length 
respectively. The mobility degradation is caused by an 
increase in the carrier mean-free-path due to the vertical 
electrical field produced by the gate voltage [7]. This 
degradation decreases the rate of increase in the drain current 
with increasing gate voltage as shown in Figure 13. 

Applying the Taylor series approximation, 

1 1 - f VG (5) 

1 + 0VG z 

and letting 0 • VT terms vanish, a linear equation is derived 
for the square root of the drain current. 
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>/io = ;-vt + vg-|vg2) 


(6) Figure 14. Temperature measured on the outer board from 
orbit 500 to 700 showing a 35 orbit periodicity. 
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CONCLUSIONS 


This experiment demonstrated the ability to obtain curve- 
tracer type transistor data from a spacecraft environment. The 
pMOSFET threshold voltages track total radiation dose and 
their damage factors are very close to the ground test value. 
This allows pMOSFETs to be used as on-chip radiation 
dosimeters. Equivalent geometry MOSFETs located in 
neighboring matrix rows experience nearly identical radiation 
degradation. Equivalent geometry MOSFETs located on 
different chips, under different amounts of shielding, have 
nearly identical damage factors. The pMOSFETs accurately 
measure total dose right up to catastrophic matrix failure. The 
nMOSFETs show annealing effects and failed at the solar flare 
onset. 
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The RADMON (Radiation Monitor) was developed as a compact device with a 4-kbit 
SRAM particle detector and two p-FET total dose monitors. Thus it can be used as 
a spacecraft radiation alarm and in situ total dose monitor. This paper discusses 
the design and calibration of the SRAM for proton, alpha, and heavy ion detection. 
Upset rates for the RADMON, based on a newly developed space particle flux 
algorithm, are shown to vary over eight orders of magnitude. On the STRV (Space 
Technology Research Vehicle) the RADMON 's SRAM will be used to detect trapped 
protons, solar flares, and cosmic rays and to evaluate our ability to predict 
space results from ground tests. 


1. MISSION OVERVIEW 

The STRV (f, ace-Technology Research Vehicle) consists of two spacecrafts as shown in 
Figure 1. each vehicle is approximately a cube 44 cm (17 in.) on a side .'".at will be 
launched into a geosynchronous transfer orbit (GTO) with 200-km perigee, a 36,000-km 
apogee, and an 11-hour period. This orbit provides the opportunity to profile the Earth's 
radiation belts. 

The STRV contains a number of radiation measuring instruments including the JPL 
RADMON (Radiation Monitor), the DRA CREDO (Cosmic Environment Dosimeter), and the 
ESA REM (Radiation Environment Monitor). The CREDO and REM use PIN diodes or 
surface barrier (SB) detectors with pulse-height detection circuitry to 
discriminate particle energy. The area of the CREDO and REM d^ectors is between 
22 to 400 mnr; whereas, the area of the RADMON is only 0.17 mm. Thus the upset 

rates are much higher for the CREDO and REM than for the RADMON. 

The purpose of this paper is to present the methodology for calculating the 
RADMON's upset rate using a newly developed algorithm for the space environment 
and the RADMON's device parameters. After launch, the upset rate for the CREDO, 
REM and RADMON will be compared to determine the validity of the ground 
predictions and the relative sensitivity of the three particle detectors. 

The space particles relevant to RADMON upsets are: (a) trapped protons (TP), (b) 
galactic cosmic rays (GCRs), and (c) solar flare particles (SFP) . Electrons in 
the electron belt will not upset the RADMON. Electrons provide a total dose which 
causes a slow drift in the RADMON baseline which is constantly monitored. 


27 



2. SEU SRAM 


The RADMON, shown in Fig. 2, is a 28-pin chip with a 4-kbit SRAM for SEU (single- 
event upset) detection and two p-FETs [1] for total dose monitoring. The SRAM 

contains a six-transistor memory cell which is illustrated in the schematic cross 
section of Fig. 3. This figure shows that the memory cell has been modified in 
two ways to increase the cell's sensitivity to upsets. First, the drain area, <7 m , 
of diode Dn2 has been increased to increase its particle capture cross section. 
Second, the diode Dn2 is biased with an offset voltage, V Q> which is applied 

through p-FET Mp2. This allows the cell to be operated near its metastable point 
which in turn allows the cell's upset sensitivity be adjusted via V Q . 

The interaction of high-energy particles with the RADMON is described in terms of 
the various layers the particle passes through as seen in Fig. 4. In passing 
through the layers, the particle loses energy and in general increases its LET 

(linear energy transfer) or its ability to deposit charge. Thus knowing the 
thickness of the shield and device overlayer is key to calculating particle upset 
rates . 

The response of the RADMON to protons is illustrated in Fig. 5 where 0.55, 1.0, 
and 2.0 MeV protons upset the cells. The figure also includes the spontaneous 
upset curve which has a mean value of V ps „. The spontaneous curve defines the 

metastable point for each of the SRAM cells. This curve represents the baseline 
for the device. 


In operation the RADMON is operated at a V Q greater than V 0? „ or <5V 0 = V op „ - 
where V 0DU is the mean particle offset voltage defined by the a m /2 crossing point 
shown irrFig. 5. For the STRV, the RADMON is operated with three <5V 0 values 
corresponding to a proton threshold, an alpha particle threshold, and a heavy ion 
threshold. These thresholds are listed in Table 1. 


The vertical axis of Fig. 5 displays the particle cross section. It is calculated 
using: 

(1) 0 - K/W m -VN t ) 

where <p m is the measured particle flux, t„ is the measurement time, and N t is the 
number of bits in the SRAM in this case N t = 4096. This equation is applicable 
for N « N t which is satisfied when the number of flipped bits is less than 10 
percent of N t . This equation indicates that the cross section is proportional to 
the number of flipped bits, N. The proportionality factor, the denominator, 
normalizes the data with respect to flux and measurement time. 

In order to predict the RADMON's upset rate, a^, <5Q C , and <5X4 must be known. The 
device cross section is taken from the physical layout and in this case o m = 42.1 
//m . The critical charge is related to V Q through: 

(2) SQ C = C U -6M 0 

where C u is the upset capacitance. The collection depth, <5X4, is shown in Fig. 3. 
Both C u and <5X4 are determined experimentally and in this discussion C u = 28.9 
fC/V and <5X4 = 2 ^m are used. 

The RADMON's proton response, shown in Fig. 5, reveals a non-monotomc behavior 
between the particle energy and <5V Q . This behavior is explained by analyzing the 
proton charge depth curves [2] shown in Fig. 6. The vertical lines indicate the 
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width of the overlayer and the width of the collection depth; that is, 6 X 3 = 5.6 
pm and <5X 4 = 2 /urn. The figure shows that the 0.55 MeV proton stops in 6 X 4 and the 
1.0 and 2.0 MeV protons stop beyond <5X a . The 0.2 MeV proton stops in the 
overlayer so it can not be detected. The explanation for the non-monotonic energy 
versus <5V Q behavior follows from the charge deposited ir. the collection layer, 
6 Q c , as indicated in Fig. 6 . This charge is proportional to <5V Q as given in Eq. 2 
and thus explains the non-monotonic behavior. 

It is concluded that the SRAM detects protons in the 0.5 to 2.0 MeV range. Higher 
energy protons can be detected only if their energy is brought in to the detection 
range by using a shield which drops the proton energy into the 0.5 to 2.0 MeV 
energy range. The STRV contains 16 RADMONs; half are shielded by a 2.0-mm A1 

shield and the other half are shielded by a 5.5-nun A1 shield. The effect of 
shielding on the detection of protons is discussed in a later section. 


3. GALACTIC COSMIC RAY AND SOLAR FLARE PARTICLE UPSETS: 

Galactic cosmic ray (GCR) and solar flare particle (SFP) upsets are illustrated in 
Fig. 7. These curves were obtained by fitting Adam's curves [3]; the GCR are for 
solar minimum conditions. The curves show that shielding has a significant effect 
on the SFPs but relatively little effect on the GCRs. 

The number of upsets per SRAM is given by: 

(3) U = N ■ 0 • t • a m • F 

where N is the number of bits per SRAM, 0 is the solid angle, t is the observation 
or stare time, a m is the area sensitive to upsets, and F is the integral particle 
flux. The Petersen Equation [4] is commonly used to calculate SRAM upsets. The 
flux expression is: 

(4) F(l/pm 2 • ster • sec) = 4.4xlO”* 2 /[L(MeV-cm 2 /mg )] 2 

where L is the linear energy transfer. The conversion factor for L in silicon is 
L(MeV-cnr/mg) = L(fC/pm)/10.25 where L(fC/pm) = <5Q c (fC)/<5X 4 (pm) . The Petersen 
equation holds for a 0.63-mm A1 shield in a 10 percent worse case environment 
(10%WC) . 

The effect of shielding on the particle flux was obtained by fitting the Adam's 
data [3] . The flux is: 

(5) F(l/pm 2 -ster-sec) = F Q - (L/L 0 ) a+b ' L0G ( 1+D/D °) 

where D is the shielding thickness in g/cm 2 . The shielding thickness is given by 
<5X1 = D/p where p is the shield density in g/cm . The values for F 0 , L 0 , a, and b 
are listed in Table 2. 

The upsets expected for the RADMON are shown in Fig. 8 where the three LET 
thresholds are labeled. For the SFPs the RADMON will experience numerous upsets 
during the 600 sec stare time. For the GCRs, the RADMON will experience less than 
one upset per 600 sec stare time. The results are summarized in Table 1. 
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4. PROTON BELT UPSETS: 


The profile of the proton belt is shown in Fig. 9. It shows that the STRV remains 
in the vicinity of proton belt peak for several hundred seconds. Before reaching 
the RADMON the protons must pass through the A1 shields. From the proton range 
curves [2], it can be shown that the protons must have Eg > 18.5 MeV to pass 
through a 2.0-mm A1 shield and Eg > 32.5 MeV to pass througn a 5.5-mm A1 shield. 
The proton energy must be lowered by the shield so that it has an energy between 
0.5 and 2.0 MeV when it enters the RADMON 1 s collection in order to be detected. 
This effect is shown graphically in Fig. 10 where two solid horizontal lines have 
been drawn at the bounds of the RADMON’ s detection range. The problem now is to 
determine 6Eg, the incremental incident proton energy that can be detected by the 
RADMON. 

The number of protons that upset the RADMON during the stare time, t, is: 

(6) N = <V* d -t-N t -$E 0 

where <p ^ is the differential incident proton flux which is determined from the 
peak j)f the proton belt curves shown in Fig. 9. For Eq < 30 MeV, tfu = 600 
p/(cnr -sec-MeV) and for Eg > 30 MeV, = 200 p/(cnr -sec-MeV) . The <5Eg is given 
by: 

(7) <5E 0 = (dEg/dEjjSEj 

where <5E^ is the incremental exit proton energy. The dEg/dEi is determined from 
the proton range curves shown in Fig. 10 where dEg/dEj = 0.36 for both the 2.0- 
and 5.5-mm A1 shields. 

A value for 6E^ is obtained by assuming that most of the proton energy is lost in 
the shield and virtually none is lost in the chip overlayer. Thus <5Ei = <5E? = <5Eo 
= 2.0 - 0.5 = 1.5 MeV. Using Eq. 7, <5Eg = 0.54 MeV for both A1 shields. Thus the 
RADMON will detect protons through the 2.0-mm A1 shield with E 0 = 18.77 ± 0.27 MeV 

and through the 5.5-mm A1 shield with E- 0 = 32.77 ± 0.27 MeV. This shows 
that the RADMON is really an energy detector for high-energy protons. 

The number of trapped-proton induced upsets, calculated using Eq. 6, is listed in 
Table 3. It shows that several hundred upsets can be expected near the peak of 
the proton belt. The above calculation does not include particle scattering in 
the A1 shield. Including this effect will increase the number of upsets. 


5. DISCUSSION: 

The number of upsets estimated for the RADMON varies over eight orders of 
magnitude and depends on the type of the space environment encountered. The upset 
estimates depend on the accuracy of the environmental models and device parameters 
such as the collection depth. A simple algorithm was developed for estimating 
space particle flux that includes the effect of shielding. The flight of the 
RADMON on the STRV along with the CREDO and REM will provide data to quantify the 
sensitivity of the three particle detectors and validate ground predictions. 
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Table 1. RADMON Particle Thresholds and LETs. 


THRESHOLD 

VDD 

V 

Vos/u 

V 

sv o 

V 

LET o 

MeV-cirr/mg 

PROTON 

3 

1 

0.61 

0.086 

ALPHA 

3 

1.139 

0.250 

0.353 

HEAVY 

5 

1 


3.280 

4.628 


Table 2. GCR and LSF Flux Parameters 


PARAMETER 

UNITS 

GCR 

LSF 

F 0 

l/(^m^-ster-sec) 

lxlO'5 

lxlO‘2 

n° 

D o 

a 

MeV-cm^/mg 

lxlO -4 

lxlO -J 

mg/cnr 

78.36 

0.5990 

unitl ess 

-1.92 

-1.425 

b 

unitl ess 

-1.35 

-0.974 


Table 3. RADMON Anticipated Upset Rates 


ENVIRONMENT 

PROTON BELT 
U(upsets/SRAM) 

SOLAR FLARE 
U(upsets/SRAM) 

10% WORSE CASE 
U(upsets/SRAM) 

GCR 

U(upsets/SRAM) 

A1 -SHIELD 

2.0-mm 5.5-mm 

2.0-mm 5.5-mm 

0 . 63-mm 

2.0-mm 5.5-mm 

m 

3.3E+2 1.1E+2 

1.7E+3 5.4E+2 
1.5E+2 3.5E+1 
2.0E+0 2.3E-1 

1.93E-1 

1.15E-2 

6.68E-5 

7.5E-3 7.2E-3 
5.0E-4 4.7E-4 
3.6E-6 3.3E-6 
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Figure 1. STRV orbit relative to the earth's proton belt. 



Figure 2. STRV-RADMON chip 2.57 x 2.70 mm 2 . 
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Figure 7. Integral LET i nterpl anetary spectra for various aluminum shields [3] 
due to galactic cosmic rays (GCR) , 10 percent worse case environment (10%WC), and 
solar flare particles (SFP) . 
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Figure 8. STRV RADMON SEU SRAM upset rate due to galactic cosmic rays (GCR), 10 
percent worse case environment (10%WC), and solar flare particles (SFP). 
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Figure 9. STRV traverse of the earth's proton belt. 
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Figure 10. Proton exit energies for various aluminum thickness and incident 
proton energies. For X A i = 2 mm, E 0 = 20 MeV and for X A1 =5.5 mm, E n = 35 MeV. 
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ABSTRACT 

On-chip p-FETs were developed to monitor the 
radiation dose of n-well CMOS ICs by monitoring 
threshold voltage shifts due to radiation induced 
oxide and interface charge. The design employs 
closed geometry FETs and a zero-biased n-well to 
eliminate leakage currents. The FETs are 
operated using a constant current chosen to 
greatly reduce the FET's temperature sensitivity. 
The dose sensitivity of these p-FETs is about -2.6 
mV/krad(Si) and the off-chip instrumentation 
resolves about 400 rad(Si)/bit. When operated 
with a current at the temperature-independent 
point, it was discovered that the pre-irradiation 
output voltage is about -1.5 V which depends only 
on design-independent silicon material 
parameters. The temperature sensitivity is less 
than 63 pV/°C over a 70°C temperature range 
centered about the temperature insensitive point. 

1. INTRODUCTION 

The use of FETs (Field-Effect Transistors) as 
dosimeters was pioneered by Holmes-Siedle [1]. 
A number of these devices have flown on Earth- 
orbiting satellites [2-4], 

In recent years, p-FET dosimeters have been 
developed with specially-grown thick-gate oxides 
which have a large number of oxide traps. 
Sensitivities of >10 mV/rad(Si) [5] have been 
achieved. The sensitivity to radiation can be 
enhanced by applying a large positive bias during 
radiation which forces more of the positive oxide 
charge to the interface. The sensitivity to 
temperature can be minimized by operating the p- 
FET with a current at the temperature- 
independent point [6, 7], 

In this work, a p-FET dosimeter is developed 
under the constraint that the dosimeter be useful 
in predicting the radiation dose of an 1C fabricated 
with a non-radiation hardened 1.2-pm CMOS 
process. As shown in Fig. 1, two p-FETs were 
fabricated on the RADMON (Radiation Monitor), 
which also includes an SEU-SRAM for monitoring 
particle upsets. This appears to be the first time a 


p-FET dosimeter was fabricated in a fine-line, 
thin-oxide semiconductor technology. 

The p-FET is biased to about -1.5 V during 
measurement and is unbiased when not being 
measured. This approach to biasing is intended 
to provide a known bias environment at all times. 
In certain applications, the availability of 
spacecraft power is unpredictable. Thus, being 
unbiased during irradiation provides a known bias 
scenario. 

On-chip dosimetry provides the advantage that 
the dose is measured directly next to the 1C. This 
reduces the uncertainty inherent in dosimetry 
calculations which are complicated especially for 
highly shielded electronics. 



Figure 1. STRV-RADMON chip 2.6 mm x 2.7 mm. 


2. p-FET DESIGN 

The RADMON, shown in Fig. 1, contains two p- 
FETs. The geometry of p-FET4 is W/L = 182/4 
pm/pm and p-FET8 is W/L = 182/8 pm/pm where 
W and L are the channel width and length, 


©IEEE 1993 


37 









respectively. The layout of the p-FET, shown in 
Fig. 2, features a closed geometry design where 
the drain completely surrounds the source. The 
closed geometry eliminates the bird's beak or 
channel-edge leakage encountered in a linear 
FET. 


The schematic cross section of the device, shown 
in Fig. 3, indicates that the n-well and source are 
separated so they can operate at slightly different 
biases required by the operational amplifier, U1. 
This allows all the forced current to flow through 
the p-FET channel. Drain-to-well leakage is 
shunted to ground. Grounding the n-well is a 

departure from normal CMOS circuit operation 
where the n-well is normally connected to VDD. 



Figure 2. p-FET, MP4, layout where the gate 
length is L = 4 pm and width is W = 182 pm. 



In operation, all terminals of the p-FET are 
operated near ground except the drain which 
operates near VO = -1.5 V. During irradiation, the 
device is biased in the off state. The two RL 
resistors are used to bleed off any charge 
remaining on the p-FET. The instrumentation is 
calibrated by positioning the switch to CAL shown 
in Fig. 3. This replaces the p-FET with calibration 
resistor RC. 

The p-FET constant drain current, ID, is 
established by VR and Rl or ID = VR/RI. The 
output of amplifier U1 is scaled by U2 so that it 
makes maximum use of the range of the 8-bit ADC 
(Analog-to-Digital Converter). As seen in Fig. 3, 
the gate is connected to the drain. This insures 
that the p-FET is operated in saturation. 

3. p-FET MODEL EQUATIONS 

In saturation the p-FET drain current [7] is: 

3 (-VO + VT) 2 

(1) ID= 

2 1 + G(-VO + VT) 

where VO is the p-FET output voltage, 15 = 
KP W e /L ei KP = p C 0 and VT is the p-FET 
threshold voltage. W e = W - DW and L e = L - DL 
are the effective channel width and length, 
respectively, p is the zero-field channel mobility, 
C 0 is the gate oxide capacitance per unit area, 
and 0 is the mobility electric-field degradation 
parameter. 

The above equation was simplified by taking its 
square root. Then the term was linearized using 
Taylor Series expansion: 

(2) NiF = ^'IftFj-VO + VT)[1 - ^ (-VO + VT)] 

2 

Expressions for the temperature and dose 
dependence of VT, (5, and 0 are given below. The 
equations are expanded about the reference 
temperature T 0 . In these equations, the 
temperature and dose effects are assumed to be 
independent. That is, the equations do not 
include a dose temperature product term. 

The threshold voltage is described by: 

(3) VT = VT 0 + VT t (T - T 0 ) + VTq D 

where D is the dose, VT 0 is the threshold voltage 
at T 0 and D = 0, VTj = dVT/3T| T -> To and VT D = 

5VT/3D|q->o- The temperature [7] and dose 
dependence [8] of 15 is given by: 
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(4) (S = IJ 0 (T/T 0 )-" + (J D D 

where G 0 is (1 evaluated at T 0 and D = 0, Gq = 
3G/3D|q g, n characterizes the mobility temperature 
dependence, and T is the absolute temperature. For 
D T = 3D/3T = 0: 

(5) Gy = 3G/3T = -n G/T 

The temperature and dose dependence of 0 is given by: 

(6) 0 = 6 0 + 0j (T - T 0 ) + 0 D D 

where 0 O is 0 evaluated at T 0 and D = 0, 0j = 
30/3T|y jo and 0 q = 30/3D|q q. 

4. p-FET TEMPERATURE EQUATIONS 

(INCLUDING 0) 

In this section the p-FET equations are derived 
including the 0 parameter. This allows an accurate 
analysis at the operating or measurement 
temperature, T m . This temperature is usually 
different from the reference temperature, T 0 . The 
measurement temperature, T m , might be the mean 
spacecraft operating temperature. 

The p-FET output voltage follows from Eq. 2. The 
solution requires solving the equation which is 
quadratic in VO: 

(7) VO = VT - 1(1 - Nl - 0 N[8 ID/G]) 

0 

The sign of the square root is negative which can 
be verified by evaluating this equation in the limit 
0 -> 0. The value of VO at T m is calculated from: 

(7a) VO m = VT m --1(1 - Nl - 0 m N[8 ID m /G m )) 

0 

The current, ID m , at the measurement 
temperature is found from the above equation by 
setting 3VO/3T|j j m =0 and solving the 
resulting quadraticTequation for ID m : 

(8) NlD m = + A . c ± E-N[1 + B (-C + A/4)] 

2 

where 

(9) A = a 2 b/d 2 

(10) B = b 

(11) C = c/d 

(12) E = a/d 

(13) a = VTj + 0j/0 m 2 

(14) b = -0 m '(8/G m ) 

(15) c = 0 T /9 m 2 

(16) d = -[(20 T /0 m ) - (n/T m )]/V(2G m ) 

The sign of the square root in Eq. 8 is positive for 
0 > 0 and negative for 0 < 0. 


5. p-FET TEMPERATURE EQUATIONS (0 = 0) 

The equations derived in Section 4, include the 0 
parameter. Since 0 is small, it is neglected in this 
section to gain physical insight into the meaning 
of the equations. 

The p-FET IV characteristics are plotted in Fig. 4 
using Eqs. 1, 3, and 4 for D = 0 and the 
parameters [7] listed in the Fig. 4. This figure 
shows that the so called "temperature 
independent" point is in fact ill-defined when 
viewed in detail. 

In this analysis, the most important p-FET 
parameter is VO. Its temperature dependence is 
analyzed by expressing Eq. 2 as follows: 

(17) VO = VT - 2ID/G 

The VO is calculated for D = 0 using Eqs. 3 and 4 
and plotted in Fig. 5. These curves show that 
there is a point at which VO is independent of 
temperature. This point is determined by 
differentiating Eq. 17 with respect to temperature 
and setting the result to zero. This leads to the 
simple expression for the measurement current: 

(18) ID m = 2G m 3(-VT T /G Tm ) 2 

The value for ID m = 19.2 pA was calculated using 
the parameters shown in Fig. 5. The five curves 
were plotted with ID values that vary by one 
percent. This shows the effect of missing the 
target current of ID m . As seen in the figure, the 
peak in the curve moves to higher temperatures 
as ID increases. The value of VO changes by 
0.25 percent for a one percent change in ID. 

The temperature sensitivity of VO can be 
expressed simply by combining Eqs. 4, 5, 17, 
and 18 which leads to: 

(19) VO T = 3VO/3T = VTj[1 - (T m /T)^' n/2 ] 

This shows that VOy = 0 at T “ T m for any n or 
VOy = 0 for n = 2 for all T. 

The effect of temperature variations on VOy is 
analyzed with the help of the horizontal line shown 
in Fig. 5. This line is 1 mV below the target value 
of VO m = -1.5 V. The value of 1 mV was chosen 
because it is the same magnitude as the dose 
effects to be measured. The VOy = -63 pV/°C at 
point "a" and 63 pV/°C at point "c" seen in Fig. 5. 
The temperature range between "a" and "c" is 
more than 70°C. This means that VOy is less than 
±63 pV/°C over the 70°C temperature range 
centered at the measurement temperature, T m . 
This is a significant improvement in the 
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uncompensated temperature sensitivity of VTy = 2 
mV/°C used in calculating Figs. 4-6. 

The solution for n = 2 means that the curves 
shown in Fig. 5 have zero slope. For this case, 
the p-FET IV characteristics display a true 
temperature insensitive point as shown in Fig. 6 . 



Figure 4. Expanded p-FET IV characteristics 
showing that the "temperature independent" point 
is not a point. Note n * 2. 



Figure 5. The temperature and current 
dependence of the p-FET output voltage using Eq. 
17. 


Next, ID m is analyzed for design purposes by 
substituting IS m = KP m W e /L e into Eq. 18; this 
leads to: 

(20) ID m = 2KP m (VT T T m /n)2w e /Le 

This equation shows that ID m depends on the 
silicon parameters, KP m , VTy, and n, the 
measurement temperature, T m , and the FET 
geometry. Thus, once the silicon parameters and 


T m are known, the designer is free to choose ID m 
by adjusting W and/or L. 

Finally, VO is analyzed for design purposes by 
substituting Eqs. 5 and 18 into Eq. 17: 

(21) VO m = VT m - 2VTy T m /n 

This equation shows that VO m depends on the 
silicon parameters, VT 0 , VTy, and n and the 
measurement temperature, T m . Note that VO m is 
independent of FET W/L geometry. Thus the 
designer has no control over this value. In this 
study, the value for VO m is about -1.5 V. With 
irradiation, VO m becomes more negative. The 
initial VO m value is important from an 
instrumentation stand point. Since the ADC 
shown in Fig. 3 has an input voltage range 
between 0 to 5 V, the pre-radiation value for VZ in 
Fig. 3 is set to zero volts by compensating VO m 
by setting RO = -RA VR/VO m . 



Figure 6 . Expanded p-FET IV characteristics 
showing a true temperature independent point 
for n= 2 . 


6 . p-FET TEMPERATURE DATA ANALYSIS 

The p-FETs were measured in packages in an 
oven using an hp4062 parametric test system. 
The temperatures are estimated to be accurate to 
within ±1°C. The measurements were obtained by 
forcing VO = 2 V and measuring ID 5 . Then four 
additional currents were forced at v/l Di = 
O. 2 VID 5 , 7ID2 = O. 4 VID 5 , VID 3 = O. 6 VID 5 , and 
V/ID 4 = O. 8 VID 5 . Once these values are 
determined, they are used throughout the rest of 
the measurements. 

The experimental data, shown in Fig. 7, was fitted 
using the method of least squares. In the 
analysis, the following parabolic equation was 
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used. The coefficients of the equation were 
related to the parameters in Eq. 2: 

(22) \/|D = a 0 + a-! VO + a 2 V0 2 
where 

(23) a 0 = 7(13/2) VT (1 - 0 VT/2) 

(24) a-\ = -v/(U/2) (1 - 9 VT) 

(25) a 2 = -v/((i/2 )• 6/2. 

The "a" parameters are used to obtain VT, 13, and 
0 for each IV curve. The VT solution was 
obtained by setting VO = VT at ID = 0. This leads 
to a quadratic equation whose solution is: 


(26) VT = a 1 /(2a 2 ) (-1 +V^1-4a 0 a 2 /a 1 A 2) 

The sign of the square root is positive. This can 
be verified by setting 0 = 0 or a 2 = 0. The 
solution for li is: 

(27) II = 2(a 1 + 2a 2 VT) 2 
The solution for 0 is: 

(28) 9 = 2a 2 /(a-| + 2a 2 VT) 

The IV points, shown in Fig. 7, were fitted using 
the above procedure. The li, VT, and 0 values are 
listed in Table 1 for three temperatures. 



Figure 7. 1.2-pm CMOS p-FET IV temperature 

response for flight chip number W12P4C26. 


The temperature parameters for VT and 0 were 
extracted by least squares fitting the data listed in 
Table 1 by using Eqs. 3 and 6. The temperature 
parameter for li was extracted using Eq. 4 after it 
was linearized by taking the logarithm. The 
temperature parameters are listed in Table 2. 


Values for ID m and VO m were calculated using 
Eqs. 7a and 8 and the parameters listed in Table 
2. The results are listed in Table 2 and plotted as 
the vertical and horizontal lines in Fig. 7. 


Table 1. Flight p-FET4 Parameters (W12P4C26) 


T 

°C 

VT 

V 

a 

mA/V 2 

0 

1/V 

30 

-0.874 

1.100 

0.056 

75 

-0.800 

0.891 

0.041 

125 

-0.703 

0.708 

0.027 


Table 2. Flight p-FET4 Parameters 
(T 0 = 300 K, T m = 10°C, W12P4C26) 


PARAM 

UNITS 

MEAN STDEV 

*0 

nm 

21.8* 

W 

um 

182 

L 

um 

4 

DL 

um 

0.3459* 

VT 0 

V 

-0. 882110.0051 

vt t 

mV/°C 

1.8002+0.0800 

&o 

mA/V 2 

1.124010.0110 

KP 0 

pA/V 2 

24.704110.2423 

n 

— 

1.614010.0528 

^To 

(pA/V 2 )/°C 

-6.0471 


1/V 

0. 057110.0013 

0T 

(1/kV)/°C 

-0.315310. 0210 

VT m 

V 

-0.9125 

mA/V 2 

1.2339 

®m 

I^Tm 

1/V 

0.0624 

(pA/V 2 )/°C 

-6.6384 

**-*m 

pA 

244.2 

vo m 

V 

-1.554 


* MOSIS supplied parameter. 


7. p-FET DOSE DATA ANALYSIS 

The p-FET dose dependence was determined 
using Cobalt-60 irradiation. The devices were 
irradiated with their lids on, at room temperature, 
at 1 rad/sec, and at zero bias. The p-FETs were 
measured at room temperature within 15 minutes 
after the Cobalt-60 irradiation. The p-FET were 
also annealed at room temperature for times 
measured after the final irradiation. 

The p-FET irradiation results, shown in Fig. 8, 
were fitted using Eq. 22. This produced a set of 
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VT, IS, and 0 values for each dose value. These 
values are plotted in Figs. 9 to 1 1 for four p-FETs. 
The radiation results are listed in Table 3 for one 
of the p-FETs. 

The VT values, plotted in Fig. 9, show a high 
degree of linearity during irradiation and a slight 
recovery with annealing. The group average 
slope of the VT vs dose curve during irradiation is 
VT d = -1.698±0.038 mV/krad(Si). The shift in VT 
with radiation is consistent with the build up of 
positive oxide charge and interface states. The 
slight recovery of VT during room temperature 
annealing is consistent with the slight loss of 
oxide charge. The interface state density is 
stable as seen by the flat response of the mobility 
during annealing. 

The li values were converted to zero-field hole 
channel mobility using: 

(29) p = (L - AL) IS/(W C 0 ) 

In this analysis, the value for C 0 was calculated 
using the MOSIS supplied value for the gate-oxide 
thickness, X Q = 21.8 nm. 

Table 3. Ground Test p-FET Cobalt-60 Radiation 
Parameters for Chip Number W12P4C05 


PAR. 

UNITS 

MEAN STDEV 

vt d 

mV/krad(Si) 

-1.674±0.004 

(iD 

pA/V 2 /krad(Si) 

-3.585±0.164 

KPd 

pA/V 2 /krad(Si) 

-0.079±0.004 

0D 

1/kV/krad(Si) 

-1.184±0.018 

V0Dm 

mV/krad(Si) 

-2.591±0.018 



Figure 8. 1.2-pm CMOS p-FET IV dose/anneal 

response for chip number W12P4C05. 



W1HJ7IUIC DOSE/ANNEAL, D(krad(Si))/A(hr) 

Figure 9. Four 1.2-pm CMOS p-FET threshold 
voltage dose/anneal responses. 



PM8L171IJac DOSE/ANNEAL, D(krad(Si))/A(hr) 

Figure 10. Four 1.2-pm CMOS p-FET hole 
mobility dose/anneal responses. 



Figure 11. Four 1.2-pm CMOS p-FET theta 
dose/anneal responses. 
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The data, shown in Figure 10, are clustered 
according to the p-FET's W/L ratios. This ratio is 
designated by the "P" number given in the wafer 
number listed in the figure. That is, for P4 the 
ratio is W/L = 182/4 pm/pm and for P8 the ratio is 
W/L = 182/8 pm/pm. This dependence on the W/L 
ratio does not affect the output voltage, for VO is 
geometry independent as seen in Eq. 21. 

As seen in the Fig. 10, p is stable after annealing. 
The build up of positively charged interface states 
during irradiation degrades the mobility due to the 
increase in channel scattering centers [9]. Since 
the mobility does not anneal at room temperature, 
the interface state density is stable. 

The 0 values, shown in Fig. 11, are clustered 
according to the p-FET W/L ratios. The 
characteristics show no annealing, which 
indicates that the interface state density is stable 
during annealing. 

The 0 parameter represents the curvature term for 
the parabolic equation given in Eq. 2. The 
curvature is difficult to see in Fig. 8 because 0 is 
very small. As seen in Fig. 11, the maximum 
value for 0 is 0.06 1/V. The influence of 0 in Eq. 2 
is less than 3 percent because it enters as e/2. 
This small value for 0 justifies the use of the 
Taylor Series expansion in simplifying Eq. 2. 

As shown in Fig. 11, the charge is positive at low 
dose and negative after high dose and annealing. 
The sign of 0 indicates the direction of the ID vs 
VO parabola given by Eq. 2. This point is 

illustrated in Fig. 12. For 0 > 0 the parabola 
points down and for 0 < 0 the parabola points up. 
For 0 = 0, the \/lD vs VO relationship is a straight 
line. 



WJ.I916.KIC OUTPUT VOLTAGE. VO (V) 


Figure 12. Plot of Eq. 2 for different 0 values. 


8 DOSIMETRY 

The device dose sensitivity is derived in this 
section. Since the dose is measured at the 
constant current, ID m , the radiation dose 
sensitivity of VO is greater than VTq. This is 
evident in Fig. 8 where the spread in the curves is 
wider at ID = ID m than at ID = 0 due to the dose 
dependence of the mobility. 

In this section the output-voltage dose sensitivity, 
VOq, for 0 = 0 is calculated by differentiating Eq. 
17 with respect to dose at ID = ID m . Then ID m , 
given in Eq. 18, is substituted into the result. The 
output-voltage dose sensitivity is: 

(30) VO Dm = aVO/0D| Tm = VT d - VT T IJ D /li Tm 

A value for VOpm = -2.59 mV/krad(Si) was 
calculated using the values listed in Tables 2 and 
3. This result is considerably larger than VTq = - 
1 .67 mV/krad(Si). 

The output voltage for the four p-FET samples is 
shown in Fig. 13. This plot was obtained from 
data sets like those shown in Fig. 8 where the VO 
values were obtained at ID m . The data shows a 
nearly linear rise in VO with dose and a slight 
annealing effect. 



Figure 13. Four 1.2-pm CMOS p-FET output 
voltage dose/anneal responses determined at the 
current, ID m . 

The similarity of the data is remarkable 
considering the p-FET samples came from 
different wafers from the same run. This means 
that the sample set can be assumed to be uniform 
and that the results can be used to calibrate the 
flight parts which, of course, can not be irradiated 
on the ground. 

For the circuitry shown in Fig. 3, which uses an 8- 
bit ADC to span 100 krad(Si), the resolution is 
100k/256 = 390 rad(Si)/bit. This means that the 
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p-FET dosimeter can easily resolve a dose of 1 
krad(Si). 

9. CONCLUSION 

The use of on-chip p-FET dosimeters has been 
established and provides a radiation sensitivity of 
-2.6 mV/krad(Si) for the 1.2-pm CMOS used in 
this study. The temperature dependence is less 
than ±63 pV/°C over a temperature range of 70°C 
centered about the p-FET temperature 
independent point. At this point, the pre- 
irradiation output voltage was determined to be 
about -1.5 V and this value is design independent. 
The use of on-chip p-FETs provides a direct 
measure of the radiation dose experienced by the 
associated CMOS 1C. 
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ABSTRACT 

An inverter matrix test circuit was designed for 
the Clementine space mission and is built into the 
RRELAX (Radiation and Reliability Assurance 
Experiment). The objective is to develop a circuit 
that will allow the evaluation of the CMOS FETs 
using a lean data set in the noisy spacecraft 
environment. As will be shown, only nine data 
points are needed to acquire ten CMOS FET 

parameters . 


INTRODUCTION 

An inverter matrix was designed for the RRELAX and 
included on the Clementine space mission. Two 
RRELAX units were fabricated. One unit is included 
on the Clementine spacecraft which will orbit the 
moon and fly-by the asteroid Geographos. The 
second unit is included on the Inter-Stage Adaptor 
(ISA) which goes into a translunar orbit. The 
Clementine is scheduled for launch on January 25, 
1994. 

The RRELAX, shown in Fig. 1, is by 4" x 4” x 1.5" 
(10 x 10 x 3.8 cm), weighs 22 oz (0.62 kg) and requires 
2.5 watts peak power. 
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The unit contains an 80C86 microprocessor operating 
at 2.4 MHz and data acquisition circuitry to 
measure the experimental devices. The data is 

stored in the RRELAX memory and the results sent to 
the spacecraft through a serial port for down- 
linking. 

The labels in Fig. 1 point to the experimental 
devices including SRAMs for Cosmic Ray detection, 
p-FETs for total dose measurement, and a CCD for 
radiation effects analysis. The locations of the 
two inverter matrices are shown in the figure. One 
matrix has a 7-mil (0.18 mm) A1 equivalent shield and the 
other a 37-mil (0.94 mm) A1 equivalent shield. [0.18 mm 
and 0.94 mm] 

The inverter matrix was developed to evaluate CMOS 
FETs fabricated at a rad-soft rapid-prototype CMOS 
foundry in order to evaluate its space worthiness. 
The FETs included in the inverter matrix are 
typical of FETs found in integrated circuits. In 
addition the matrix was designed to require a small 
data set in order to minimize spacecraft down-link 
capacity. As will be shown, only nine data points 
are needed to characterize ten n- and p-FET 
parameters. 

The layout of the inverter matrix, shown in Fig. 2, 
is similar to the inverter matrix developed in the 


RO TIE 



INV23A22 PCX 


REL23731 PCX 

Figure 1. RRELAX showing the arrangement of test Figure 2. Inverter matrix with 64 test inverters 

devices. The inverter matrices are labeled INVT11 arranged in four identical quadrants. The chip 

(Z08) and INVT12 (Z07). size is 1.63 x 2.58 mm 2 . 
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middle 80's for evaluating process statistics [1] 
and [2]. The matrix consists of 64 inverters 
arranged in four identical quadrants. In two of 
the quadrants, the gates are biased high during 
radiation and in the other two quadrants, the gates 
are biased low during irradiation. The dimensions 
of the inverter matrix FETs are listed in Table 1. 


The key to this approach is the dependence of the 
inverter threshold voltage, VT-j , on geometry [3]: 


( 1 ) 


VT i = 


VDD + VT n VB^ - VT p 
1 + 


where VT n is the n-FET threshold voltage, and VTp 
is the magnitude of the p-FET threshold voltage. 
The dependence of VT^ on geometry, shown in Fig. 3, 
indicates how VT^ depends on the geometry of the 
inverter, B r : 

B n KP n (W n - AW n )(L p - AL p ) 

(2) B r — 

8 p KP p (W p - AW p )(L n - Al n ) 


(7) B p = KP p (Wp - AW p )/(Lp - AL p ) 
and 

(3) 0 p = 0Lp/(Lp - ALp) 

The above FET models were chosen after numerous 
alternative models were tried. It was found that 

giving VT a spatial description led to unstable 
solutions. This is understandable since the data 
points, as seen in Fig. 5, are not located near 
the VT extraction points. Also 0 is described by 

an L* dependence only. When the offset term. Oo. 
was introduced into the model, 0 parameters were 

too unstable and fluctuated about zero. When 0W 
was introduced into the model, these values also 

fluctuated about zero, the net result is a simple 
model that fits the data reasonably well and 
extracts reasonable parameters. Also the 
extraction procedure is robust. 


DATA FITTING PROCESS 


where KP = /i 0 *C 0 ; fi 0 is the zero-field channel 

mobility and C 0 is the gate capacitance per unit 

area. The FET design width is W, and L is the 
design length. The actual width and length are 
described using AW and AL. Notice that for B r * 0, 
VTj * VDD - VTp and for B r + VT 1 * VT n . 

The transfer curves for the 16 inverters listed in 
Table 1 are shown in Fig. 4. The inverter 

threshold voltage is located at the intersection of 
transfer curve and the positive sloping Vout-Vin 

line. 


FET MODEL 


The above equations were combined into the 
following equation set and fitted to the data shown 
in Fig. 5. The FET geometries for the cardinal 
points shown in Fig. 5 are given in Table 3. The 
equations for the n-FETs are: 


(VTi-VT n )2 

(9) Vl n « a 0 + a 1 *VT i + a 2 

(L n -Al n ) 

where 


(id 

( 12 ) 

(13) 


2-ID n (L -AL n )/(W n -AW n ) 

* VT n *»K p n 

VKP " , , 

-0L n -VKP n /2 


The analysis uses the fact that at the inverter 
threshold both n- and p-FETs are in saturation. 
Thus-the following simple square-law drain-current 
expression adequately describes the n-FET behavior 

[4]s 


(3) ID n 
where 


Bn(VT i - VT n ) 2 
2 [1 + 0 n (VTi - VT n )] 


Using the above "a" coefficients, the n-FET 
parameters are: 

(14) VT n * -a 0 / ai 

(15) KP n = ai 2 

(16) 0L n - -2a 2 /a! 

The above equations are plotted as curves #4, #5, 
and #6 in Fig. 5. 


(4) B n = KP n (W n - AW n )/(L n - AL n ) 


The equations for the p-FETs are: 


and 

(5) 0 n = 0L n /(L n - AL n ) 

The p-FET equations are: 

B p (VDD - VT, - VT p ) 2 

(6) I D p « ? * 

2[1 + 0 p (VDD - V^ - VTp)] 

where 


(17) Vl p = 
where 


( VDD-VT i -VT 0 ) 2 

bo + bi‘VTi+b 2 - 

(Lp-DLp) 


(18) 

(19) 

( 20 ) 
( 21 ) 


2*1 Do ( Lp-DL p ) / (Wp-OWp) 
VTn'VKP n 


h P = 

D 0 “ * 'n p 

b! = -ViCP n H 


b 2 = 


0LpVKP p /2 
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Using the above "b" coefficients, the p-FET 
parameters are: 

(22) VT p - VDD + ao/ai 

(23) KP p = a L z 

(24) 0Lj = 2a 2 /a! 

The p-FET curves are plotted as curves #1, #2, and 
#3 in Fig* 5. 

The solution for the FET parameters, VT KP, and 0L 
were obtained by using a least squares fitting 

procedure to determine the "a" and "b" coefficients 
of the linear equations, Eqs. 9 and 17. The least 
squares procedure was embedded inside a Simplex 
solver. The Solver changed DL, DW, and VT in order 
to optimize the least squares coefficient of 
determination with the constraint that the 
estimated VT value equals the least squares VT 
value within a precision of 0.01. 

The fit of the data to the cardinal points is shown 
by the curves given in Fig. 5.^ The fit is 
reasonable considering the simplicity of the FET 
model given in Eqs. 9 and 17. By design, the 
curves converge to a single value of VT n and VTp. 
The matrix contained additional FET geometries and 
so these were included in the analysis as "bonus" 
points which help stabilize the extraction of VT n . 
As will be shown, VT p has considerable variability 
because the data points are far from the VTp 
extraction point. 

The test results for the data plotted in Fig. 5 is 
listed in Table 2. The listing indicates that 5 p- 
FET and 5 n-FET parameters have been extracted from 
the data. The RRELAX extracted results are 

compared to results published by MOSIS. The RRELAX 
results are surprisingly close to the MOSIS 
results. The discrepancies can be explained by 
differences in extraction procedures and wafer-to- 
wafer parameter variations. The MOSIS parameter 
extraction occurs in the FET linear region; 
whereas, RRELAX parameter extraction occurs in the 
FET saturation region. 


CIRCUIT DESCRIPTION 

The inverter matrix, shown in Fig. 6, consists of 
an 8-by-8 array of test inverters fabricated in a 
1.2-um n-well CMOS process. The chip has enables, 
El and E2‘, that allow two or more chips to be 
bussed together; power, VDD ( 5 V ) and ground, GND 
for decoder logic; row address, RO - R2 and column 
address CO - C2 for selecting the test inverter; 
test inverter power, HI, ground, LO, input, IN, and 
output, OUT; and logic input, TIE, that tristates 
IN and connects OUT to the addressed test inverter 
input. The 64 test inverters are arranged in 4 
identical quadrants. The quadrants contain 16 test 
inverters having unique n- and p-FET geometries 
listed in Table I. When the chip is deselected (El 


and/or E2‘ inactive), the inverter inputs are 
biased low (0 V) in QUAD#! and QUADI3 and high (5 
V) in QUAD#2 and QUAD#4. 

The inverter matrix circuit diagram in Fig. 6 shows 
one test inverter (MN1 and MP1) and supporting 
select transistors. The dashed lines indicate 

connections to the remaining seven rows and seven 
columns of the matrix. Row and column decoders are 
used to select one of the 64 test inverters for 
measurement. These decoders are deselected (no 
selected outputs) when one or both chip enables are 
inactive, putting all test inverters in a known 
biased condition. Two of the quadrants have the 
test inverter inputs I* connected to D as shown in 
Fig. 6 while the other two quadrants have inputs I* 
connected to U. 

Fig. 7 is an equivalent circuit of an enabled 
inverter matrix and the support electronics used to 
obtain measurements. The RRELAX inverter matrix 
measurements use a 12-bit Analog-to-Digi tal 
Converter (ADC). The inverter output voltage is 
measured when the output switch is in the V 
position and the current through the inverter (IDD) 
is measured when the output switch is in the I 
position. OP-amps 0P1 and 0P2 function as a 

current-to-vol tage converter which holds the LO 

point at a virtual ground. All current 
measurements are corrected by a baseline leakage 
measurement made on the matrix when it is 

deselected. When the matrix is deselected, the 

matrix baseline leakage current is the sum of 
leakage currents through 32 off n-FETs and 32 off 
p-FETs. 

When TIE is low the Digi tal -to-Analog (DAC) drives 
the inverter input and allows the measurement of 
inverter transfer characteristics. The data 
analyzed in this paper was all obtained with TIE 
high, where the internal feedback connection is 
made, bringing the circuit to the inverter threshold 
condition . 


DISCUSSION 

Representative ground test results from the 

inverter matrix, shown in Figs. 8 to 10, were 

obtained from the Clementine spacecraft. Test 
results, shown in Figs. 11 to 13, were obtained 

from the Inter-Stage Adapter (ISA). The first 

results, shown in these figures, we^ obtained at 

JPL (Pasadena, CA) using an HP 4062 parametric 
tester to measure inverter matrix packages. 
Subsequent results were obtained from the inverter 
matrices mounted in the RRELAX. The second results 
were also obtained at JPL where the RRELAX was in a 
noise free environment. The third and fourth 
results were obtained with one RRELAX unit 

integrated into the Clementine spacecraft and the 
other RRELAX unit integrated into the Inter-Stage 
Adapter. The third result came from NRL (Naval 
Research Laboratory, Washington DC); and the fourth 
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result came from VAB (Vandenberg Air Force Base, 
Lompoc, CA). 

The symbols used in Figs 8 to 13 are explained in 
Table 4. The symbols indicate (a) the bias on the 
inverter gates during irradiation and (b) the 
amount of shielding. Thus test results are 

expected to diverge into four trends once the 
devices receive sufficient radiation. 

The results shown in Fig. 8 to 13 indicate that the 
laboratory results have least variability. On 
occasion results from RRELAX can be noisy. The 
source of the noise is under investigation. 

The measurement has the interesting attribute that 
it requires only a single input voltage namely VDD 
and the output is an easily measured voltage that 
ranges between the VTn and VDD - VTp. This makes 
the measurement extremely simple and robust. These 
are the kind of attributes needed in the noisy 
environment of a spacecraft. 


CONCLUSION 

The inverter matrix developed for the RRELAX meets 
the design goal of providing FET parameters from a 
lean data set. The quality of the data could be 
improved by designing the inverter matrix FETs with 
geometries closer to the VT extraction points. 
This is important because it would reduce the 
uncertainties in the data extraction process and 
provide parameters with less variations. 

The noise encountered in the RRELAX is due to the 
spacecraft operating environment. The noise enters 
through the power supplied to the RRELAX. Also the 
RRELAX experimental devices are lightly shielded to 
allow a maximun exposure to radiation. This allows 
RF from spacecraft operations to be picked up by 
the RRELAX electronics. Hopefully, these noise 
sources will be small enough so that radiation 
effects will be clearly observed in the data. 
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Table 1. Inverter matrix inverter geometries 


DIMENSIONS (//m) 

Wp Lp Wn Ln 

7.2 

1.2 

9.6 

1.2 

28.8 

1.2 

2.4 

2.4 

28.8 

1.2 

2.4 

1.2 

28.8 

1.2 

9.6 

1.2 

7.2 

2.4 

2.4 

1.2 

7.2 

2.4 

9.6 

1.2 

7.2 

1.2 

2.4 

2.4 

7.2 

1.2 

2.4 

1.2 

2.4 

2.4 

12.0 

1.2 

14.4 

1.2 

2.4 

2.4 

12.0 

1.2 

2.4 

2.4 

7.2 

2.4 

2.4 

2.4 

2.4 

3.3 

2.4 

1.5 

1.8 

2.4 

2.4 

1.2 

19.2 

2.4 

19.2 

2.4 

2.4 

2.4 

14.4 

1.2 


Table 2. Inverter matrix test results (CZ071vabl) 
for the data shown in Fig. 5 and from MOSIS. 


PARAM 

UNITS 

— RRELAX — 
n-FET p-FET 

— MOSIS 

n-FET p-FET 

DL 

/Ml 

0.470 

0.686 

0.227 

0.116 

DW 

fjm 

0.616 

0.768 

0.524 

0.620 

VT 

V 

0.690 

0.952 

0.728 

0.903 

0L 

pm/L 

0.136 

0.118 

— 

... 

KP 

pA/V 2 

65.324 

20.160 

86.400 

26.600 

CC 

-- 

0.998 

0.995 

— 

— 
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Table 3. FET Dimensions for the curves shown in 
Fig. 5. 


CURVE 

W(/im) 


#i 

7.2 

2.4 

#2 

7.2 

1.2 

#3 

28.8 

1.2 

#4 

2.4 

2.4 

#5 

2.4 

1.2 

#6 

9.6 

1.2 


Table 4. Symbols for Figs. 8 to 13. 



SHIELD 

GATE-HIGH-BIAS 

GATE-LOW-BIAS 

DEVICE 

mil A1 

QUAD#1 QUAD#3 

QUAD#2 QUADI4 

Z07-INVT12 

37 

SQUARE 

PLUS (+) 

Z08-INVT11 

7 

DIAMOND 

CROSS (x) 



0r™^r>/0p 

Figure 3. Spatial dependence of the inverter 
threshold voltage. 



INPUT VOLTAGE, Vm (V) 


Figure 4. Inverter matrix transfer curves for one 
of the four quadrants. 



CL14141 12JCLS INVERTER THRESHOLD VOLTAGE, VTi (V) 


Figure 5. Inverter matrix data analysis showing 
the nine cardinal points (squares) and six bonus 
points (triangles). (CZ071vabl) 


| TO OTHER ROWS ] 



Figure 6. Inverter matrix circuit diagram showing 
one test inverter (MN1 and MP1) and signal steering 
transistors. 


VDD - 5V 



Figure 7. Test inverter measurement circuitry. 
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SEU/SRAM as a Process Monitor 

Brent R. Blaes and Martin G. Buehler, Member . IEEE 


Abstract — The SEU/SRAM is a 4-kbit static random access 
memory (SRAM) designed to detect single-event upsets (SEU’s) 
produced by high-energy particles, either in space or on the 
ground. This device was used to determine the distribution of 
the memory cell spontaneous flip potential. The variance in 
this potential was determined to be due to the variation in the 
threshold voltage of one of the n-FET's found in the memory cell. 
For a 1.2-//m CMOS process, the standard deviation of the n-FET 
VT was found to be 8 mV. This structure provides substantial 
statistical data; 4096 data points are measured. Using cumulative 
distribution and residual plots, outlier cells, such as stuck cells 
and nonnormally distributed cells, are easily identified. 

I. Introduction 

T HE use of matrixed test structures has been shown to be 
an effective approach to collecting statistical data with 
respect to inverter threshold voltages [1], metal steps [2], 
linewidths [2], and contact resistances [3]. Such structures 
require digital-to-analog converters (DAC) to stimulate the 
structure and analog-to-digital converters (ADC) to measure 
the required values. In this paper, the SEU/SRAM is used to 
obtain analog information about the cells without the use of 
an ADC. In this case, a DAC is used to force voltages and the 
on-chip latches (memory cells) are used to perform the ADC 
function. This allows rapid measurement of the SEU/SRAM. 

The structure used in this study is the radiation monitor 
(RADMON), shown in Fig. 1. The structure consists of a 
4-kbit SRAM for the detection of particles that cause single- 
event upsets in the memory cells and two p-FET’s used to 
detect total dose radiation [4]. This structure was fabricated 
in a 1.2- /mi CMOS process which is a scaled version of a 
RADMON fabricated using a 1.6-/mi CMOS process [5]. 

The size of the RADMON, as shown in Fig. 1, is 7.0 mm 2 . 
The memory cell layout is shown, in Fig. 2. The dimensions of 
the MOSFETs in the cell are listed in Table I. The SRAM was 
fabricated with a 1.2-/zm n-well CMOS process at a MOSIS 
brokered foundry. 

The SRAM cell schematic, shown in Fig. 3, has a six- 
transistor memory cell with an offset voltage V 0 that is used to 
evaluate the spontaneous cell flip potential. The two internal 
nodes of the cell are labeled Vl and V2. These nodes will 
be used subsequently to describe the transfer curves for the 
inverter. 

The timing diagram for the operation of the cell is shown 
in Fig. 4. This diagram shows that the memory cell has three 
states of operation: Read, Write, and Stare. In the Read and 
Write cycles, V„ = 5 V. Initially, all the cells are written into 
an initial state which is described in Fig. 3. Then the Stare 

Manuscript received October 12. 1993; revised April 18, 1994. This 
research was supported by the Strategic Defense Initiative Organization. 

The authors are with Jet Propulsion Laboratory, Pasadena, CA 91109. 
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Fig. i. KAUMUN cmp. - o nun x 1.7 mm. was iaoncaied in 1.2-/«mCMOS 



Fig. 2. Memory cell: 33.6 //m x 36.0 //m. The n-well and n- and p-select 
layers are omitted. 


TABLE I 

Dimensions and Drive Capability of SEU/SRAM MOSFET’s 
Used in the SRAM Memory Cell Shown in Fig. 3 


FET 

L (/mi) 

IF (/an) 

Ad {/nil 2 ) 

.JlfiA/V*) 

Mnl 

1.5 

2.4 

11.16 

1 10.4 

Mn2 

1.2 

2.4 

42.12 

138.0 

Mpl 

3.3 

2.4 

8.64 

16.7 

Mp2 

2.4 

1.8 

6.84 

17.3 

Mtl 

1.2 

1.8 

— 

103.5 

Mt2 

1.2 

1.8 

— 

103.5 


cycle voltages are applied to the cells. In this cycle V 0 is 
gradually lowered to a potential well above the spontaneous 
flip point. In this state, ionizing particles that deposit sufficient 
charge flip individual cells. 

If V 0 is lowered sufficiently, the memory cells will flip 
spontaneously. This is the behavior that is analyzed in this 
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Fig. 3. SRAM memory cell circuit biased into the Stare state. 


READ CYCLE fW»0. Vo = VD0«5V) 

s r~ 


WRITE CYCLE (P*0. 0«MiZ. Vo*VDD-5V) 
S ! 


STARE CYCLE (VDD 

s — 


5V) 


5V 




/' Vo (stqrg) 


Fig. 4. SRAM memory cell timing diagram. 


paper. It will be shown that the spontaneous flip potential is a 
measure of the uniformity of the threshold voltage of inverter 
1, VTi i as seen in Fig. 3. As such, the SEU/SRAM provides 
data that is useful for process control and lot acceptance. 

II. Memory Cell Model 

The cell spontaneous flip behavior is explained by the 
SRAM transfer curves shown in Fig. 5. These curves were 
generated using a simple model for the MOSFET drain current 
which does not include channel length modulation [6]. The 
inverter's input voltage is V- n and its output voltage is Vo Ut . 
These voltages are related to the memory cell voltages: that 
is, V in is V2 for IN VI and is VI for INV2. Likewise, V out is 
VI for IN VI and V2 for INV2. 

The CMOS inverter transfer curve is divided into five 
regions [7], These five regions are described in Table II with 
respect to VT n , the n-FET threshold voltage, VT P the p-FET 
threshold voltage, and VT\ the inverter threshold voltage. The 
inverter threshold voltage appears at the output of an inverter 
when its output is connected to its input. 

The MOSFET drain currents used in calculating the curves 
in Fig. 5 are listed below [7]. In the following expressions I n [[ n 
and 7pii n describe the n- and p-drain currents in the linear 
region of operation and / nsat and / psat describe the n- and 
p-drain current in the saturation region of operation. 

/„ .in = Pn(V ln ~ VT n - V' out /2)V out (1) 

/nsat =(/W2)(V- in - VTn) 2 (2) 

/pun =0 P (VDD - V in - VT P - (VDD - V OM )/2) 

■ ( VDD - V oat ) (3) 



Fig. 5. SRAM transfer curves showing the variation in the stable point, 
indicated by circles, as the offset voltage Y 0 is lowered to the spontaneous 
Hip point. 


TABLE II 

t in Regions and MOSFET Currents Used to 
Calculate the Transfer Curves Shown in Fig. 5 


Region 

V™ 

Current 

I Mill 

I 

o < v; n < rr„ 

In = / P = 0 

Ldd 

11 

vr„ < r in < it, 

Ais&i — fplin 

— 

III 

Ln — Linl — I T\ 

7m»ai “ Ip * ai 

IT, 

IV 

VTi < Ln < \'DD- VTp 

^nlin = / p*a» 

— 

V 

VDD-VT P < \ i n < VDD 

o 

II 

Q. 

II 

C 

0 


TABLE III 

MOSFET Model Parameters (Run N26D. T t> 

= 2i)°C) 

Parameter 

Units 

Mean 

Standard 

Deviation 

n-FET Results 

VT 0 

V 

O.G9± 

0.0101 

KP 0 

l> A/V- 

G9.00± 

1.2000 

AW 

//III 

0.4G± 

0.0200 

AL 

//111 

0.4G± 

0.0116 

p-FET Results 

VTo 

V 

0.95± 

0.0087 

KPo 

/'•A/V 1 

23.00± 

0.5200 

AW 

//111 

0.30± 

0.0310 

AL 

//111 

0.33± 

0.0180 


/psat - ((3 P /2)(VDD - V; n - VT P ) 2 (4) 

where 3 = KP • W e /L e and VT is the threshold voltage. For 
p-FET' s VT is the magnitude of the threshold voltage. Also 
KP — /zC ox « where fi is the channel mobility, and C ox is 
the gate oxide capacitance/area. Finally W e — W - AW and 
L e = L- AL where W and L are the as-drawn channel width 
and length, respectively, AW and AL are the channel width 
and length correction factors, respectively. The MOSFET 
values used in the following analysis are shown in Table III. 

III. Spontaneous Cell Flip Behavior 

The key to the measurements presented here is the sponta- 
neous flip behavior described by the transfer curves shown in 
Fig. 5. As seen in the figure, the memory cell has two stable 
states and one metastable state. One stable state is inscribed 
by a triangle and is located at VI = 5 V and V2 = 0. The 
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1.70 1.72 1-74 1.76 

OFrSET VOLTAGE. V 0 (V) 

Fig. 6. SEU/SRAM spontaneous flip response from 1 1 chips taken from 
1.2-^m CMOS wafer 1. 

other stable state is identified by the circles. The metastable 
states are inscribed by squares. A stable state is a bias state 
to which the inverter will return if momentarily biased away 
from by noise. If an inverter is debiased from a metastable 
state, it will not return to the metastable state but will proceed 
to a stable state. 

The initial state of the cell for V 0 = 5 V, is located at the 
upper left-hand comer of the chart. As V 0 decreases from 5 
V, the stable point follows a path described by the circles and 
the metastable point follows a path described by the squares. 
When the circles meet the squares, the cell flips to the stable 
point described by the triangle. This occurs when V 0 = 1.5 V 
and is termed the spontaneous flip voltage, V oa . At this point 
the cell has only one stable point described by the triangle. 

In essence, the transfer curve of INV2 is used to probe the 
transfer curve of INV1. In this process the metastable point 
for the cell is determined. 

The above process takes place by operating the SRAM 
through the Write, Stare, and Read cycles. That is, with V Q = 5 
V, the cells are written into the circle state. The V a is carefully 
lowered somewhat and returned to 5 V. Then the cells are Read 
to see how many have flipped. Then the process is repeated; 
that is, V 0 is lowered to an even lower bias and returned to 5 
V. Again the cells are Read and the number of flipped cells 
noted. This process is repeated until all the cells have flipped. 

In terms of measurement technology, the external stimulus 
is a DAC that provides the voltage, V Q . The ADC action is 
provided by each cell which serves as a 1-bit comparitor. This 
means that the ADC action is taking place on-chip within 
each cell. 

IV. Experimental Results 

The memory cell V r 0 distributions are shown in Fig. 6 for 
11 chips, from Wafer 1. As seen in Fig. 6, the SRAM's have 
a distribution of offset voltages at which the cells flip. As 
mentioned above, the data is acquired by lowering the offset 
voltage, V Q . and counting the number of flipped cells at that 
V Q value. The memory is then reset and the offset voltage 
lowered to a new V Q that is 1 mV lower than the previous 
value and again the number of flipped cells determined. This 
process is repeated until all 4096 cells flip. Note that these 



Fig. 7. Cumulative distribution plots for the chips shown in Fig. 6. 

curves are completely deterministic within the limits of noise. 
That is for a given V 0 . the same cells flip. The noise limit for 
these measurements was about 5 cells. 

The distributions shown in Fig. 6 are Gaussian in nature and 
can be characterized by the normal distribution with a mean 
of and a standard deviation of V Q9a * The cumulative 
distribution plots for the chips shown in Fig. 6 are shown in 
Fig. 7. This data is characterized by the cumulative probability 
function using: 

P(V ol > Vo) = 100 • (N - 0.5)/.Vf (5) 

where N is the number of flipped cells at V 0 and for this 
memory N t — 4096. The analytical formula that describes the 
cumulative distribution is 

N = N t {l - erf [(V 0 - V <HL )/V oa y/2]}/2 (6) 

where erf is the error function. The result of a least squares fit 
to each of the curves shown in Fig. 7 is listed in Table IV. The 
entire range of data was fitted. Notice that the curves with the 
largest standard deviations, namely, chips 2 and 4. have cells 
that deviate significantly from the main distribution. For chip 
4, the problem is a stuck column and for chip 2, the problem 
is a stuck bit. 

The standard deviations for these devices is the tightest 
observed to date being about 8 mV. Previously observed 
standard deviations for a 1 .6-/nn CMOS process were about 
10 mV [5], 

Selected chips are examined in detail in Figs. 8 to 10 
where the cumulative and residual distributions are shown. 
The cumulative distribution allows a critical examination of 
the tails of the distribution. The residual distribution allows 
a critical examination of the deviation of the response from 
the fitted behavior. An example of a stuck column is shown in 
Fig. 8 and a stuck cell is shown in Fig. 9. Acceptable behavior 
is shown in Fig. 10. 

The scale for the vertical axes in Fig. 7 is expressed in terms 
of probabilities. This is equal to the scales used to describe the 
vertical axes shown in Figs. 8 to 10 which are expressed in 
terms of sigmas. The number of sigmas is (V 0 - V 0fX )/V 0(r . 
The relation between the number of sigmas and probability is 
listed in Table IV [8]. 


53 






1.69 UO 1.71 1.72 1.73 1.74 1.75 1.76 
W1P8CXXA.XLS OFFSET VOLTAGE. Vo (V) 


Fig. 8. Results from chip 4 (Fig. 6) showing a stuck column. The squares 
are the cumulative probability and the triangles are the residual. 



1.69 1.70 1.71 1.72 1.73 1.74 1.75 1.76 

W1P8CXXA.XLS OFFSET VOLTAGE, Vo (V) 

Fig. 9. Results from chip 2 (Fig. 6) showing a stuck bit. The squares are the 
cumulative probability and the triangles are the residuals. 



1.69 170 1.71 1.72 1.73 1.74 1.75 1.76 

W1P8CXXA.XLS OFFSET VOLTAGE. Vo (V) 

Fig. 10. Results from chip 1 (Fig. 6) showing acceptable behavior. The 
squares are the cumulative probability and the triangles are the residuals. 

The results were simulated using a normally distributed 
sample data set. As seen in Fig. 11, several data points fall 
slightly below the fitted line in the tails of the distribution. 
This same behavior was also observed in the experimental 
data shown in Fig. 10 and so is considered normal for the 
tails of the distribution. 

A summary of the results from ail the wafers included in 
this study are shown in Fig. 12. These samples came from 
four wafers and are tightly clustered. The mean offset voltage 
has a span of 30 mV and the standard deviation of the offset 



1.69 1.70 1.71 1.72 1.73 1.74 1.75 1.76 


OFFSET voltage. v 0 (V) 

CM/SJCI* 

Fig. 11. Ideal cumulative distribution plot for \ = 1 . T2 V and \' 0 <r =8 

mV. 



1.69 1.70 1.71 1.72 1.73 1.74 

OFFSET VOLTAGE MEAN, V (V) 

*o*jc7J>*e au 

Fig, 12. Standard deviation vs. mean offset voltage for chips fabricated on 
four wafers. The lines are linear regression fits to data from each wafer. 


TABLE TV 

Cumulative Distribution Probabilities 
Related to the Number of Sigmas [8] 


Number of Sigmas 

Cumulative Probability 

4 

0.99997 

3 

0.99865 

2 

0.97724 

1 

0.84135 

0 

0.50000 

-1 

0.15865 

— 2 

0.02276 

-3 

0.00135 

-4 

0.00003 


voltage varies from 7 to 9 mV. This is considered excellent 
cell distributions for cells located 1) within a chip, 2) between 
chips, and 3) between wafers. 


V. Data Analysis 

The interpretation of results follows from observing the 
nature of the transfer curves shown in Fig. 5. A close exami- 
nation of this figure reveals that the spontaneous flip point is 
determined when V Q reaches the threshold voltage of inverter 
1, VT\\. The CMOS inverter threshold voltage is determined 
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TABLE V 

SEU/SRAM \\, Results 



Fig. 13. Inverter threshold voltage dependence on the MOSFET geometry 
factor J r . 


by the conditions given in region III (see Table II): 


VDD + VT^M-^VTp 

1 + \/0r 


(7) 


where VT n is the n-MOSFET threshold voltage, and VT P is 
the magnitude of the p-MOSFET threshold voltage. The Beta 
factor is 


0n = KP n (W n - AVF n )(L p - AL P ) 
0 P KP p (W p - AlTp)(L n - AL n ) 


The inverter threshold equation is plotted in Fig. 13 and shows 
that for 0 r — 0. VT\ = VDD — VT p and for j3 r — ► oc, VT\ = 
VT n . 

Using propagation of error analysis, the variance of the 
inverter threshold voltage is 


(1 + s/JTr) 2 (l + VFr) 2 

PrjVDD - VT n - VT P ) 2 
4(1 + \//^r) 4 


(9) 


where 


W 2 IT~ /2 

— *n<7 , y pg , 

^ Uf 2 ^ W2 f [2 

KV ep L en 


( 10 ) 


These equations show that for j3 r = 0, VT\ = (VDD - 
VT pm ) ± VTpjj and for 0 r oc, VTj = VT nfi ± VT na . 

Using (7) and (8), the spontaneous flip voltage was calcu- 
lated to be 1.48 V. This was determined by first calculating 
0 ri — 10.7 and 3 r2 = 15.1. Then VT\i = 1.48 V and 
VT [2 = 1.39 V were calculated and the assignment was made 
that V Q3 = VTii = 1.48 V. This value compares with the 
experimentally observed average value for V 0 of 1.72 V as 
displayed by the data shown in Fig. 12. 

The discrepancy between 1.48 V and 1.72 V is explained 
by the simplistic MOSFET model given in (l)-(4). If channel 
length modulation was included in the MOSFET model, then 
the transfer curve for INV 1 would have a finite slope at the 
midpoint of its transfer curve. This effect will increase the 
modeled spontaneous flip point from 1.48 V and bring the 
result closer to the experimental average value of 1.72 V. 


Chip 

WAFER 1 


r,„. 


1 

1.7246 

0.0076 

T 

1.7202 

0.0091 

3 

1.7226 

0.0081 

4 

1.7319 

0.0092 

5 

1.7209 

0.0078 

6 

1.7224 

0.0075 

7 

1.7247 

0.0077 

8 

1.7099 

0.0089 

9 

1.7236 

0.0084 

10 

1.7108 

0.0079 

11 

1.7205 

0.0078 


The standard deviation of V oa . V osrT . is due to the variation 
in the n-FET threshold voltage. This can be verified by 
introducing the FET model parameters, listed in Table III, into 
(9) and (10). This leads to 

VTl = 0.055 • VT 2 a + 0.586 • VT 2 a + 0.008G (11) 

for inverter 1. A close examination of this equation reveals 
that the coefficients are small except for the VT na term. Thus 
it is the dominant parameter. 

Now the results in Fig. 12 and Table V can be summarized 
as follows. The mean offset voltage is determined by the 
inverter threshold voltage of inverter 1. that is, 

~ V Tn ( 12) 

and the standard deviation of the offset voltage is determined 
by the variation in the threshold voltage of the n-FET in 
inverter 1, that is, 


Vos* ~ VT na \. (13) 

The result, given in (12) and (13), is a function of the layout of 
INV1. Thus by changing the layout of the cell, other features 
of the cell can be sensed. In this case, 0 r \ = 10.7 is relatively 
large which emphasizes VT n i. If 0 r \ were made small, then 
VT pl could be monitored. 


VI. Conclusion 

The SEU/SRAM provides data on the uniformity of a 
CMOS process. The 4-kbit SRAM memory cell offset voltages 
were found to be normally distributed. The offset voltage de- 
pends on the threshold voltage of inverter 1 and its distribution 
depends on the variation in VT n \ . Cumulative distribution 
plots reveal SRAM’s with stuck bits and other anomalies 
which appear in the tails of the distribution. 

The observed threshold voltage variances were about 8 ± 1 
mV for all “good” chips. Since this variance applies to within a 
chip, chip-to-chip, and wafer-to-wafer values, it is considered 
excellent threshold-voltage control. This result provides a 
measure of excellence to be met by future CMOS foundry 
runs. 
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ABSTRACT 

The Clementine RRELAX radiation monitor chip 
consists of a p-FET total dose monitor and a 4- 
kbit SRAM particle spectrometer. Eight of 
these chips were included in the RRELAX and 
used to detect the passage of the Clementine 
(S/C) and the innerstage adapter (ISA) through 
the earth's radiation belts and the 21-Feb 1994 
solar flare. This is the first space flight 
for this 1.2-jim rad-soft custom CMOS radiation 
monitor. This paper emphasizes results from 
the SRAM particle detector which showed that it 

(a) has a detection range of five orders of 
magnitude relative to the 21-Feb solar flare, 

(b) is not affected by electrons, and (c) 


detected microflares occurring with a 26. 5 -day 
period. 

1. INTRODUCTION 

The Clementine spacecraft was launched January 
25, 1994 (day 25) from Vandenberg Air Force 
Base on a Titan IIG booster. The spacecraft, 
shown in Fig. 1. contains a number of advanced 
light weight cameras used to photograph the 
moon during its 2.4-month lunar orbit. The S/C 
and ISA also contain a RRELAX (Radiation and 
Reliability Assurance Experiment) which was 
located on the S/C and ISA as seen in Fig. 1. 


ISA s/c 

RRELAX RRELAX 


SOLAR PANEL 



BOOSTER CLEMENTINE 


Figure 1. Launch configuration for the Clementine spacecraft (S/C) and the inner stage adapter (ISA) 
with a RRELAX box mounted on each. The S/C RRELAX is actually located on the opposite panel. 
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The total dose profile for the S/C is shown in 
Fig. 2. These measurements were obtained from 
p-FETs which were measured with a constant 
current set to the temperature -insensitive 
point [1]. The figure shows the passage of the 
S/C through the earth's electron belt followed 
by moon orbit which began on day 50. Shortly 
after beginning moon orbit, a solar flare was 
encountered on day 52. After leaving the moon 
on day 123 the S/C again passed through the 
earth's radiation belt and shortly thereafter 
ceased operation. 


The total dose profile for the ISA is shown in 
Fig. 3. The ISA was in a 2.1-day earth orbit 
from Feb 4 (day 35) until it re-entered the 
atmosphere on about May 8 (day 128). During 
this time it experienced essentially a linear 
increase in total dose as shown in Fig. 3. The 
solar flare was not visible in the ISA p-FET 
data because the ISA is more heavily shielded 
than the S/C. 


2. RRELAX 



00 S 2 

clm@4219.xls TIME, t (days from Jan 1,1994) 


Figure 2. Overview of the S/C p-FET dose 
profile for lightly shielded device Z2. 



Figure 3 Overview of the ISA p-FET dose 
profile for lightly shielded device Z2. 


The RRELAX is a 624 g. 2.4 watt microprocessor 
controlled experiment, shown in Fig. 4, that 
contains 166 test devices designed to evaluate 
the effects of radiation and the thermal 
environment on the performance of these 
devices. This paper discusses results from six 
RADMONs (Radiation Monitors). Z2. Z3. Z5. Z4. 
XI. and Y1 mounted on X-. Y-, and Z- faces of 
the RRELAX. 


CFETHl 


RADMON(Z2) 

'NVril CCD2 

I INVT12 


CFET81 

I 

RAOMON(Z3) 

RADMON(2S) 

RADMON(Z4) 



RAD MON (X2) 
RAD MON (XI) 


RADMON(Y1) 
RADMON{Y2) 


REL34719 PCX 

Figure 4. RRELAX components arranged in 
10.2 cm x 10.2 cm x 3.8 cm box. 


3, SRAM PARTICLE SPECTROMETER 

The RADMON chip, shown in Fig. 5, was 
fabricated in 1.2-pm radiation-soft CMOS. The 
chip is a full custom designed chip that 
contains seven components. In this paper only 
the results from the p-FET and the 4 -kbit SRAM 
are presented. The p-FET was used to determine 
the total radiation dose and the SRAM was used 
to determine the particle flux. 
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The SRAM as a particle detector enhances the 
single -event -upset (SEU) effect that is 
detrimental to SRAMs when used as memory 
elements [2], The SRAM has a six-transistor 
cell with an offset voltage, V 0 , to adjust the 

sensitivity of the cell to particle upsets 
induced by for example protons [3] or alpha 
particles [4], In addition, the sensitive 
drain of the n-FET was bloated to A 0 = 42.1 pm 2 

so the 4-kbit SRAM active area is 0.0017 cm 2 - 



MOS34424 PCX 


Figure 5. RADMON chip 2.6 x 3.4 mm 2 fabricated 
in 1.2-/jm CMOS contains both a p-FET total dose 
monitor and an SRAM particle detector. 


from the range-energy curves [5], These 
calculations require the overlayer thickness, 
DX3, and the collection layer thickness DX4. 
Values for these thicknesses are shown in the 
figure. The deposited energy is expressed as 
the difference between the energy at the start 
of the collection layer, E3, and the energy at 
the end of the collection layer, E4. During 
the rising portion of the curve, the ions stop 
in the collection layer and so E4 - 0. 
Finally, the offset voltages are stated 
relative to the spontaneous offset voltage V os ^ 

= 1.148 V. 

The SRAM was operated with two threshold 
voltages which were 0.15 V and 1.0 V above the 
spontaneous upset point. V os ^. As seen in Fig. 
6 for DV 0 = 0.15 V, protons in the 0.55 to 0.8 

MeV range can upset the SRAM and al pha 
particles with energies above 1.8 MeV can upset 
the SRAM. When operated with this threshold, 
the SRAM is said to be sensitive to particles 
with atomic numbers of Z>1. For DV 0 - 1.0 V 

only heavy ions such as oxygen or iron can 
upset the SRAM. 



CR0S4727 XLS 


ION ENERGY, E2 (MeV) 


Figure 6. Calibration curves and operational 
thresholds for the SRAM particle detector. 


The calibration of the cell, shown in Fig. 6, 
used both protons and alpha particles. The 
experimental data points were obtained at 
various incident ion energies, E2, for an 
offset voltage that corresponds to a cross 
section of Aq/ 2. The cross section is the 

upset rate divided by the particle flux times 
the number of cells which is 4096 in this case. 
The solid lines were fitted to the experimental 
data using the deposited ion energy derived 


In operation [6], the SRAM is written with ones 
with V 0 = VDD = 5 V. Then V 0 is lowered to a 

value DV 0 above V os ^ for a stare period of 100 
s. Then V 0 is returned to VDD and the number 

of zeros or particle flipped cells read. In 
addition to lowering V 0 , the VDD for the chip 

was also lowered to 3 V to reduce the chip 
power consumption. During satellite operation. 
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each SRAM was operated for 100 s with DV 0 * 
0.15 V and then for 100 s with DV 0 - 1 V. This 

cycle was continued for 1 hour. The number of 
upsets reported is the sum of the upsets 
recorded during the 18 x 100 s - 30 minutes 
that the SRAM spent in each state. 

The RADMONs were shielded by different lids as 
listed in Table 1. The devices were packaged 

Table 1. Shields in mm and Normal Incident 
Proton Energy in MeV. 


■* + 


DEV 

NO. 

S/C 

S/C 

S/C 

ISA 

ISA 

ISA 

NO. 

LIDS 

mm 

mil 

MeV 

inn 

mil 

Mev 

Z2 

0 

0.025 

1 

1.3 

0.178 

7 

4.6 

Z3 

1 

0.787 

31 

11.4 

0.940 

37 

12.7 

Z5 

2 

1.549 

61 

16.9 

1.702 

67 

17.8 

Z4 

3 

2.311 

91 

21.1 

2.464 

97 

21.9 

XI 

1 

0.787 

31 

11.4 

0.940 

37 

12.7 

Y1 

1 

0.787 

31 

11.4 

0.940 

37 

12.7 


■* + 


with an integral number of 0.25 -urn Kovar lids. 
A conversion factor of three was used to 
convert Kovar thickness to A1 thickness. That 
is 1 mm of Kovar is equivalent to 3 mm of A1 . 
The 0-lid device was covered with 13-pm Al- 
coated Kapton and the RRELAX box was similarly 
covered with this Kapton. Thus the 0-lid 
device is listed as having a 25 -pm Al shield. 
The equivalent thickness of each shield in mm 
and mils of Al is listed in Table 1. These 
numbers represent the thickness of the lid. box 
covering, and, for the ISA RRELAX, 0.152 -mm 
thermal blanket. 

The mean proton energy needed to penetrate the 
shield and upset the SRAM when operated with 
DV 0 - 0.15 V is also listed in Table 1. The 

values shown in the table were obtained from 
proton range-energy curves [5]. The shields 
slow the protons so their energy is reduced to 
between 0.55 and 0.8 MeV when they reach the 
collection region of the SRAM and thus can 
cause an upset. 


shown as 0.2 upsets. The 21 -Feb solar flare is 
clearly visible at day 52 in Figs. 7 to 9. 

1E+5 

| 1E+4 

= 1E+3 
w 

tn ie+2 

co 
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=> 1E+1 

T“ 
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n 1E +° 

IE-1 

45 65 85 105 125 145 

PRQ1 

ciW 4 ?i 9 xis TIME, t (days from Jan 1 , 1 994) 

Figure 7. Overview of the Zil upsets for S/C 
12 . Z3, Z5, Z4, XI. and Y1 SRAMs. 

The data shown in Fig. 8 was heavily censored 
for there were a great many system glitches of 
unknown origin on the ISA which caused massive 
SRAM upsets. These upsets are most likely due 
to noise on the power supplied to the RRELAX. 
Thus, a lesson learned is that the power supply 
must be stable especially during the 100 s 
stare period. 


SOLAR 

FLARE 



— o— Z2, 0.025 mm 
— Z3, 0.787 mm 
— a— Z 5, 1.649 mm 
— o — ZA, 2.311 mm 
-x- XI ,0.787 mm 
—f-Y 1, 0.787 mm 


■ f - A " t 



Figure 8. Overview of the Z^l upsets for ISA 
Z2, Z3. Z5. Z4, XI. and Y1 SRAMs. 


4. TEST RESULTS 

An overview of the S/C and ISA SRAM particle 
detector data is presented in Figs. 7 to 9 for 
Zsi and heavy ion particles. Zero upsets are 


As shown in the following figures, most upsets 
follow a pattern where the lightly shielded 
SRAMs have more upsets than the heavily 
shielded SRAMs. As seen in Figs. 7 to 9, this 
pattern is broken during quiet times when the 
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background radiation can upset the more heavily 
shielded SRAMs more often than the lightly 
shielded SRAMs. 



Figure 9. Overview of the heavy ion upsets for 
S/C 12 . 13 . 15 . Z4, XI. and Y1 SRAMs. 


The effect of shielding on the 21 -Feb solar 
flare is shown in Fig. 10. The shielding 
thickness is listed in the legend and the 
corresponding protons energies are listed in 
Table 1. The figure illustrates how the upsets 
decrease with shielding thickness. 



Figure 10. Effect of shielding on the upset 
response of the S/C Z2. Z3. Z5. and Z4 SRAMs 
during the 21 -Feb solar flare. 


The energy dependence of proton flux at the 
solar flare peak is shown in Fig. 11. This 
data was corrected for the reduction in the 


memory cell count during the stare cycle using 
Nq - -N t -ln(l - N/N t ) where N t - 18-4096 and N 

is the number of upsets shown in Fig. 10. The 
field-of-view was assumed to be re sr. As seen 
in Fig. 10 for 12 . N - 56,503 and N 0 - 107,202. 



Figure 11. Peak proton flux energy dependence 
for the 21 -Feb solar flare for S/C and ISA. 


As seen in Fig. 11. the flare was more intense 
for the S/C than for the ISA. In comparison 
with other worst case flares [7], this flare 
was several orders of magnitude lower in its 
intensity. 

The effect of orientation is shown in Figs. 12 



52 53 

ms TIME, t (days from Jan 1 , 1994) 


Figure 12. Effect of orientation on the upset 
response of the S/C Z3, XI. and Y1 SRAMs during 
the 21 -Feb solar flare. 
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Figure 13. Effect of orientation on the upset 
response of the ISA Z3, XI, and Y1 SRAMs during 
the 21 -Feb solar flare. 


and 13. The ordering of the data follows 
approximately the field-of-view of each SRAM as 
seen in Fig. 1. 

For the S/C, the field-of-view or solid angle 
for the Z-axis is approximately 2n; whereas, 
for the X- and Y- axes, it is closer to 7t. 
Thus the Z-axis should have the most upsets but 
the curves for the X- and Y- axes should be 
closer. For the ISA. the field-of-view is 
largest for the X-axis and smallest for the Y- 
axis. The ISA contains a large spherical tank 
the occludes the RRELAX field-of-view. The ISA 
results fit intuition better than the S/C 
results. 



Figure 14. Effect of thresholding on the upset 
response of the S/C Z2 SRAM during the 21 -Feb 
solar flare. 


The effect of thresholding on upsets is shown 
in Figs. 14 and 15. These figures indicate 
that the number of upsets decreases as the 
threshold increases. The upper curve is due to 
Z>1 particles and the lower curve is due to 
heavy ions. This shows how the SRAM is used to 
electronically discriminate between light and 
heavy particles. In Fig 14, there is a 
significant shift between the light and heavy 
particle responses. 



50 

OFF I 


51 52 53 

TIME, t (days from Jan 1, 1994) 


Figure 15. Effect of thresholding on the upset 
response of the ISA Z2 SRAM during the 21 -Feb 
solar flare. 

A comparison of the p-FET and SRAM response is 
shown in Fig. 16 during the passage of the S/C 
through the earth's radiation belts. First the 
SRAM is upset by protons in the proton belt and 
then the p-FET is dosed by electrons in the 
electron belt. The p-FET responds to total 
ionizing charge: whereas, the SRAM responds 
only to ion strikes and not electrons. 

The 21-Feb solar flare response of the S/C Z3 
SRAM was compared to the GOES -7 particle 
detector in Fig. 17. During the flare, the S/C 
was in moon orbit and the GOES -7 in 
geosynchronous orbit. This explains the one- 
hour shift in the responses. 
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Figure 16 Response of the S/C Z2 p-FET which 
detected the electron belt passage and the SRAM 
which detected the proton belt passage. 



Figure 17. Comparison of the GOES-7 particle 
spectrometer with the S/C Z3 SRAM detector. 


Two microflare series are shown in Fig. 18 for 
the 12 SRAM. This SRAM is very lightly 
shielded with a 25-/;m A1 shield which means 
that it responds on the average to 1.3 MeV 
protons. Note the data shown in the figure was 
not corrected for memory cell depletion during 
the stare cycle. 

The microflare at day 105 is apparently 
associated with the 21 -Feb solar flare. This 
feature was also seen in the SAMPEX's [8] MAST 
[9] data as seen in Fig 19. In this figure the 



Figure 18. Microflares appear in two cycles 
with a 26.5-day period. 
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Figure 19. A microflare appears at day 105 in 
this data from the SAMPEX MAST detector. This 
feature is displaced by two 26. 5 -day periods 
from the 21 -Feb (day 52) solar flare. 

Series #2 upsets are due to both protons and 
helium ions. Notice that there is no 
micro flare after the first 26. 5 -day period at 
day 78.5. The data shown i-n this figure was 
not corrected for detector dead time which will 
increase the flux near the solar flare peak. 

The Series #1 microflares, seen in Fig. 18. 
were also observed by the SAMPEX's LEICA 
instrument [10]. This series is composed of 
0.5 - 6.6 MeV/nuc helium atoms. 
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5. DISCUSSION 

The upsets, shown in Fig. 19. indicate that the 
SRAM has a dynamic range of five orders of 
magnitude relative to the 21-Feb solar flare. 
The response in passing through the earth's 
radiation belts, shown in Fig. 16. indicates 
that the SRAM does not respond to electrons but 
responds only to protons. 

As seen in Fig. 7. the SRAM has most upsets 
when it is virtually unshielded. The active 

area of the SRAM is 0.00172 cm 2 which is 
smaller than centimeter -size pin diode and 
surface-barrier detectors used in particle 
spectrometers for measuring space background 
radiation particles. Because the active area 
of the SRAM is small, its usefulness lies in 
analyzing higher density particles found in 
solar flares and the earth's radiation belts. 

The spectrometer aspects of the SRAM are 
determined by the shielding and V 0 . Shielding 

filters the external particle spectrum so that 
only those particles that fall into the energy 
detection window can upset the SRAM; see Fig. 

6. The offset voltage. V 0 . sets the upset 

threshold as depicted in Fig. 6 so the SRAM can 
discriminate between light and heavy particles. 
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ABSTRACT 

A custom designed SRAM was used to 
measure the thickness of integrated circuit 
over layers and the epi-layer thickness using 
alpha particles and a test SRAM. The over 
layer consists of oxide, nitride, metal and 
junction regions. 

INTRODUCTION 

It is more difficult to measure the thickness 
of integrated circuit layers than the width. 
Traditional methods such as ellipsometry 
require areas that are large compared to 
integrated circuit features, and cross 
sectioning techniques are time consuming. 
In this paper, alpha particles are used to 
measure the over layer and collection layer 
using a custom SRAM. 



Figure 1. STRV-RADMON chip 2.57 x 2.70 
mm 2 . 


This structure was designed as a radiation 
monitor, RADMON, and its chip-level layout 
is shown in Fig. 1. The device is currently 
flying on the Space Technology Research 
Vehicle (STRV-lb) in a Geosynchronous 
Transfer Orbit. This is a joint US-UK project 
designed to evaluate the effects of space 
radiation on advanced technology. 

The device is a custom designed 4-kbit 
SRAM with a traditional six-transistor cell 
[1]. A cross-section through the two 
inverters which form the memory cell is 
shown in Fig. 2. As 


voo 



Figure 2. Cross section of the RADMON 
memory cell showing a particle track 
creating hole-electron pairs beneath drain 
Dn2. 


indicated in the figure, particles create hole- 
electron pairs as they pass through the 
silicon and if enough charge is deposited in 
the collection layer, DX4, the cell will change 
state. In the experiments described here, the 
alpha energy is initially low enough so the 
alphas do not reach DX4. As the energy is 
increased the alphas enter DX4 and reach 
end-of-range before exiting DX4. At higher 
energies, the alphas pass through DX4 and 
reach end-of-range in the substrate. The 
amount of charge deposited in DX4 will vary 
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according to the Bragg curve. In general, 
the deposited charge will increase and then 
decrease as the alpha energy increases. 

The cell features an offset voltage, V p , 
connected to inverter #2. This voltage is 
used to increase the sensitivity of the SRAM 
to particle upsets by biasing the inverter 
relative to the cell's metastable point. When 
biased close to the metastable point, a 
particle that creates a small amount of 
charge in DX4 can upset the cell. 

In operation the SRAM is written with all 
ones. Then V Q is lowered and the SRAM is 
placed in the stare state for a period of time 
during which the particles flip the cells. 
Then, V c is raised and the number of zeroes 
read. 

As indicated in Fig. 2, the drain of the n-FET, 
Mn2, is connected to V p through p-FET Mp2. 
The drain, Dn2, of Mn2 is enlarged to 
increase its cross section to particles. 



Figure 3. Memory cell layout showing the 
enlarged drain, Dn2, in bold outline. 


This is shown in greater detail in the cell 
layout shown in Fig. 3 by the area outlined 
with a bold line. A close inspection of this 
region indicates that it can be subdivided 
into three areas each with a different over 
layer thickness. A cross section of the 
drain, Dn2, is shown in Fig. 4. There are 
three distinct areas: (a) no metal (MO), (b) 

metal-1 or metal-2, (M1&M2) or (c) metal-1 
and metal-2 (M1+M2) and these are listed in 
Table 1. The over layer consists of the 
insulator layers including the deposited 
oxides and nitrides, the metal layers, and the 
diffused junction. The thickness of the three 


regions and the epi-layer is the subject of 
this paper. 


NO 

METAL 
AREA, MO 


INSULATOR 


SINGLE 

METAL 

AREA, M1&M2 


METAL-1 


METAL-21 


DOUBLE 
METAL 
AREA, M1+M2 


METAL-21 


METAL- 1 


n-IMPLANT/ACTIVE/DlFPUSION 

p-EPITAXIAL REGION (COLLECTION LAYER) 


p-SUBSTRATE 

KAOMSm.MK 


Figure 4. Cross section through sensitive 
diode Dn2. 


Table 1. Over Layer Thickness 


REGION 

AREA 

EST 

ALPHA 

ALPHA 


pm 1 

DX3 

DX3 

DX4 



pm 

pm 

pm 

MO 

5.6 

4.0 

4.3 

3.3 

METAL-1 

7.5 

4.8 



METAL-2 

12.4 

5.4 



M1&M2 

19.9 

5.3 

5.2 

3.6 

M1+M2 

16.6 

6.3 

6.5 

3.6 

TOTAL 

42.1 






ALPHA ENERGY, E(MeV) 


Figure 5. Alpha particle range-energy curve 
for silicon and aluminum [2]. 


The range-energy dependence of alphas in 
silicon and aluminum is shown in Fig. 5. 
This graph indicates that there is little 
difference in the range-energy response for 
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alphas in silicon and aluminum. However, 
for other materials the range-energy relation 
must be computed as discussed below. 

The analysis for the over layer thickness 
used the energy-range curve for alphas in 
silicon. An algorithm was fitted to the data 
shown in Fig. 5. The algorithm is given in 
the Appendix. 

The CMOS over layers consist of a number 
of different materials and the range of alphas 
in these materials was converted to an 
equivalent range in silicon using the Bragg- 
Kleeman rule [3]: 

R1/R2 = (p2/p1)VA1/A2 (1) 

where p is the density and A is the atomic 
weight of material 1 and 2. In the analysis, 
material 2 is silicon and the non-silicon layer 
is material 1. The average thickness of each 
layer was obtained from the silicon broker, 
MOSIS. The estimated equivalent silicon 
thickness for each layer is listed in Table 1. 

The nomenclature used to describe the 
thickness and energy of the layers is shown 
in Table 2. Experiments were performed in a 
vacuum of less than 10‘ 4 Torr. Gold 
scattering foils were used to adjust the 
energy of the alpha particles. The energy of 
the alpha particles was determined using a 
pin diode detector and particle spectroscopy 
equipment. The energy of the alphas as they 
enter the silicon chip is given by E2. 


Table 2. Layer Notation 


Yfi. 

PARTICLE SOURCE 

EO 

AU* 

Y1. 

SCATTERING FOIL 

DX1 

El 

A 1 
YO. 

VACUUM 

DX2 

E2 


SRAM OVER LAYER 

DX3 

E3 

A 
Y A 

SRAM COLLECTION LAYER 

DX4 

E4 

AH 

SRAM SUBSTRATE 

DX5 


where N( = 4096 cells, ^ is the particle flux, 
A 0 is the memory cell cross section. In the 
limit where the number of upset cells is 
small compared with Nj, the cross section is: 

A c = R 0 /(<J>N t ) (3) 

where R 0 = dn/dt|{ q. For this SRAM, A 0 
is 42.1 pm z . The physical meaning of the 
cross section is the number of upset cells, 
dn, for a given 4>, Nfc and dt, the stare time. 
The stare time is the time the SRAM is 
exposed to the alpha particles. 

In these experiments, the alpha sources 
were Americium-241 and Gadolinium-148. 
The incident energy was adjusted using gold 
scattering foils with various thicknesses. 
The Am-241 flux was 145 alphas/cm 2 -sec with 
a stare time of 1000 sec. The Gd-148 flux 
was 87.8 alphas/cm 2 -sec with a stare time of 
1000 sec. 

The results seen in Fig. 6 indicate that the 
SRAM has a spontaneous (SPON) response 
when V Q is lowered below 1.143 V. The 
SPON curve has a normally distributed set of 
values due to the spread in cell parameters 
[1]. This curve is characterized by the mean 
value, V 0 su- When alpha particles strike the 
cell, the cells flip at higher V Q values which 
depend on the flux and energy of the 
particles. 
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Figure 6 Upset cross section of the SRAM 
as a function of V 0 and various alpha 
particles. 


EXPERIMENT 

The RADMON was fabricated in a 1.2-pm 
CMOS process and test results are shown in 
Fig. 4. These results were derived using the 
detector equation which describes the rate 
at which cells are upset: 


The Am-241 data is shown in greater detail in 
Fig. 7 and the Gd-148 data is shown in Fig. 8. 
This data shows that the cross section is 
significantly greater than A 0 = 42.1 pm 2 . 
This is due to the capture of charge that 
strikes active regions, such as M0 or M1&M2 
or M1+M2. 


dn/dt <t>-A 0 (Nt - n) (2) The results were analyzed at A 0 /2 and 

displayed as the square data shown in Fig. 9. 
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At A 0 /2, half the particle flux has sufficient 
energy to flip the memory cells due to 
energy dispersion of the alpha particles. 
The data is shown as the differential offset 
voltage given by the relation: DV 0 = V Q - 
V D sp- 

The data plotted in Figs. 7 and 8 shows four 
distinct regions if plotted as in Fig. 10. A 
three dimensional representation is shown in 
Fig. 11. Three distinct peaks are shown on 
the left-side of the graph and these have 
been correlated to the three areas of the 
memory cell. 



1.7 1.8 1.9 

OFFSET VOLTAGE Vo (V) 


Figure 7. SRAM cross-section versus offset 
voltage response to Am-241 alpha particles. 



Figure 8. SRAM cross-section versus offset 
voltage response to Gd-148 alpha particles. 


The analysis for the layer thickness is shown 
by the three solid lines given in Fig. 9. 
These curves were obtained assuming 
values for DX3 and DX4. These parameters 
were adjusted manually until a reasonable fit 
was obtained. The energy deposited in the 
over layer, DX3, is E3 - E4. The range-energy 


algorithms listed in the Appendix were used 
in these calculations. The layer thicknesses 
are listed in the inset of Fig. 9 and in Table 1. 
The thickness values listed in Table 1 are in 
good agreement with the thickness 
estimates derived from the manufacturer's 
target values. 

The solid curves shown in Fig. 9 have three 
regions. In the first region, particles have 
insufficient energy to reach DX4 and stop in 
the DX3. Particles with these low energies 
will not upset the cells and thus DV Q = 0. 
Particles with sufficient energy will stop in 
DX3 and deposit sufficient charge to upset 
the cells. This is seen as the rising portion 
of the curves in Fig. 5. Finally, particles with 
more energy pass through DX4 and stop in 
DX5. This is seen as the falling portion of 
the curves in Fig. 5. The point at which the 
alphas just reach DX5 is the point at which 
the slope of the rising portion of the curve 
begins to decrease near the top of the curve. 



1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 
cros4C02°kls INCIDENT ENERGY. E2 (MeV) 

Figure 9. Particle physics fitted to 
experimental data showing fitted values for 
the over layer thickness, DX3, and the 
collection depth, DX4. 


An additional layer, which was measured by 
this technique, is the epi-layer noted as DX4. 
Two values, 3.3 and 3.6 pm, are listed for 
DX4 in the inset to Fig. 9. 


DISCUSSION 

The use of a custom SRAM is shown to have 
the potential for measuring the over layer 
and epi-thickness of integrated circuit 
layers. This method is effective in 
measuring the thickness of both conducting 
and insulating layers. The method requires 
a knowledge of the alpha particle range- 
energy dependence and the conversion of 
the thickness of non-silicon layers to an 
equivalent silicon thickness. 

The memory cell used in this study has 
regions with different over layer thicknesses 
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and suffers from the collection of peripheral 
charge. A simpler structure should be 
designed with a single over layer thickness 
on each memory cell and the cells should be 
designed with a channel stop to prevent the 
collection of peripheral charge. The 
memory could be designed as a scribe-line 
monitor and used for end-of-line over layer 
and epi-thickness measurements. 



Figure. 10. SRAM cross-section versus 
alpha energy and V 0 values showing three 
major peaks and one minor peak. 



Figure 11. SRAM cross-section response 
due to three memory cell regions: no metal 
(MO), Metal-1 and metal-2 (M1&M2), and 
metal-1 plus metal-2 (M1+M2). 


APPENDIX 

The range-energy curves [2] for alpha 
particles in silicon shown in Fig. 4 was fitted 


with the following parabolic equation where 
E has units of MeV: 

Y = log 10 (E) = aX 2 + bX + c (A1) 

where X = log-jo(R) and R has units of pm. 
The coefficients of this equation are: 

a = -0.2867 
b = 1.4395 
c = -0.7151. 

The inverse equation was determined by 
solving the above quadratic equation. 

X = log 10 (R) = A(-1 + V[1 - B (C - Y)] (A2) 

where 

A = b/2a = -2.51046 
B = 4a/b 2 = -0.55343 
C = c. 

It is essential that the inverse equation be 
derived from the initial a, b, c fitting 
parameters using this procedure. 
Attempting to fit the inverse relation, that is 
the logR vs logE curve, leads to different 
parameters that ruin the analysis. 
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Clementine Dosimetry 
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and 
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Thus paper presents radiation dosimetry results from the radiation and reliability assurance experiments on the 
Clementine spacecraft and Interstate Adapter Satellite. The dosimetry instruments utilize Ion -dose- response />• 
channel field -effect -transistor and proton -sensitive static-RAM dosimeters. These dosimeters successfully demon- 
strated an order-of-mapnitude decrease in instrument weight and power from previous systems. The data confirm 
prelaunch predictions of the total-dose effects on charge-coupled-device dark current and on complementary- 
m eta b-oxide-sem icon duct or 1.2- and technology design and process parameters in the space environment. 

A solar proton event dunng the eart> phases of the mission allowed a comparison between the Clementine mea- 
sured proton fiuence and energy spectra and the Geostationary Operational Environmental Satellite 6 measured 
proton fiuence and energy spectra. These data confirmed the proton-spectrometer design and operation. 


Nomenclature 

a„ = second-order polynomial coefficients 

C y = upset capacitance 

D< - space-measured radiation dose 

D(T) = temperature -corrected space radiation dose 

E f = silicon bandgap energy 

F(E) = external proton frequently 

/ e = proton internal environment fraction 

= SRAM on-time fraction 
G = p-channel field-effect transistor (/>FET) 

temperature -independent circuit gain 
l D = complementary metal-oxide-semiconductor (CMOS) 

drain current 

l d ~ charge-coupled-device (CCD) dark -current damage 

factor, n A/cm 2 ■ krad 
lg = prelaunch CCD2 dark current 

K = charged-particle-produced ionization per unit energy 

loss in silicon, fC/MeV 
Kp = planetary magnetic index 

k - Boltzmann’s constant 

N = SRAM -measured number of counts per hour 
Nr = number of pixels per SRAM chip 
T s = space-measured temperature of z DUT board 
V f = CCD2 current -to- voltage conversion factor 
V DD - SRAM chip bias voltage 

V d = pFET temperature-independent threshold voltage 
damage factor. mV/krad 
V c - CMOS gate voltage 

V g = ground- measured voltage 

V p = SRAM adjustable offset voltage 

V, = space -measured voltage 

V T - CMOS threshold voltage 

V Td = CMOS threshold-voltage damage factor, mV/krad 

V Ti = CMOS space-measured threshold voltage 
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V Tl = CMOS threshold-voltage temperature factor. mV/C 
V,, = CCD2 extended pixel voltage mean 

A E = static-RAM (SRAM) proton energy window width 

AV r = SRAM offset voltage above metastable point 
a = SRAM pixel sensitive area or “cross section" 

£2 = Clementine SRAM field of view 


Introduction 

T HE mission goal of the Clementine Engineering Experiments 
Program was to evaluate the effects of the space environment 
on advanced Ballistic Missile Defense Organization (BMDO) tech- 
nologies being flight-tested on the main Clementine spacecraft and 
the Interstage Adapter Satellite (ISAS) The radiation and reliabil- 
ity assurance experiment (RRELAX) instrument contributed to the 
flight qualification of advanced BMDO technologies by perform- 
ing radiation dosimetry in the space environment. The RRELAX 
dosimetry' instruments utilized Jet Propulsion Laboratory (JPL) 
low-dose -response />-channel field effect transistor (pFET) and 
proton-sensitive static RAM {SRAM) integrated-circuit chip (1C- 
chip) dosimeters. These dosimeters, collectively called Radiation 
monitors (RADMONs). demonstrated an order-of-magnitude de- 
crease in instrument weight and power from previous systems. 

Dosimetry data for the Clementine spacecraft, on Feb. 21, 1994 
and April 27. 1995. and ISAS, throughout its mission, are presented 
in this paper. The ISAS pFET data confirmed prelaunch predic- 
tions of the total-dose effects on charged-coupled device (CCD) 
dark current and on Hewlett-Packard (HP) 1.2- and 0.8-^m tech- 
nology design and process parameters in the space environment. A 
solar proton event during the early phases of the mission allowed 
a comparison between the Clementine RRELAX measured proton 
fiuence and energy spectra and the Geostationary Operational En- 
vironmental Satellite 6 (GOES-6) measured proton fiuence and en- 
ergy spectra. These data confirmed the SRAM proton spectrometer 
design and operation. 

RRELAX Experiment Design 

Two identical RRELAX experiments flew on the Clementine 
missions — one on the Clementine spacecraft and one on the 
ISAS. The RRELAX experiments are autonomous microprocessor- 
controlled data handling and analog test systems. They were con- 
tained on three device under lest (DUT) boards. These boards, 
termed the x, y, and z DUT boards, are shown in Fig. 1 The z 
DUT board, which contains most of the RRELAX experiments, is 
mourned on the top of the RRELAX box normal to the z axis. 
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Fig. I RRELAX experiment layout in 4 x 4 x 1.5-in. box. 
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Fig. 2 Chip dosimeter 
(pFET on SRAM) and 
operating-mode drcuit 
diagram. 


The data presented in this paper were generated by the z DUT 
board RADMONs and CCD2. Each DUT board was covered with 
a 0.5-mil (A1 equivalent) aluminized kapton dust cover, and RAD- 
MON(Z2), which had no Kovar lid. was covered with an additional 
0.5-mil aluminized kapton dust cover. CCD2 had a 15-mil Kovar 
lid, RADMON(Z3) a 10-mil Kovar lid, RADMON(Z5) two 10- 
mil Kovar lids, and RADMON(Z4) three 10-mil Kovar lids. The 
ISAS RRELAX was also covered by a 6-mil (A! equivalent) ther- 
mal blanket. [0.5, 6, 10 mil = 13, 15, 25 pm] 

The z-DUT-board RADMONs were fabricated through MOSIS 
at HP Corvallis in 1 .2-pm technology, and the JPL CCD2s at Loral. 
A pFET dosimeter and SRAM proton counter was contained on 
each 1.2-Aim RADMON IC chip (Z2, Z3, Z5, Z4) on the : DUT 
board, shown in Ftg. 1. 

Design of the pFET Dosimeter System 

The pFET dosimeter system was designed to measure the total 
dose from all radiation types. The sensors were placed behind four 
different shielding thicknesses to produce dose-depth curves for the 
Clementine and ISAS experiments 
The dosimeter consisted of a pFET that was operated in a 
constant-current mode (Fig. 2). This allowed the temperature effects 
on channel mobility and threshold voltage to cancel each other. The 
objective of this design was to minimize temperature effects and em- 
phasize dose dependence. The pFETs have a closed geometry gate, 
eliminating the source -to-drain bird’s-beak leakage path. Dosime- 
try measured via the pFET threshold- voltage shifts was influenced 
by two second-order effects: 1 ) the temperature sensitivity of the 
transconductance factor and the threshold voltage V T and 2) the 
source-to-drain bird's beak leakage current. The goal in developing 
the advanced pFET dosimeter was to minimize or eliminate these 
effects. 1 

The pFET dosimeters were only powered during readout and 
accumulated dose while in an unpowered stale. The four pFET 
dosimeters were read once per hour, and their data from the previous 
hour overwritten. Because pFETs are integrating dosimeters and 
have a larger radiation damage factor when unpowered, this mode 
of operation is satisfactory. 


Design of the SRAM Proton Spectrometer System 

The SRAM proton spectrometer system was designed to measure 
the proton energy spectrum and fluence for use in computing proton 
fluence and energy at the Clementine experiments. Four shielding 
thicknesses were used to pick ofF the external proton environment 
at four different energies, allowing the energy spectrum to be mea- 
sured. T\vo energy bins were used to allow heavy-ion- and proton- 
induced nuclear reactions to be subtracted from the proton data. The 
threshold of the proton-sensitive bin was set just below the proton 
injection peak to optimize the energy measurement. 

The SRAM detector was fabricated in 1.2-pm, /i-well, double- 
metal CMOS. A schematic diagram of the SRAM cell, or pixel, 
is shown in Fig. 3. This cell differs from that of a standard six- 
transistor SRAM cell in three ways: 1) the source of the pFET, 
Mp2, is connected to an adjustable offset voltage V 0 instead of V DD , 
to provi de co ntrol of the cell's upset capacitance; 2) the drain area 
Dn2 of/iFET Mn2 has been enlarged by a factor of 4 over minimum 
to enhance upset rates, thus reducing the measurement time; and 3) 
the cell is unbalanced by widening Mn2 over minimum to enhance 
its SEU sensitivity vs V a . 

In operation, all the memory cells were written into a sensitive 
state, where Mn2 was turned off and Mp2 was turned on, connecting 
V 0 to the bloated drain, Dn2. V DD was then lowered to 3 V and V 0 
was lowered below V DD = 3 V, allowing the SRAM to accumulate 
upsets at a given V 0 value. Thereafter V Q and V DD were returned to 
5 V and the cells read to determine the number of upsets. 2 

Four SRAM chips were mounted to face along the z axis of the 
experiment, two facing the x axis and two facing the y axis, for a 
total of eight chips, as shown in Fig. 1 . The SRAM portion of each 
chip was controlled using a single V DD and V 0 for all eight chips. All 
z-axis SRAMs, and one SRAM on the x axis and one on the y axis, 
were fabricated in 1 .2-pm technology. One SRAM on the x axis 
and one SRAM on the y axis were fabricated in 0.8-pm technology. 

The SRAMs were made sensitive to upset by applying 3-V V DD 
and V a voltages that have the following A V p values above the spon- 
taneous flip voltage: 


AV p = 0.150 V 
AV P = 1.0 V 

The spontaneous flip voltage is the value of V n where the SRAM 
cell becomes metastable and spontaneously flips. The offset voltage 
was set above the spontaneous flip voltage by an amount AV, to 
sensitize the SRAM cell to protons, heavy ions, and nuclear reactions 
(AV p = 0.150 V) or to heavy ions and nuclear reactions only 
{AV p a 1.0 V). The integration period for these sensitive states 
was a constant of 100’s. At the end of each integration period , 
each of the eight SRAMS was evaluated for bit changes, or upsets, 
and the total number for each chip stored in the data structure. The data 
structure and software flow are shown in Fig. 4. 

When the evaluation of SRAM response was complete, software 
accumulated the upset totals. Each chip had its own accumulation, 
or running total, for each value of A V p .The SRAM contents were 
then restored to their starting state and the second value of A V, set 
for the next integration time. The experiment software alternated 
between the two A V p values on a 1:1 ratio. This was done by the 
SWITCH VOLTAGE box in Fig 4 
The 100-s integration records together with their time stamp were 
called the high-resolution SRAM data record. Software maintained 
each record for one hour before disposing of it. This generated 18 
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Fig. 3 Schematic diagram af the 
SRAM ceil showing the placement 
of V» and the Moated * -drain, 
DnZ. 
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Fig. 4 Software flow diagram showing SRAM data -band ling strategy. 


consecutive sets (18 periods per hour) of 16- word records (one for 
each SRAM and AV'^ value), which were stored in memory These 
records formed the highest- resolution data set per SRAM and per 
AV' r value 

A high-resolution SRAM record that was about to be dis- 
carded was accumulated into a medium-resolution record The 
medium- resolution record accumulated, or summed, the 18 high- 
resolution records for each SRAM and A V p value Each completed 
medium -resolution record holds integration results for a one-hour 
period for each SRAM and AV, value. The medium resolution 
records time entry represented the lime when the record was formed 
The software retained medium resolution records for 36 h for a total 
of 576 medium-resolution records This group of records formed 
one medium-resolution data set per SRAM and per & V p value 

Each SRAM, for each A V p value, measured the total fluence over 
the whole mission. This process used memory locations that were 
only added to and never overwritten The data set was called the 
low-resolution data set. 

The SRAMs were calibrated during each medium resolution data 
set by measuring the numbers of stuck zeros and ones. The temper- 
atures of the analog board and the jc, y, and z DUT boards and a 
temperature circuit calibration measurement were then stored, over- 
writing the previous calibration data record. 

Ground Test Data 

Selected pFET dosimeters from the flight fabrication run were 
calibrated with the JPL ^'Co source The W) Co irradiation was done 
at room temperature to a total dose of 1 00 krad(Si ). Tbe pFETs were 
tested at dose levels of 20, 40, 60, 80. and 100 krad and at -30. 20, 
and 50 : C at each dose level The pFETs were measured using an 
HP 4062 parametric test system in an oven that was also cooled 
with liquid nitrogen so that total-dose effects could be measured 
as a function of temperature. A temperature- independent thresh- 
old voltage damage factor V d of 3.68 mV/krad was measured The 
pFET was operated in the saturation region, which was ensured 
by connecting the gate to the drain as shown in Fig 2, and at the 
temperature-independent point. 1 

A JPL CCD from the flight fabrication run was ground-tested 
to 4.8 krad(Si) in the JPL w> Co source at room temperature and 
annealed at 25 C for 72 h. The CCD was characterized in the 
JPL CCD laboratory, and a dark -current damage factor I d of 7.38 
nA/cm 2 krad at 25 C was measured / 

An SRAM from the flight fabrication run w'as calibrated using 
the Caltech Tandem Van de Graaff proton accelerator with pro- 
tons at 0.75. 1.0. and 2.0 MeV. The SRAM ground-test response 
was analyzed using the TRIM (transport reactions in matter) com- 
puter code. 4 Calibration data are shown in Fig 5. The model used 
in the calibration assumed a 5-pm overlayer (dead layer) and a 
7- Jim charge collection depth (depletion depth). The proton cal- 
ibration data points are the mean values of the offset voltage, 
AV,, where one-half the protons hitting a cell sensitive volume, 



PROTON ENERGY. E 2 (MeV 

Fig. 5 Clementine 1 2-^iro- technology proton calibration curv es show- 
ing the flight energy threshold and the 1-MeV-wide energy window. 


Dn2. cause the cell to flip. The proton data calibrated the 1 2-pm- 
technology SRAM at an upset capacitance, C„ y K — 2 667 MeV/V. 
where K =44.2 fC/MeV for charged-panicle-produced ionization 
per unit energy loss in silicon. The measured energy A E 4 . shown in 
Fig. 5. is given by A E 4 = ( C u /K ) A V p . The 1 -MeV energy window, 
shown in Fig 5, is the measure of the SRAM response to protons 
in the space environment 2 

Space Data 

The ground-test-measured M, Co damage factor. V d = 3 .68 
mV/krad. was used to compute the space- measured pFET total dose 
The downlinked numbers were bits, with one least-significant bit 
equal to 1 .22 mV, from a 1 2-bit ADC measuring between 0 and 5 V: 

V, - V f 

- — krad(Si) (1) 

where V , (mV) is the space- measured pFET voltage. (mV) is the 
ground- measured voltage, and G = 32.7 is the circuit gain The zero- 
dose voltages V g (mV) and aluminium-equivalent shield thickness d 
(mils) for each flight RADMON on the ISAS RRELAX at 7 = 25 C 
and on the Clementine spacecraft RRELAX at T - -2 4! C are 
listed in Table 1. 

The p-FET total-dose flight results for the Clementine spacecraft 
on April 27. 1995, and the ISAS on April 20, 1994 (ISAS 1994 day 
1 10), are shown in Fig. 6. The Clementine spacecraft total dose after 
457 days in space is less than the ISAS over its mission life. On the 
Clementine spacecraft, the RADMON dosimeter-circuit dynamic 
range was exceeded for pFET 72 because G = 32.7 for that circuit 
and a total dose of 142.95 krad (Si) was extracted with the CMOS 
experiment circuit for that pFET. The ISAS pFET-dosi meter mis- 
sion dose profile is shown in Fig. 7. The ISAS is in a 2. 13-day lunar 
transfer orbit that penetrates the Earth’s radiation bells, causing the 
rapid increase in total dose. 

The ISAS RRELAX CCD experiment was shielded by 52 mils (1.3 
mm) (A 1 equivalent). Power-law dose-depth curves woe fit to the 
pFET data as a function of shield thickness. The pFET-measured dose 
at CCD2 was then extrapolated from these curves as shown in Fig. 8. 

The CCD2 sensor chips were 512-pixel line arrays manufactured 
at Loral in /i-channel metal-oxide system (nMOS) buried-channei 
technology. The RRELAX CCD experiment measured the dark- 
current -generated voltage in the last 8 pixels, pixel 505 through 51 2. 
and the extended pixels with no photosensors. 5 1 3 through 5 1 9 The 
data presented here had a photosensor dark-current integration time 
of 1 s. ISAS CCD2 data showing the prelaunch pixel voltage mean 
V g and the extended pixel voltage mean (for ground and space data) 
Vo are presented in Fig. 9. The ISAS CCD2 had a prelaunch dark 
current measured on the flight CCD 3 of l g = 1 4 nA/cnrr at 25°C 
and a voltage conversion of V ( = I f /(V g - V () ) =11.7 nA/cm 2 V. 

The 25~C dark -current damage factor l d = 7.38 nA/cm 2 krad 
was used to compute the space-measured total dose: 

( 2 ) 

Id 


73 













Table i Zero-dose pFET grmmd vefcagm V g md afczmioium 
equivalent shield thickness 4 far tnck flight RADMON oa the 
ISAS RRELAX at r = 25°C Md m the CkMDdBe 
spacecraft RRELAX at T = -2.41°C 


RADMON 

pFET 

Spacecraft 

ISAS 

> 

E 

4, mils 

V v mV 

d , mils 

22 

129.9 

1 

621.8 

7 

Z3 

74.8 

31 

3354 

37 

Z5 

274.8 

61 

592.2 

67 

Z4 

381.0 

91 

494.0 

97 



SHIELD THICKNESS, d (mils(AI)) 

Fig. 6 Total radiation dose on RADMON pFET 22 as a function of 
ahuninhun equivalent shield depth for the Gementine spacecraft on 
April 27, 1995, and the ISAS on April 20, 1994. 



Fig. 8 ISAS dose-depth curves showing pFET dose at CCD2. 



CCD PIXEL NUMBER 



ISAS TIME, GMT (days 1994) 

Fig. 7 Mission dose profile for the ISAS RADMON pFETs in a 2.13- 
day lunar transfer orbit and the associated temperature variation. 


Fig. 9 A sample of the RRELAX ISAS CCD2 data showing the 
prelaunch, zero-dose pixel voltage mean V g and the extended pixel volt- 
age mean V# for ground and space data. 



ISAS TIME, GMT (days 1994) 


where V t (V) is the space- measured CCD voltage mean for pixels 
with photosensors, pixels 505 through 512. The tool dose at 25 : C. 
0(7*). in kilorads (Si), is given by 


DiT) = 




( 3 ) 


where Eg= 1 .21 eV, T s is the z-DUT-board temperature, D s is the 
space -measured dose from Eq. (2). and T = 25 C. lemperacurc- 
corrected CCD dose data arc shown in Fig. 10. The temperature 
correction given by Eq. (3) is due to temperature changes in 5 the 
intrinsic carrier density associated with bulk current. 

The CCD2 dark-current measured dose is compared with the 
pFET measured dose at the CCD in Fig. 1 1 . Two prominent geo- 
physical events were observed during the mission. The first, a solar 
proton event on day 52 of 1994 (Feb. 21,1 994). had a planetary mag- 
netic index Kp of 7+ and was observed at a number of spacecraft 
(the proton observations will be discussed shortly). The increased 
dose rate after the second major event (the sudden commencement 
of a storm followed by a brief but intense geomagnetic storm on day 
107) is due to trapped -electron- belt heating on days 106 and 107 of 
1994 (April 16-17. 1994). The increase implies more and higher- 
energy trapped electrons along the ISAS orbit Consistent with this 
interpretation. Kp = 8+ for the event, indicating a large increase in 
the trapped ring current and subsequent electron-belt heating. There 
was. however, no large solar proton event associated with the event 74 


Fig. 10 CCD temperature corrected [7X71] and uocorrected [D t \ dose 
data as a function of temperature. 

The ISAS pFET dosimeter, operated at its temperature- 
independent drain current, measured a total dose of 1 1.83 (krad) 
on day 109.99. The ISAS pFET complementary metal-oxide- 
senv conductor (CMOS) experiment, using the same pFET oper- 
ated at the five different drain currents shown in Fig. 12. measured 
a total dose, Z>(7*) given by Eq. (6), of 1 1.47 krad on day 109.99. 
The space- measured threshold voltage V Ts is computed with Eq. (5) 
from second-order polynomial coefficients generated by fitting flight 
data 1 on the drain current //> vs the gate voltage, V c with Eq. (4) as 
shown in Fig. 12: 


y/To — do + d] Vc — 02 

(4) 

+FZ:) 

(5) 

v n + v Tl(Tl - T) 
V T , 

(6) 


where V-n = 1 .24 mV/ : C, T = 25 r C, and V r = 883.7 mV measured 
on the ISAS pFET prior to launch, and V Td — 1 .674 mV/krad mea- 
sured on WI2P4C05 for the HP 1.2-pm technology. 1 The CMOS 
experiment pFET dose verifies the dosimeter design and ground 
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Fig. 11 CCD dose D\T \ and pFET dose at the CCD, showing in- 
creased dose rate after day 106. 
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Fig. 1 2 ISAS CMOS HP 1 .2- pm p FFT gate voltage shift on day 109.99, 
April 20, 1994. 
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Fig. 13 Clementine CMOS HP 12-^m pFET Z2 gate voltage shift on 
April 27, 1995. 


Table 2 Externa) environment energy windows and internal 
environment fractions f t for Clementine spacecraft, 
as a function of Kovar shielding 


Kovar shields 


RADMON 

pFET 

Thickness, 
mils mm 

f min- 
MeV 

£ mu • 

MeV 

Mean f t 
<0 7-17 MeVi 

Z2 

0 0 

2 16 

3 16 

0 261 ± 0 031 

Z3 

10 0.25 

1 I 81 

12.81 

0 072 s 0 009 

Z5 

20 0.5 

16 % 

17 96 

0,04" ± 0 00*" 

Z4 

30 0.75 

2! 04 

22. 04 

0 037 T 0 006 
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Fig. 14 Comparison of Clementine data with GOES-6 data during 
Feb. 21, 1994, solar proton event. The proton spectrometer sensitivity 
of one count per hour is shown. 



PROTON ENERGY, E (MeV) 

Fig. 15 Clementine spacecraft and GOES-6 external environment 
proton energy spectra on Feb, 21, 1994. 

pixel sensitive area has an as-drawn cross section a of 42 12 /<nr. 
The spectrometer hourly frequently £(£ ). in p/errr MeV h is 
given by 


calibration methodology The same equations are used to extract 
the pFET Z2 total dose on the Clementine spacecraft on April 27, 
1995, shown in Fig 13. where T s = 3.32 C and T = — 2.41 C. 

Vy = 863.9 mV measured on the Clementine pFET poor to launch, 
and the other variables are the same CMOS radiation damage fac- 
tors for n- and p-channel transistors in HP 1.2- and 0.8-pm tech- 
nologies are being utilized in radiation test- structure designs for the 
active-pixel image sensors (APS) 

SRAM detector proton data from the Feb. 21. 1994, solar pro- 
ton event have been compared with GOES-6 proton data. First, the 
Novice code was used to compute the proton environment as a func- 
tion of energy inside the SRAM shields Table 2 lists the external 
environment energy windows. £ mir , to E^, in the 1 -MeV- wide 
energy window measured in Fig 5 The mean values f c of the inter- 
nal environment fraction inside the shields are also listed in Table 2 
Next the environment fractions were computed with the Novice code 
from a 27r-sr omnidirectional fluence of 1.96 x 10 y p/cm* MeV at 
all energies outside the shields The proton fluence was then reduced 
by the factor f c inside the shields 
The SRAM spectrometers are sensitized to protons for 100 s, 
every other 100-s period, for one hour, giving an on-time fraction 
f m of 0.5. The energy window width A E (Fig. 5) is 1 MeV The 
instrument was designed with a 2rr-sr field of view Q. The SRAM ^ 


F{E)= ll— J »L —) | , 7 ) 

cQAEf t f m \N t - N/18 ) 

Each pixel can only count one proton in each 100-s proton-sensitive 
period, and there are 18 sensitive periods each hour There are N T = 
4096 pixels on each SRAM chip, and N is the measured number 
of counts per hour in each chip The Clementine spectrometer data. 
F(E) from Eq. (7). for the Clementine SRAM Z2 and Z3 proton 
energy windows and the GOES-6 hourly fluencies for the proton 
energy windows PI and P3 are plotted in Fig 14 for comparison 
The GOES-6 external environment proton energy windows. £ mir 
and £ m „, are PI : 0.6-4. 2 MeV; P2 4 2-8.7 MeV. P3: 8 7-14 MeV; 
and P4: 15-44 MeV The Clementine and GOES-6 energy spec- 
tra for the total measured fluence on Feb 21. 1994. are shown in 
Fig. 15 The data-point energies are taken at the center of the en- 
ergy windows. ( E mir 4- £ma* >/ 2. listed in Table 2 for the Clementine 
instrument and at the beginning of this paragraph for the GOES-6 
instrument 

Dosimetry Conclusions 

The primary objective of the RRELAX experiment was to monitor 
the radiation environment for the Clementine mission Two conclu- 
sions have been presented about this objective First, the ISAS CCD 



and CMOS experiment measurements were shown independently 
io agree well with the pFET dosimeter measured total dose This 
demonstrates that the results of the ground radiation tests can be 
used to accurately predict spacecraft sensor and electronic compo- 
nent radiation degradation. The test- structure approach employed 
also allowed direct measurement of radiation effects in space and 
verified our predictive understanding of the space environment’s 
radiation effects on spacecraft sensors and electronic components. 
Secondly, the SRAM proton spectrometer successfully measured 
'he proton fluence and energy spectrum. The data agreed quite well 
with GOES-6 data for the same event We therefore feel confident in 
stating that the dosimetry system, RAD MON. has accurately tracked 
the proton fluence and energy and total dose for the Clementine En- 
gineering experiments. 
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Abstract 

This paper describes an integrated space environmental moni- 
toring system for geosynchronous satellites. The system provides 
measurements of surface charging, internal charging, and total dose 
radiation from the space environment. It is installed on seven of 
Martin Marietta’s commercial communications spacecraft, which 
are scheduled for launch in the near future. There are two compo- 
nents: an internal charge monitor, and a surface charge detector 
Total weight is less than one pound; power consumption is about 

0.5 W. Charging flux information from these detectors will 
permit evaluation of the effects of internal charging on electronic 
PC boards and will lead to the development of reliable, low- 
weight methods to detect on-orbit spacecraft charging events and 
to control charging and discharging. 

1. INTRODUCTION 

The uncertainties in predicting the effects of energetic par- 
ticles on satellite systems remain a concern to all spacecraft de- 
signers and users. To assist the development of effective mitiga- 
tion techniques, a lightweight, low-power monitoring system has 
been designed and built by Martin Marietta Astro Space (now a 
division of Lockheed Martin) and Jet Propulsion Laboratory. The 
system consists of two types of instruments — (1) a radiation 
dosimeter calibrated to measure internal charging and (2) a sur- 
face charging detector. Both instruments will produce continu- 
ous telemetry data throughout the mission. 

The system will provide a means of monitoring external and 
internal charging events aboard satellites. It will be very useful 
in determining the effectiveness of surface ESD control mea- 
sures, such as conductive coatings, and of protective measures, 
such as grounded metal shielding, developed for mitigating in- 
ternal charging effects on PC boards and cables. These evalua- 
tions will point the way to future design solutions and to the 


development of realistic design requirements for control of sur- 
face charging and internal charging. 

The surface charging monitor was designed to measure charg- 
ing by magnetospheric plasma electrons with energies from 5 
keV to 20 keV and with worst-case current density from 0. 1 to 1 
nA/cm 1 2 . This range was chosen because it has a probability of 
occurrence of 99% [1]. The charge monitor will indicate when 
such an environment is present at the spacecraft’s orbital loca- 
tion. Two of these units are being mounted in different places on 
seven of Martin Marietta’s commercial communications satel- 
lites scheduled for launch in the near future. 

The internal charging detector was designed to monitor flux 
of, and total dose imparted by, electrons with energies from 200 
ke V to 6 MeV. This detector will indicate when such an environ- 
ment is present at the spacecraft’s orbital location and will mea- 
sure the flux of these electrons at various locations inside the 
spacecraft structure. The key component is a dosimeter designed, 
fabricated, and tested by the Center for Space Microelectronics 
Technology, Jet Propulsion Laboratory of the California Insti- 
tute of Technology [2, 3]. Two of these units are being mounted 
on seven Martin Marietta satellites. 

2. DESIGN DESCRIPTION 
2.7 Surface Charging Detector 

The surface charge detector has two components - the charg- 
ing plate and the associated control electronics. The charging 
plate is mounted to the spacecraft structure outside the space- 
craft and is exposed directly to the space environment. The con- 
trol electronics unit operates from the +/- 15-volt power lines. It is 
located inside the spacecraft with the charging plate connected to 
the input of the instrumentation operational amplifier through a 
high-input-impedance buffering RC circuit as shown in Fig. 1. The 
operational amplifier was configured to provide a 0- to 5-volt out- 
put compatible with the standard telemetry format. 
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Three critical choices during the design of the detector were: 
charging plate surface material, charging plate substrate, and 
mounting adhesive. The selection parameters for the surface 
material and mounting adhesive included electrical conductiv- 
ity, high voltage leakage current, thermal expansion coefficient, 
thermo-optical properties, and space environmental stability in 
the operational temperature range. Surface materials that were 
considered for the charging plate included S 1 3GLO white paint, 
Z306 Chemglaze black paint, black anodize, and alumina. Ad- 
hesives included RTV 566, CV2946, CV2943, and alumina- 
loaded Solithane 1 13. All candidate surface materials and ad- 
hesives were tested in a monoenergetic electron beam for their 
charging potentials, charging and discharging time constants, 
and compatibility with the readout circuits. A detailed descrip- 
tion of the test procedure and test facilities is given in the refer- 
enced documents [4-6]. 

The result of the selection process is a charge collector in the 
form of a 2” by 2” (5x5 cm) aluminum plate painted with 
Chemglaze Z306 dielectric paint bonded to the spacecraft with 
nonconductive CV2946 adhesive. The plate is charged by 
impinging electrons. The charging potential on the painted 
surface can vary from slightly positive to several hundred volts 
negative, with a corresponding potential on the underlying 
aluminum plate-that is, at the input to the buffer assembly in 
Fig. 1— varying from +5 to -25 volts. 

The sensor was calibrated with 5-keV, 10-keV, and 20-keV 
electron beams (beam current densities from 0.25 nA/cm 2 to 
1 nA/cm 2 ) over a -40°C to 60°C temperature range. Typical cali- 
bration curves are shown in Fig. 2. Because in space spacecraft 
ground floats with respect to ambient, these calibration curves 
are only directly applicable to the case of shadowed detectors 
while the rest of the spacecraft is in the sun. For this condition 
spacecraft ground floats only about a few volts positive. 


For the conditon of the spacecraft in eclipse or for the case 
of the sunlit sensor, NASCAP analysis is required to develop 
predictions for sensor response. Several qualification tests were 
performed on all sensors prior to their integration into the space- 
craft. These tests included radiation testing (400 megarads total 
dose) and thermal cycling (5000 cycles from +60°C to -60°C) 
on the charge plate, plus radiation testing to a total dose of 30 
kilorads on the readout circuit. ESD and DC overvoltage tests 
were also conducted on the readout circuit. 



Fig. 2 Surface charging sensor calibration (T=25°C) 


2.2 Dosimeterllnterrtal Charge Monitor 

The dosimeter within the charge monitor comprises four iden- 
tical radiation-sensitive p-FETs, each shielded with a different- 
thickness Kovar lid. The thinly shielded transistors correspond 
to components near the surface of a spacecraft; the thickly 
shielded transistors correspond to components deep within tfte 
spacecraft. Shield thicknesses were selected by trading FET 
lifetime vs. sensitivity to dose. The directly sensed quantity is 
ON-channel voltage drop at the temperature-independent 
current. This is related by calibration to total dose. Total dose 
can then be related by calibration to incident electron fluence. 
Differentiating the fluence with respect to time gives charging 
flux. The outermost transistors have a higher charge rate, which 
leads to better resolution, but, since the ON-channel voltage 
drop saturates at a certain total dose, this also leads to a shorter 
lifetime. That is, the outermost transistors will probably 
saturate before the end of the satellite's mission. 

The microchip layout developed on this effort is shown in 
Fig. 3. It consists of seven components listed in Table 1 includ- 
ing p-FETs P0, PI, P2, andP7, each with different shields (Note: 
1 mil = 0.00254 cm). These shields were fabricated in the Kovar 
lid by chemical etching. The microchip, shown in Fig. 3, was 
fabricated in a 1.2-mm CMOS rad-soft process. The chip was 
packaged in a 16-pin flat pack. The equivalent circuit for the 
chip, shown in Fig. 4, has a multiplexer and decoder used to 
place each component selectively in the feedback loop of the 
operational amplifier circuit shown in Fig. 5. The circuitry was 
designed to force a current at the temperature independent point 
through the p-FET during measurement. When not operating, 
the p-FETs are unbiased. 
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ICM547Q1.PCX VDD OPl GND2 


Fig. 3 Microchip with four p-FETs, PO, PI, P2, and P7. 
Each P-FET is located behind a different shield whose 
thickness is listed in Table 1. 


Table 1 

Microchip Devices 


Device 

Type 

Shield 
mil (^tm) 

P0 

Gate Oxide p-FET 

6(150) 

PI 

Gate Oxide p-FET 

0 

P2 

Gate Oxide p-FET 

2(50) 

P7 

Gate Oxide p-FET 

10(250) 

N3 

Field Oxide n-FET 

not used 

N5 

Gate Oxide n-FET 

not used 

R4 

Metal- 1 Resistor 

10(250) 

E6 

External Device 

NA 



10W4713PLT 

Fig. 4 Microchip schematic 


Dose levels and shielding for each of the p-FETs are shown in 
Table 2. In addition to the Kovar lid, the outside dosimeter is 
shielded by the thermal blanket and the inside dosimeter by a 
spacecraft panel. The equivalent shield thicknesses are shown 
in Table 2. 


Table 2 

The 250 krad time for each p-FET 
in the two dosimeters 


Device 

PI 

P2 

P0 

P7 

Outside Dosimeter 





Lid Thickness (in mil KOVAR) 

0 

2 

6 

10 

Shield Thickness (in mil KOVAR) 

2 

2 

2 

2 

Total Thickness (in mil KOVAR) 

2 

4 

8 

12 

250 Krad Time (months) 

0.7 

1.6 

5.5 

13.3 

Inside Dosimeter 





Lid Thickness (in mil KOVAR) 

0 

2 

6 

10 

Shield Thickness (in mil KOVAR) 

10 

10 

10 

10 

Total Thickness (in mil KOVAR) 

10 

12 

16 

20 

250 Krad Time (months) 

8.9 

13.3 

26.4 

45.9 


[SI units on p. 83] opi ext cpo vo 



The calibration of the p-FET must account for both dose and 
temperature effects. The analog voltage, VA, shown in Fig. 5 is: 

VA (T, D) = VAoo + VAp (T - To) + VO D G D (I) 

where VAoo is the initial output voltage determined from ground 
tests, VA t is the temperature coefficient determined from ground 
test data, VO D is the dose coefficient determined from Co-60 
ground tests, G = Gain = R13/R7, and To = 300K. D is the dose 
in kilorads; T is the temperature in K. The dose is calculated 
from: 

D = [VA - VAoo - VAp (T - To)]/(VO D G) (2) 

The p-FET output voltage, VO (Figs 4 & 5), follows from 
the square law behavior of the FET operating in saturation [2]: 

VO = VT-(2-ID/p) 1/2 (3) 

where VT is the threshold voltage, p is the transconductance, 
and ID is the drain current chosen at the temperature indepen- 
dent point. This expression is used to determine the VO-dose 
calibration curve given VT, p, and ID as discussed below. 

The temperature-dose expression for the VT is: 

VT = VToo + VT T (T- To) + AVT[1 - exp(-D/Do)] (4) 
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where VToo is the threshold voltage at To and D = 0, VT t is the 
threshold voltage temperature coefficient, Do is the VT dose 
coefficient, and AVT is the maximum change in VT when the 
radiation reaches infinity. When D= Do and at constant T, the 
VT increases by AVT[l-l/e]. The expression for VT indicates 
that VT is rate limited, which depends on the filling of a finite 
supply of gate-oxide hole traps [7]. The VT equation was solved 
using the least-squares method to determine VToo, VT t , and 
AVT. The parameter Do was found using an optimization tech- 
nique that maximizes the least-squares correlation coefficient. 

The dose-temperature expression for p is : 

P = Poo(T/To)" n /[l+ (D/Dm)(T/To)' n ] (5) 

where Poo is P evaluated at To and D=0, n is the P temperature 
coefficient and Dm is the P or mobility dose coefficient. For- 
mula (5) is equivalent to 2.9.1 in [8] with an added term ex- 
pressing radiation dependence. When D = Dm at T=To, then p 
is reduced by 50%. This equation was formulated by combin- 
ing the P temperature and dose dependence as the sum of recip- 
rocals. 

The p-FET drain current at the temperature independent op- 
erating point, IDm, is found by differentiating the VO expres- 
sion with respect to temperature. Setting the resulting expres- 
sion equal to zero at the measurement temperature, Tm (which 
is the expected on-orbit operating temperature), and D=0 leads 
to (6): 

IDm = 2(prno) 3 (-VT T /p Tmo ) 2 (6) 

where p Tm0 = (-n/Tm)pmo and pmo = (Tm/To)" n . 

The VT, p, and VO results are presented graphically in Figs. 
6 through 8. The curves were obtained by fitting the data shown 
in Fig. 6 to Eq. 4 and the data shown in Fig. 7 to Eq. 5. Thus, the 
p-FET results can be characterized by seven parameters, namely 
VToo, VT t , AVT, Do, Poo, n and Dm. These parameters are 
listed in Table 3 for transistor P0. The device-to-device unifor- 
mity was excellent. 

VO values were obtained from Eq. 3 using Eq. 4 for VT, Eq. 
5 for p, and Eq. 6 for IDm. The VOmo values for each of the 
four p-FETs, are within 0.1% of each other. The VO values are 
plotted in Fig. 8 for two different values of IDmo. The curves 
for the P0 design value of IDmo = 88 pA show that VO is fully 
temperature compensated at D=0. The curves for the operating 
value IDmo=100 jiA show that temperature compensation 
occurs near 100 krads and that VOmo=l .507 V. These parameters 
were used as the target design parameters, thus the 100-pA curve 
is designated as the calibration curve. The samples from the 
flight lot of dosimeters were radiation tested and calibrated using 
a Co-60 gamma source at -25C, +10C and 45C. A bum-in test, 
centrifuging and temperature cycling testing were also performed 
on the flight devices. 

A similarly fabricated p-FET was flown on STRV-lb, 
launched June 17, 1994 into a 10.59-hour geosynchronous trans- 
fer orbit with an apogee of 36,000 km, a perigee of 300 km and 
an inclination of 7° south. Four p-FETs were located across the 
bottom of the JPL experiment box. Two p-FETs termed exposed 


(E) had an 80-mil (200 pm) A1 shield and two p-FETs termed 
shielded (S) had a 200-mil (560 pm) A1 shield. The p-FET 
temperature coefficient VA T was determined from the slope of 
the flight data as shown in Fig. 9 where VA t = 3.5 mV/°C. The 
dose for the STRV-lb is shown in Fig. 10; the data shown in Fig. 
10 was fitted with a linear line using the least squares method. 
The dose rate for exposed p-FET P4E is 96.0 rads/day and P8E is 
75.6 rads/day and the dose rate for the shielded p-FET P8S is 8.0 
rads/day and for P4S is 6.2 rads/day. p-FETs P4E and P8S were 
located near the outside of the box and P4S and P8E were located 
near the center of the box. Notice that the dose rate for the center 
p-FETs is 78 percent of the outside p-FETs. This shows the dose 
is different across the JPL experiment box. 

The p-FET dose rates determined from the data shown in Fig. 
10 are compared to the NASA models [9] as shown in Fig. 
1 1 . The data falls between the electron and proton data and has 
a slope that is closer to the proton data. There is some concern 
that the experimental data does not follow either the proton or 
electron predicted behavior. New experiments are being pre- 
pared to acquire space data to clarify this point. The Martin 
marietta dosimeter is expected to have very little temperature 
dependence, but if any is observed it can be used to adjust the 
value of VA t determined from ground calibration. Additional 
information on the STR-lb experiment is provided in [10,1 1]. 


Table 3 

p-FET Design Parameters 


Parameter 

Description 

Units 

Value 

To 

Reference value of temperature 

°C 

27 

Tm 

On-orbit value of temperature 

°C 

10 

VToo 

Threshold voltage @ zero dose & To 

V 

-0.8598 

vt t 

5VT/3T at constant dose D 

mV/°C 

1.7595 

AVT 

Saturation value of change in VT 
by dose D 

V 

-0.4998 

Do 

Value of D at which [VT] has 
increased by AVT[ 1-1/e] 

kilorads 

118.91 

poo 

Channel transconductance 

HA/V 2 

474.73 

Dm 

Dose coefficient of P in eq. (5) 

kilorads 

588.5 

n 

Temperature coefficient of p in eq. (5) 

- 

1.7196 

IDmo 

Design value of temperature - 
independent drain current 

uA 

prelim.: 88.36; 
final: 100 

VOmo 

Output voltage at ID^ 

V 

prelim.: -1.47 
final: -1.507 


3. SPACECRAFT IMPLEMENTATION 

As stated above, two surface charge monitors and two do- 
simeters are being installed on each of seven Martin Marietta 
commercial communications satellites. Total weight of the 
system is less than one pound (500 g) and the system requires 
approximately 0.5 watt of power. 

Surface charge monitors were installed on the north and anti- 
earth panels. The sensor located on the north panel will be con- 
tinuously shadowed for a period of six months; then it will be 
continuously exposed to sunlight during the other six months; it 
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Fi* 6 p-FET threshold voltage dose-temperature fit to Eq. 4 Fig. 9 STRV-lb p-FET VA- temperature dependence, 

G ' VAT = 3.5 mV/°C 



Fig. 7 p-FET beta dose temperature fit to Eq. 5 



Fig. 10 STRB-lb dose profile for devices behind 
80-mm and 220-mm A1 shields 



Fig. 8 p-FET output voltage dose dependence 
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Fig. 1 1 Comparison of STRV RADMON results 
with NASA radiation models 
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will provide data on seasonal variations of spacecraft charging 
potential. The sensor mounted on the anti-earth panel will be mov- 
ing in and out of the sunlight on the daily basis; it will provide 
data on daily variations in spacecraft charging potential. 

In addition, two dosimeters will be used to monitor internal 
charging fluxes and total dose accumulation. One dosimeter was 
mounted on the north panel internal to the spacecraft structure 
with the sensor facing the panel. In addition to its Kovar lid, this 
sensor is shielded from external radiation by the spacecraft panel. 
This location is typical for most spacecraft components. The sec- 
ond dosimeter is located outside the spacecraft structure on the 
anti-earth panel. In addition to the Kovar lid, this sensor is shielded 
from external radiation by a ten-layer Kapton thermal blanket. 
Thus a total of eight p-FETs with different shield thicknesses will 
provide electron flux count in eight energy channels. 

Support logic is combined with the sensors in order to meet 
the spacecraft’s interface requirements. The sensors were tested 
for functional performance following their installation on the 
spacecraft and during spacecraft-level thermal testing. 

4. CONCLUSIONS 

Development of mitigation techniques for spacecraft charg- 
ing control requires thorough knowledge of the environmental 
conditions and of their impact on spacecraft components. An 
integrated space environmental monitoring system has been 
developed, built, and installed on a geosynchronous communi- 
cations satellite. Its low weight and power requirements allow 
easy incorporation into the spacecraft design. Satellite perfor- 
mance can be correlated with in-situ charging and dose rate. 
Until now, the only environmental data available was from re- 
search satellites, in different time zones and along different flux 
tubes. We envision that in the future, such systems may become 
a standard part of the bus telemetry as common as temperature 
sensors for all geosynchronous spacecraft. 

The surface charging monitors must be mounted in carefully 
selected locations so that the degree of differential charging can 
be measured. In particular the differences in charging between 
the shaded and sunlit areas must be detected as the 3-axis-stabi- 
lized satellite gradually passes through the stages in its orbit. 
Fig. 12 shows one possible arrangement to accomplish this. 
Surface charge monitor # 1 , on the north panel, is always shaded 
between the fall and spring equinoxes and always illuminated 
the other half of the year. Charge monitor #2 is alternately illu- 
minated and shaded during each orbit. During this time, at any 
time except during eclipse, entry into a charged plasma cloud 
would give a transient differential reading between a shaded 
sensor and an illuminated sensor. A differential reading will also 
occur on emergence from eclipse when the satellite carries a 
net charge and one of the sensors is facing the sun. Thus the 
charging monitors will collect a significant datum that is not 
readily available from scientific satellites that carry particle 
detectors. The particle detector output readily reveals only the 
net charge on the satellite and not the differential surface charg- 
ing. Particle detectors, to provide net charge, would be desir- 
able but not feasible for a commercial communications satellite 
where weight and power are at a premium. 


When an in-situ indication of differential surface charging is 
available, anomalies can be correlated with this condition and 
with the deep-charging flux measured by the pair of dosim- 
eters. The instruments on CRRES provided similar correlation 
opportunities [12,13]. System anomalies are discussed in [12] 
and discharges, as registered by the Internal Discharge Monitor 
(IDM), are discussed in [13]. The system described in this pa- 
per is a monitoring tool, rather than a research tool, and is there- 
fore much less sophisticated than the collection of instruments 
available on CRRES. This simplicity and ease of accommoda- 
tion, however, is necessary to allow widespread use on com- 
mercial spacecraft. 
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Fig. 12 Internal and surface charge monitoring system. 
Internal charge monitors are identified as “dosimeters” 


Table 2 

The 250 krad time for each p-FET 
in the two dosimeters (SI units) 


Device 

PI 

P2 

PO 

P7 

Outside Dosimeter 

Lid Thickness (in (am K0VAR) 

0 

5 

150 

250 

Shield Thickness (in nm K0VAR) 

50 

50 

50 

50 

Total Thickness (in fim K0VAR) 

50 

100 

200 

300 

250 Krad Time (months) 

0.7 

1.6 

5.5 

13.3 

Inside Dosimeter 

Lid Thickness (in urn K0VAR) 

0 

50 

150 

250 

Shield Thickness (in (am K0VAR) 

250 

250 

250 

250 

Total Thickness (in nm K0VAR) 

250 

300 

410 

510 

250 Krad Time (months) 

8.9 

13.3 

26.4 

45.9 
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INTEGRATED CHARGE MONITOR 


June 4, 1996 


Martin G. Buehler 


Microdevices Laboratory 
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1. INTRODUCTION: This effort is a follow-on effort to develop the Integrated 
Charge Monitor, ICM, for Lockheed Martin for use as dosimeters on their 
communication satellites. The original ICM1 was delivered in September 15, 
1994. This effort differs from the first effort in that the first effort developed only 
a 1.2-um device; whereas, in this effort both 1.2- and 0.8-pm devices were 
fabricated. Also, the maximum radiation level was 160 krad of Co-60; whereas, 
in this effort the maximum radiation level was krad of 1-MeV electrons. Finally, 
the field oxide n-FET was replaced by a diode to facilitate on-chip theimometry. 


In this effort, 25 1.2-pm CMOS chips and 25 0.8-pm CMOS chips were designed 
and fabricated through MOSIS at the Hewlett-Packard facility, Corvallis, OR. The 
chips were packaged and tested at Halcyon Microelectronics, Irwindale, CA. The 
components were then electron irradiated and tested at JPL. The initial chip set 
consists of 25 1.2-pm CMOS devices designated ICM212 and 25 0.8-pm CMOS 
devices, designated ICM208. The ICM212 devices work extremely well; whereas, 
the ICM208 devices did not operate properly. Thus only 12 ICM2 devices were 
shipped to Lockheed Martin on May 20, 1996. 

2. ICM2 CHIP: The chip is 2.48 mm x 1.76 mm and contains the eight devices 
listed in Table 1. A schematic layout of the chip in the package is shown in Fig. 1 
and a photomicrograph is shown in Fig. 2. The chips were packaged in 16-pin 
flat packs using conductive epoxy and gold wire bonds. The wire bonding for the 
chip is shown in Fig. 1. Additional wire bonds were added, as shown in Fig. 3, to 
ground both the lid and the platform. The lid was attached to the package with 
conductive epoxy. The lid was altered to allow different shielding over each of 
the four p-FETs. The target shield thickness for each p-FET is listed in Table 1. 
The lid is thinned in two places over devices PGO and PG2 with a 25-mil (0.64-mm) 
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dimple. A 25-mil (0.64-mm) hole was formed over device PG1. Because the 
lid contains a hole, the packages are not hermetic. 

3. TEMPERATURE TESTS: The ICM2’s were tested electrically at room 

temperature prior to lidding so that rework, if needed, could be easily done. All 
25 ICM212s passed initial functional tests. Six ICM2 were measured at three 
different temperatures and a sample of the data is shown in Figs. 4 to 7. After 
electrical tests the devices were burned-in without bias at 125°C for 168 hours 
and centrifuged at 5000 g's. 

4. VT DISTRIBUTIONS: Electrical tests were again conducted at room 
temperature and the results are listed in Table 1. The threshold voltage, VT, 
distribution is shown in Fig. 8 The VT distribution is very tight so that the 
devices can be considered identical. The standard deviation for the VT is 5 
millivolts. This is important because it means that batch radiation calibration can 
be used. 

At this point 12 devices were shipped to Lockheed Martin on April 29 and 
electron irradiation began at the JPL dynamatron facility. In the course of this 
test it was discovered that the lids of the ICMs were not grounded. This caused 
the devices to fail. Thus the 12 devices previously shipped to Lockheed Martin 
were recalled and additional bond wires added as shown in Fig. 3. The devices 
were retested and shipped to Lockheed Martin on May 20. The devices shipped 
to Lockheed Martin are designated LM in Table 2. 

5. ELECTRON RADIATION CALIBRATION: The electron radiation test results are 
shown in Figs. 9 and 10. The tests were conducted at the JPL Dynamatron 
facility by Dr. Bruce Anspaugh. The results show the effect of the KOVAR lid in 
reducing the dose received by the p-FETs. The target Kovar lid thickness over 
each p-FET is given in Table 1. Measured values from four lids are listed in Table 
3. The electron irradiation was conducted using a 1-MeV electron beam that had 
an 8-inch (20-cm) diameter. The dose rate was 10 rads/sec for dose increments 
of 40 rads and 20 rads/sec for dose increments of 80 rads. The devices were 
irradiated at room temperature in a vacuum of IE-5 torr (1.33 x lO -6 kPa). The 
results show that the devices perform satisfactorily up to 600 krads of electron 
irradiation. 

The PG1 curve shown in Fig. 9 was fitted with an algorithm developed for the 
ICM1 dosimeters [1]. The following analysis is applicable to the unshielded 
device, PG1. The algorithm is based on radiation-induced filling of p-FET gate- 
oxide traps and the generation of oxide-silicon interface states. The expression 
for the output voltage at room temperature is: 

VO = VT 0 + AVT[1 - exp(-D/D v )] + > /[(2-ID/p 0 )-(1 + D/Dfc,)] (1) 
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where VT 0 is the initial threshold voltage, AVT is the trap filling voltage, D v is VT 
dose coefficient, ID is the temperature-independent operating current (100 pA), p Q 
is the initial transconductance, and Db is the p dose coefficient. In the analysis, 
the square root term was determined from the expression for the temperature 
independent operating point; see Eq. (17) in Ref 2: 


V[2ID/p o ]=V0 o -VT o (2) 

where V0 o is the preradiation value for VO. The fitting parameters are listed in 
Table 3. 

The parameters indicate that the trap filling process is relatively small in that AVT 
= 70 mV and that the process is completed quickly in that D v = 26.3 krad. The 

main factor that governs the shift in VO after 100 krads is interface-state 
generation. The above equation is the calibration curve for the dosimeter. 
Inserting the parameters listed in Table 1 in to Eq. 1 allows the calculation of VO 
at a given dose. 

The results in Figs. 9 and 10 show that the radiation dose does not decrease as 
the thickness of the shield increases. The shielding dependence at maximum 
dose for Figs. 9 and 10 is plotted in Fig. 11. Shielding thickness measurements 
for the four shields are listed in Table 4. The values were measured using a dial 
indicator. These results can be explained by the nature of the dimples fabricated 
above PG0 and PG2. These features are about 25 mils (0.64 mm) in diameter and 
are difficult to center over the FET’s located below. In the future these features 
should be increased to a diameter of 50 mils (1.28 mm). 

Annealing results are shown in Fig. 12. Radiation testing is an accelerated test in 
that a ten-year mission dose is applied in a few hours. It is well known that 
higher dose rates lead to greater annealing of the radiation damage. The choice 
of the dose rate is driven by the economic factor that the electron beam 
cost/hour is a constant. In this case the dynamatron cost was $200/hour. Thus, 
one wants to irradiate at the highest dose rate to save money. The annealing 
behavior shown in Fig. 12 indicates a 4 percent decay in VO over a 20 hour 
period. This is acceptable and will allow a dose rate of 20 rads/sec to be used in 
the future. 

6. IV CHARACTERISTICS: Typical current-voltage (IV) characteristics for the 

ICM2s fabricated in 1.2-pm CMOS are shown in Fig. 13. The p-FET characteristics 
overlap and are separated from the n-FET characteristics. In contrast to this 
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behavior, the IV characteristics for the ICM208 fabricated in 0.8-|jm CMOS are 
shown in Fig. 14. These results are not satisfactory because the p-FET 
characteristics do not overlap. A number of explanations were examined briefly 
including latch-up and Op Amp oscillation. More measurements are needed to 
determine the cause of these results. It is important to understand the cause for 
the difficulties with the ICM208 parts because the 1.2-pm CMOS process may not 
be available in the future. 

7. DISCUSSION: The 1.2 ^m p-FET’s have very tight VT’s where the standard 
deviation is 5 mV so batch radiation calibration can be used. One device was 
irradiated to 600 krads and was still operating properly. The electron irradiation 
showed that the trap filling process is completed by 100 krads and thereafter the 
process is characterized by interface-state creation. As seen in Fig. 11, the 
radiation dose does not decrease as the shield thickness increases for devices 
PGO and PG2. The lid dimples need to be increased from a diameter of 25 mils 
(0.64 mm) to 50 mils (1.28 mm) to ensure that the p-FETs are correctly centered 
under the shields. Also the shield thickness values need to be examined. The 
data from ICM1 flown on Telstar should be evaluated to see if there is sufficient 
spread in data from the four p-FETs. As seen in Fig. 11, thinner shields are 
needed to provide a better spread in the radiation response. After redesigning 
the shields, the effectiveness of the choices needs to be verified by electron 
radiation tests. 

The packages experienced a number of handling and testing steps. This caused 
some of the package leads to break as noted in Table 2. In this effort the 
packages were mounted and demounted in the chip carriers and tested at least 
four times. The insertion of the packages into the chip carriers and subsequent 
testing places considerable stress on the package leads. It is suggested that the 
leads of the packages chosen for flight be given a rigorous visual inspection. 

The development of ICM208 is imperative because the 1.2-pm CMOS process is 
an aging process. The redesign of ICM208 should include the use of internal test 
points to aid in chip probing in case of difficulty. Also certain metal features 
need to be designed for laser cutting during problem diagnosing. The electrical 
tests must include chip current measurements to assess the occurrence of latch- 
up. 
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Table 1. ICM2 Devices and Shield Tl 

Sickness 

DEVICE 

TYPE 

W(pm) 

L(pm) 

TARGET SHIELD 
THICKNESS KOVAR 
mil mm 

PGO 

PG1 

PG2 

DP3 

RM4 

NG5 

EX6 

PG7 

GATE OXIDE p-FET 
GATE OXIDE p-FET 
GATE OXIDE p-FET 
pn DIODE 

METAL-1 RESISTOR 
GATE OXIDE n-FET 
EXTERNAL RESISTOR 
GATE OXIDE p-FET 

193.54 

193.54 

193.54 

2.4 

193.54 

NA 

193.54 

9.6 

9.6 

9.6 

24214 

9.6 

NA 

9.6 

6 150 

0 (hole in lid) 

2 50 

10 250 

10 250 

10 250 

NA 

10 250 


Table 2. Test Results After Centrifuge and Burn-in. Devices used for temperature 
tests are noted by TEMP and those used for electron irradiation are noted by E- 


RAD. Devices delivered to Lock 

heed M 

artin are denoted by 

LM. 

DEV 


VT(V) 




wmm 


n 

NUMBER 

COMMENTS 

PGO 






DP5 

icm21201 

OK -Needs lid GNDed 

0.870 

0.870 

0.870 

0.867 

0.682 

3299 

1.02 

icm21202 

E-RAD 

0.866 

0.866 

0.864 

0.867 

0.671 

3384 

1.12 

icm21203 

P7 BAD 

0.846 

0.846 

0.847 


-0.248 

2953 

2.48 

icm21 204 

BAD 

0.148 

0.184 

0.161 

0.131 

0.142 

-953 

4.61 

icm21205 

TEMP, E-RAD 

0.867 

0.868 

0.868 

0.869 

0.674 

3367 

1.12 

icm21206 

GND pin bad-useable 

0.863 

0.862 

0.862 

0.863 

0.683 

3366 

1.12 

icm21207 

TEMP, LM 

0.862 

0.862 

0.863 

0.862 

0.668 

3331 


icm21208 

TEMP, LM 

0.874 

0.874 

0.868 

0.869 

0.680 

3272 

1.06 

icm21209 

OK- to be repackaged. 

0.869 


0.869 

0.865 

0.674 

3360 

1.05 

icm21210 

LM 

0.865 

0.865 

0.862 

0.863 

0.675 

3356 

1.11 

i cm 21 211 

LM 

0.867 

0.867 

0.860 

0.866 

0.675 

3334 

1.05 

icm21212 

E-RAD 

0.870 

0.870 

0.873 

0.871 

0.681 

3293 

1.09 

icm21213 

LM 

0.858 

IB 

0.858 

0.863 

0.670 

3388 

1.12 

icm21214 

TEMP, E-RAD 

0.881 

0.874 

0.875 

0.875 

0.669 

3287 

1.02 

icm21215 

LM 

0.873 

0.872 

0.873 

0.875 

0.672 

3316 

1.02 

icm21216 

E-RAD 

0.863 

0.863 

0.872 

0.869 

0.676 

3353 

1.11 

icm21217 

BAD VDD pin-rewire 

0.860 

0.860 

0.860 

0.860 

.668 

3361 


icm21218 

TEMP, P2 HIGH, LM 

0.876 

0.876 

0.895 

0.875 

0.670 

3272 

1.02 

icm21219 

OK - no lid 

0.870 

0.870 

0.867 

0.867 

0.663 

3308 


icm21220 

LM 

0.874 

0.874 

0.864 

0.866 

0.684 

3303 

1.02 

icm21221 

TEMP, LM 

0.867 

0.869 

0.868 

0.875 

0.677 

3343 

1.05 

icm21222 

LM 

0.868 

0.868 

0.872 

0.873 

0.681 

3291 

1.02 

icm21223 

LM 

0.872 

0.865 

0.861 

0.861 

0.681 

3371 

1.09 

icm21224 

LM 

0.879 

0.879 

0.874 

0.878 

0.667 

3309 

1.12 

icm21225 

BAD 

0.866 

0.866 

0.863 

0.865 

0.660 

3380 
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Table 3. Electron Calibration Parameters for ICM21216 


PARAMETER 

UNITS 

VALUE 

vo 0 

V 

1.504 

< 

H 

O 

V 

0.868 

V[2'D/P 0 ] 

V 

0.636 

AVT 

V 

0.070 

D v 

krad 

26.3 

D b 

krad 

95.0 


Table 4. Target and Measured Shield Thicknesses 


DEVICE NO. 

THICKNESS 
PG0 (mil, mjti) 

THICKNESS 
PG2 (mil, M*n) 

TARGET 

6.0 150 

2.0 51 

12 

5.2 130 

3.0 76 

16 

4.2 110 

3.3 84 

19 

5.0 130 

2.7 69 

25 

4.2 110 

3.0 76 



0P0 GND VDD GND OPI GND2 E* A2 

1 2 3 4 5 6 7 8 

Figure 1. ICM2 chip wire bonded to a 16-lead flat package. The hole in the lid is 
directly over p-FET PG1 which can be used to locate pin one of the package. 
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Correct alignment of the lid relative to the chip can be determined by looking 
through the hole over PG1 and observing the square p-FET PG1 and the words 
“12” above “ICM2”; see Fig. 2. 



Figure 2. Integrated charge monitor (ICM212) chip 2.48 mm x 1.76 mm. 



Figure 3. Cross section of the ICM2 package showing the additional wire bonds 
used to ground the platform and the KOVAR lid. 
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Cha2FETtemp 

ICM26308.XL 


OUTPUT VOLTAGE, |VO| (V) 


Figure 4. p-FET temperature characteristics showing the temperature 
independent point near 100 pA. 



Cha2FETtemp 

ICM26308.XL 


OUTPUT VOLTAGE, |VO| (V) 


Figure 5. n-FET temperature characteristics showing the temperature 
independent point near 150 pA. The n-FETs are useful in assessing failures at 
high dose levels. 
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Cha2REStemp 

ICM26308.XL 


OUTPUT VOLTAGE, VO (V) 


Figure 6. Metal-1 temperature characteristics showing a temperature coefficient 
of 3215 ppm/°C. This metal resistor serves as an on-chip reference device that is 
presumably independent of radiation. 



0.6 

Cha2tX)DBemp 

ICM26308.XLS 


0.7 0.8 

OUTPUT VOLTAGE, VO (V) 


0.9 


Figure 7. pn diode temperature characteristics showing a temperature sensitivity 
of -1.5 mV/°C at 10 pA. This diode serves as an on-chip thermometer. 
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CM 28308 .XLS THRESHOLD VOLTAGE, VTp (V) 

Figure 8. Room temperature (22°C) p-FET threshold voltage distributions from 19 
ICM2 packages containing four p-FETs, PO, PI, P2, and P7. The average and 
standard deviation for the distribution is 0.868 1 0.005 V. The standard deviation 
is as usual very tight at 5 mV allowing batch radiation calibration. A total of 76 p- 
FETs are shown on this chart. 
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Figure 9. 1-MeV electron irradiation response of device ICM21216 measured at 
100 pA, room temperature, and in a vacuum. 
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Figure 10. 1-MeV electron irradiation response of device ICM21212 measured at 
100 pA, room temperature, and in a vacuum. 



Figure 11. KOVAR shielding results for ICM21212 and ICM21216 at the maximum 
dose. For ICM21212, Dmax = 480 krads and for ICM21216, Dmax = 600 krad. The 
devices were irradiated with a 1-MeV electron beam. 
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Figure 12. Annealing results for ICM216 measured at 100 |jA after a 1-MeV 
electron irradiation at 20 rads/sec. 
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Figure 13. Typical IV characteristics for the four p-FETs and one n-FET from 
ICm21224. Note the p-FET characteristics overlap. 
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Figure 14. Typical IV characteristics for ICM208. 
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ABSTRACT: An array of four p-FETs was used 
to obtain dose-depth curves for the 100 MeV 
proton beam at the Loma Linda Proton Therapy 
Facility. The p-FET dosimetry compares 
extremely well with ionization chamber results 
and tends to resolve the peak and distal edge 
region more sharply. 

INTRODUCTION: The use of the p-FET (p-type 
Field Effect Transistor) for radiation dosimetry 
is based primarily on the charging of the gate 
oxide through the capture of holes created by 
the passage of high energy photons or 
particles. The holes add positive charge to the 
oxide which is detected as a shift in the FET 
threshold voltage. 

For many years p-FETs have been used for 
space flight dosimetry [1]. More recently, the 
devices were used for personal radiation 
dosimetry [2]. The device used here, is termed 
RADMON (Radiation Monitor), and consists of a 
multiplexed array of four p-FETs. It was 
developed for communication satellite 
dosimetry [3, 4]. This effort is directed at 
transferring the space radiation dosimetry to 
the medical community. 

The dose-depth curve for protons has a unique 
shape as shown in Fig. 1. The uniqueness has 
to do with the peak that occurs near the 
proton's end-of-range. X-ray dose-depth 
curves have a monitonic decreasing dose- 


depth curve. The data shown in Fig. 1 was 
obtained using a conventional ion chamber. 
The ion chamber dose-depth methodology is 
time consuming and, as will be shown, lacks 
the spatial resolution of p-FET dosimetry. 



lum25C2o.xls DEPTH IN WATER (cm) 

Figure 1. Ion chamber response to a 100 MeV 
proton beam passing through water. 

The proton response of an ion chamber, shown 
in Fig. 1, is considered the reference response 
which the RADMON must match. The dose is 
expressed in MU, measurement units. The 
critical features of this curve are its shape and 
the characteristic number which describes the 
peak to zero depth ratio. This ratio is 3.5. Thus, 
the objective of this effort is to evaluate the use 
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of an array of p-FETs each with its own shield 
for use in proton dosimetry and compare test 
results with conventional ion chamber 
dosimetry. 

RADMON measurements were taken at the 
Loma Linda Medical Center (LLUMC) using a 
2.54-cm diameter, 100 MeV proton beam which 
delivered a dose of 480 rads/min. The 
experiment took four hours during which 49 
data points were acquired. 

RADMON FABRICATION: The RADMON used in 
this experiment was designed previously as the 
Integrated Charge Monitor (ICM1) [4]. The 
RADMON, shown in Fig. 2, was refabricated by 
Hewlett Packard using their 1.2-um CMOS 
process with a gate-oxide thickness of 21 nm. 
The chip consists of the devices listed in Table 
1. The devices are addressed with a 
multiplexer (MUX) that is described later. 

The important devices for this study are the 
four p-FETs, P0, PI, P2, and P7 which are 
located on 1.25-mm centers as seen in Fig. 2. 
The devices are closed geometry FETs with the 
channel width, W, and length, L, as listed in 
Table 1. 



ICM54701 pcx VDD OPI GND2 


Figure 2. RADMON chip, 1.8 mm x 2.5 mm, 
fabricated in 1 .2-pm CMOS. p-FETs located on 
1.25-mm square at P0, PI, P2, and P7. The p- 
FET oxide is 20-nm thick. 

The chip was packaged in a 16-pin flat package 
as seen in Fig. 3 with 1-mil (25-nm) Au wire 
bonds. The chip is attached to the package's 
Au-plated Kovar platform with conductive 
epoxy. Au-wire bonds are used to ground the 


platform. Another Au wire bond was used to 
ground the lid of the package. This Au wire 
connects an internal package pad to an 
external package pad. Ag epoxy is applied to 
the Au ring that surrounds the top of the 
package. The combination of the Au ring, the 
Ag epoxy and the Au wire bond grounds the 
Kovar lid. This grounding procedure prevents 
the buildup of charge on an electrically floating 
lid which could cause an electrical discharge 
that damages or destroys the chip. 


Table 1. Device Dimensions 


NO 

DEVICE 

DESCRIPTION 

W 

pm 

L 

pm 

P0 

p-FET 

Poly Gate Oxide 

193.5 

9.6 

PI 

p-FET 

Poly Gate Oxide 

193.5 

9.6 

P2 

p-FET 

Poly Gate Oxide 

193.5 

9.6 

N3 

n-FET 

Metal Field Oxide 

193.5 

9.6 

R4 

Resistor 

Metal-1 

2.4 

24214 

N5 

n-FET 

Poly Gate Oxide 

193.5 

9.6 

E6 

External 

External device 

NA 

NA 

P7 

p-FET 

Poly Gate Oxide 

193.5 

9.6 



Figure 3. RADMON wire-bonded to 16-pin flat 
package. All portions of the package grounded 
to prevent charge build-up. 

The lidded device is shown in Fig. 4. The p- 
FETs are located at the corners of a 1.25-mm 
square. p-FET P0 has a 125-pm Kovar shield, 
PI has no shield, P2 has a 75-pm Kovar shield, 
and P7 has a 250-pm Kovar shield. The Kovar 
was chem-etched to thin the lid over P0 and P2. 
In the following measurements, it is assumed 
that the proton beam is uniform over the four p- 
FETs since the beam is 2.54 cm in diameter. 
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Figure 4. RADMON shielded by a KOVAR step 
Hd. 

SURROUND CIRCUITRY: The equivalent circuit 
used to test the RADMON is shown in Fig. 5. 
The p-FETs are placed in the feedback loop of 
an Op Amp and a 100 pA current is forced at the 
temperature-independent point [3]. The 
multiplexer is represented by two resistors, 
labeled by R mux , in the fi 9 ure - This a PP roach 
is needed because solid state switches have a 
non-negligible resistance. Op Amp U2 senses 
the potential at the p-FET and thus eliminates 
the need to know the value of Rrnux* ^ the 
measurement procedure, the output of U1, VO, 
is nulled by V Q ff and U2 amplifies VO by the 
ratio of R 2 to R-| which was set to 10. The 
output of U2, V am p, is measured by the 12-bit 
ADC. 

R 1 R 2 



Figure 5. Simplified circuitry for measuring the 
p-FET where U1 provides constant current via 
Vj and Rj. The output of U1 is nulled via V D ff 
and U2 amplifies VO via R-j and R 2 . Op Amps, 
U1 and U2, used in this experiment were dual 
SOIC LMC6042. 


EXPERIMENT: The experiment was designed to 
emulate ion chamber test conditions as closely 
as possible. The beam condition was identical 
for the ion chamber and RADMON tests. In the 
ion chamber experiment, the ion detector was 
moved away from the proton nozzle in a water 
tank. In the RADMON experiment, the p-FETs 
were fixed relative to the proton nozzle and 
various polystyrene shields were placed 
between the nozzle and the RADMON. 

The proton beam was 2.54 cm in diameter, had 
an energy of 100 MeV, and delivered a dose of 
480 rads/minute. In this experiment a run was 
defined as the proton beam being on for 2.5 
minutes; thus, the dose/run was 1200 rads/run. 
Throughout the 49 runs which comprised this 
experiment, the dose/run was constant. Three 
runs constitute a cycle during which one 
absorber was used. A cycle is illustrated in Fig. 
5. Electrical measurements were taken before 
and after each run thus each cycle consists of 
four electrical measurements. The cycles were 
characterized by a mean run number, N. The 
RADMON test time/run was about 10 seconds. 
The total time/cycle was less than 10 minutes. 
After the cycle was completed the absorber was 
changed and a new cycle begun. 

The equipment consisted of a Compaq 486/33 
MHz Aero lap computer, a National Instruments 
DAQPad-1 200, and a test box. The data 

acquisition program was written in Visual 
BASIC. During each cycle, the four electrical 
tests were run autonomously by the laptop with 
dwell times allowed for beam on time. 
Progress of the electrical tests was monitored 
by writing messages to the laptop screen in 
large letters which were visible on the control 
room monitor via the beam room camera. Thus, 
from the control room, the operator saw when 
the electrical tests were completed and turned 
the beam on. After completion of a cycle, the 
absorber was changed. Also the data, written 
to a floppy disk, was taken to the control room 
for immediate analysis. This allowed absorbers 
to be chosen for subsequent cycles so that 
gaps in the dose-depth curve could be filled. 

The data for a cycle is shown in Fig. 6. These 
data were fitted with a line using the least 
squares method. For the data shown in Fig. 2, 
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the slope of the line, dVamp/dN, is 40.48 mV/run 
and the standard deviation is 1 .6 mV/run. 



LUM25C20.XLS ELECTRICAL TEST RUN NUMBER, R 


Figure 6. Measurement cycle used to determine 
the "dose" behind an absorber. In this case the 
no absorber was used. The proton beam was 
on for 2.5 minutes between each of the four 
RADMON measurements. 

DATA ANALYSIS: The data acquired during 

this experiment was influenced by two factors: 

(a) The p-FETs were "dosed" by the proton 
beam during the experiment thus the p-FET 
output voltage Is a nonlinear function of dose. 

(b) The thickness of the Kovar lid over each of 
the four p-FETs is different thus the equivalent 
Kovar thickness must be added to the 
polystyrene absorber thickness. 

The FET dosing effect is modeled as follows. 
The FET detects radiation by absorbing 
radiation induced holes in oxide traps. The trap 
filling process is described by a rate equation 
for the number, n, of holes that fill oxide traps: 

dn/dt = a-<p(N t - n) ( 1 ) 

where a is the hole cross section, <p is the 
proton flux, N( is the total number of holes in 
the oxide, and t is time. 

The above equation is transformed to the 
following model equation: 

dVO/dN = (A VO - VO)/NO (2) 


where VO is the output voltage of the p-FET 
which is proportional to n, AVO is the shift in 
the FET output voltage which is proportional to 
N{, N is the mean run number which is 
proportional to t or dose, and NO is inversely 
proportional to a*<p. The solution to this 
equation is: 

VO = VOg + AVO[1 - exp(-N/N0)] (3) 

where VOg is the p-FET output voltage at N = 0 
and is about 1.5 V at 100 pA for the p-FETs used 
jn this experiment. This equation models the 
increase in the operating voltage of the p-FET 
with increasing N (or dose) and shows a 
decrease in the rise of VO at high dose where 
hole traps are filled. 

The slope of Eq. 3 is: 

dVO/dN = (A V O/N0)exp(-N/N0) (4) 

which shows that dVO/dN decreases with 
increasing dose due to trap filling. 

The effect of trap filling was removed by 
determining AV/NO for each cycle by 
rearranging Eq. 4: 

AVO/NO = (dVO/dN)exp(N/NO) (5) 

That is, the corrected slope, AV/NO, is 
determined from the measured slope dVO/dN 
times exp(N/N0). 

The implementation of this correction requires 
NO. For fitting purposes, the above equation 
was linearized as follows: 

ln[dVO/dN] = ln(AVO/NO) - N/NO (6) 

where dVO/dN and N are the measured 
quantities and AVO/NO and NO are coefficients 
to be determined. 

Multiple data points were acquired at 
polystyrene shield thicknesses of zero and 3.64 
cm. The data, plotted in Fig. 7, clearly shows 
the effects of radiation dosing on the 
performance of the p-FET. That is the initial 
slope measured at the beginning of the 
experiment is lower toward the end of the 
experiment. 
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The data shown in Fig. 7 were fitted 
simultaneously using the least squares method. 
The result yielded NO = 75.94. The RADMON 
dose-depth data were corrected using Eq. 5 and 
this value for NO. 



shield falls at a smaller polystyrene thickness 
and the distal edge for PI with no shield falls at 
a larger polystyrene thickness. This data 
shows the possibility that an array of p-FETs 
each with its own shield can be used to obtain 
the dose-depth curve including the shape in the 
distal edge region. 
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Figure 7. Determination of NO, the trap filling 
parameter, using the trap filling model where 
NO = 75.94. 

The second correction to the RADMON dose- 
depth curves involved determining the 
equivalent polystyrene thickness for the Kovar 
lid. This was done by shifting each of the dose 
depth curves until they overlapped the distal 
edge of the ion gage response. The factor for 
this adjustment was: 

H(mm)p 0 |y S ty r ene = 50*H(mm)Kovar ( 7 ) 

The results are listed in the key of Fig. 8. The 
maximum correction occurred for thickest 
Kovar lid of 0.25 mm which is equivalent to 1 .25 
mm of polystyrene. 

The results of the data analysis are shown in 
Fig. 8. These data also show the ion chamber 
data plotted in Fig. 1 where the vertical axis of 
Fig. 1 was multiplied by 18 to bring the results 
into coincidence. 

It is worth examining the RADMON dose-depth 
curves before the application of the step shield 
thickness correction. The data are shown in 
Fig. 9. The behavior in the distal edge region 
shows that the data are shifted properly with 
respect to the Kovar shield thickness. That is, 
the distal edge for P7 with a 250-pm Kovar 


LUM25C20 xls POLYSTYRENE EQUIV. H20 DEPTH, H (cm) 

Figure 8. RADMON and ion chamber dose- 
depth curves for a 100 MeV proton beam. The 
RADMON data was corrected for gate-oxide 
trap filling and for step shield thicknesses. 
Notice that the p-FET results cluster within the 
ion chamber response thus showing a sharper 
peak response for the p-FETs. 



IJUM 25 G 20 xls POLYSTYRENE EQUIV. H20 DEPTH H (cm) 

Figure 9. Data corrected for trap filling only. 
Effects of the Kovar step shield are clearly 
visible in the distal edge region. This indicates 
the possibility of profiling the distal edge region 
with a step shield and a single radiation 
sequence. 
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DISCUSSION: A comparison of the RADMON 
with the ion chamber, shown in Fig. 8, indicates 
good agreement in that there is a close 
correspondence between the two detectors for 
the 100 MeV proton beam. In fact, the RADMON 
appears to show that the peak region is sharper 
than the ion chamber data. This is expected 
since the ion chamber, which is 1 mm thick, is 
much thicker than the p-FET gate oxide, which 
is 20 nm thick. That is, the detecting region of 
the p-FET is 50,000 times thinner than the ion 
chamber. 

CONCLUSION: This experiment showed that 

dose-depth curves for both the ion chamber 
and the RADMON are very close. The 
experiment also showed that the p-FET data 
can be corrected for the effects of proton beam 
dosing and Kovar absorber thickness. In 
addition, the data in the distal edge region as 
seen in Fig. 5 showed that an array of 
individually-shielded p-FETs can be used to 
characterize the proton beam with one radiation 
sequence. 
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INTRODUCTION: This experiment is the second dosimetry experiment 

delivered to SAMMES (Space Active Modular Materials Experiments). The 
dosimeters were delivered in September 1994 and launched in June 1995 
on the aborted Pegasus vehicle. 


This experiment consists of p-FETs used to measure the radiation dose 
and capacitors used to monitor the build-up of charge in their dielectrics. 
These devices use the electronic circuitry designed to measure the 
resistance of the SAMMES actinometers. The p-FET dosimeters are custom 
designed p-FETs that are placed behind a number of shields with different 
thicknesses to verify the dose-depth curves used to design spacecraft 
shields. Currently, the dose-depth curves are based on computer codes 
that trace their roots to the NASA AP-8 and AE-8 developed in the 70's [1]. 

The capacitors used in this experiment are surface mount capacitors that 
are commonly used in the electronic circuits. As far as the experimenter 
knows, this is the first time that such capacitors have been flown in a test 
mode in space. The capacitors are being flown behind two different 
shields to study the effects of radiation on the capacitor's electrical 
properties. The understanding of deep dielectric charging is one of the 
major unsolved problems facing Earth orbiting spacecraft. A deep 
dielectric charging event is attributed to the demise of the ANIK [2] 
communications satellite in January 1994. It is a potential failure 

mechanism of the JPL experiment box flown on the STRV-lb. Deep 
dielectric charging was studied in the IDM (Internal Discharge Monitor) 
which flew on the CRRES (Combined Release and Radiation Effects 
Satellite) in 1990 [3]. In this experiment a number of dielectrics were 
placed under test and the build-up of charge was observed as pulses. The 
dielectrics consisted of electrical cables, printed wiring boards and Teflon 
and Solithane insulators. By far the largest number of discharges occurred 
in the PWBs composed of epoxy fiberglass. 
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OBJECTIVE: To measure the space radiation dose and deep dielectric 

charging behind different shields at six locations on SAMMES using p-FET 
dosimeters and surface mount capacitor arrays. It was required that these 
devices be measured by the same electronics used to measure the 
actinometers found in the other experiments included on SAMMES. 

APPROACH: This experiment consists of two types of packages 

designated the FET variant and the CAP variant. The description of the 
FET variant is found in Figs. 1 to 4 and the CAP variant is described in 
Figs. 5 to 8 . In this experiment, four FET variants and two CAP variants 
will be flown. The lid for the FET experiment, shown in Fig. 2, has holes 
over each FET. Different shields are placed over these holes as discussed 
later. Likewise, the lid for the CAP experiment, shown in Fig. 6 , has a large 
square hole to allow a maximum of radiation exposure for one of the CAP 
experiments. The other CAP experiment has a thick shield. 

PACKAGE FABRICATION: The 16-pin platform package is shown in Fig. 1 
and the part number is listed in Table 1 . Values for the resistors and 
capacitors are listed in Table 2. The package was modified by drilling two 
holes spaced 0.72 inch (1.8 cm) apart and tapped for a 2-56 mounting 
screw. The mounting screws serve two purposes: (a) to secure the devices 
to the spacecraft and (b) to electrically ground the package. An additional 
vent hole was drilled in the package to ensure that the shields experience 
no pressure differential during launch. THIS MEANS THAT THE PACKAGES 
ARE NOT HERMETICALLY SEALED AND NEED TO BE STORED IN A CLEAN 
DRY ENVIRONMENT. THE PACKAGES AND LIDS HAVE A BLACK MARK IN 
THE VACINITY OF PIN PT+ TO ENSURE CORRECT INSERTION OF THE 
DEVICES IN THE SOCKETS. The packages are numbered on the reverse 
side of the package. 

SUBSTRATE FABRICATION: The ceramic substrate for this experiment, 

shown in Figs. 1 and 5, is 25-mil (0.62-mm) thick AI 2 O 3 and was designed 

using AutoCAD. The substrate was designed to clear both the two screw 
holes and the vent hole. The substrate consists of a top side and backside 
Au-metal layers connected by Au-vias through the ceramic. After each 
layer is screen printed, it is fired at 500°C for 1 hour. The ceramic was 
attached to the package using non-conductive epoxy which was cured at 
150°C for 1 hour. The same ceramic substrate was used for both the FET 
and CAP experiments. 

COMPONENT ASSEMBLY: Non-conductive epoxy was used to attach the 

PRT to the ceramic. Conductive epoxy was used to attach the FETs, 
resistors, and capacitors, and PRT leads to the ceramic substrate. The 
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epoxy was cured at 150°C for 1 hour. The leads of the package were wire 
bonded to the substrate pads using 1-mil (0.02 mm) Au wire. 

SHIELD ASSEMBLY: The shields used in this experiment are aluminized 

mylar, aluminum sheet, and tantalum sheet. The thicknesses are listed in 
Table 3 and the assignment to the FET and CAP variant packages is listed 
in Table 4. For the FET variant the shields were cut into 200 mil x 500 mil 
strips and attached to the outside of the lids with J. B. Weld epoxy and 
cured at 150°C for 1 hour. For the CAP variant the Al-mylar shield was 
attached to the inside of the package and the 1-mm shield was attached to 
the outside of the package. The clearance between the top of the FETs and 
the underside of the shield is about 120 mils (3.0 mm). The 0.13-mm 
aluminum shield has the same equivalent density thickness as a typical 
solar-cell cover glass. 

FET PACKAGE: The FET variant package layout, shown in Fig. 1, indicates 
that the experiment consists of three p-FETs, F2, FI, and F3 and a positive 
resistance thermometer, PRT, noted as PT. The electrical schematic, 
shown in Fig. 2, indicates that the FET is connected in a Whetstone bridge 
circuit and surrounded by resistors, R1 and R2. This means that the FET is 
unbiased except for the short time that it is being tested. This is an 
important consideration for it is well known that FET radiation sensitivity is 
bias dependent. The measurement circuitry is shown in Fig. 4 and shows 
that the voltage across the bridge is connected to a log amplifier and then 
to a 12-bit ADC. During space measurements, Vref = 10 V so the FET 
current is about 100 pA which means that the FETs operate at their 
temperature independent point. 

CAP PACKAGE: The CAP variant package layout, shown in Fig. 5, 

indicates that the package contains one FET, F2, two capacitors, Cl and 
C3, and one resistance thermometer, PT. The electrical schematic, shown 
in Fig. 7, indicates that the charging and discharging current through the 
capacitor is measured from the voltage drop across R3 when Cl is 
measured and R4 when C3 is measured. Notice that the capacitor 
experiment has a bleeder resistor, R3 = 1 Mfi, across the circuit to prevent 
the unbounded build-up of charge which might damage the measurement 
circuit. The measurement circuitry is shown in Fig. 8. 

The measurement goal was to have a 1 second time constant for each 
capacitor circuit. Two commonly available capacitor types were chosen, 
X7R and NPO, and the largest available capacitor values were used. For 
the X7R, ten 0.1-pF capacitors are connected in parallel and this 
combination placed in series with a 1 MQ resistor. For the NPO, four 1-nF 
capacitors are connected in parallel and placed in series with a 250 MO 
resistor. Results of measurements, discussed below, are given in Table 6. 
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p-FET CHIP DESIGN AND MEASUREMENT: The p-FET chip used in this 

experiment is shown in Fig. 9. The chip was fabricated through MOSIS at 
Hewlett-Packard, Corvallis, OR in 1.2-pm n-well CMOS. The six chips were arranged 
in a 2x3 array and separated using a diamond saw. The MOSIS run number is N61X 
which had a gate oxide thickness of 20.9 nm. In this FET the source is internally 
connected to the n-well and the drain internally connected to the gate. Two p-FETs 
are fabricated on chip to provide redundancy in case one FET fails. The FETs were 
epoxied to the ceramic substrate with conducting epoxy. Wire bonds connect SW to 
V+ and DG and B to IR. 

Typical current-voltage characteristics for the FETs, shown in Fig. 10, indicate that 
the FET-to-FET variation is very small. This is important for it allows batch radiation 
and temperature calibration of the devices. The FET threshold voltage values, 
shown in Table 5, again confirms that the FET-to-FET threshold voltage variation is 
very small. 

The meaning of the comments in Tables 5 and 6 are: TEMP is temperature 
measurement, E-RAD is electron irradiation, FLT is flight, FLT SPARE is flight spare, 
and HOLE is a device where the ceramic is covering the vent hole. 

The current voltage characteristics for the FET, shown in Fig. 10, was determined 
using a circuit similar to the circuit shown in Fig. 4. The FET drain current is: 

ID(pA) = 1 0[Vref(V) - V+(V)] (1 ) 

where Vref and V+ are shown in Fig. 4. In this calculation the voltage drop across R5 
was ignored for it introduced less than a 3 percent error. The threshold voltage, VT, 
was determined from a least squares fit to the FET square law relation: 

ID = (p/2)(V+ - VT}2 (2) 

where p is the FET transconductance. As seen in Fig 10, the IV points are well 
approximated by the square law equation. 

p-FET OPERATING CURRENT: Most semiconductor sensors have a significant 

temperature dependence. The p-FET dosimeters are not exception to this rule. 
Fortunately, the temperature effects in the p-FET can be canceled by operating it at 
its temperature independent operating point. Temperature curves, shown in Fig. 11, 
indicate that the temperature effects on the IV characteristics cancel for an operating 
current of ID 0 = 100 pA. 
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The measurement of ID 0 follows by measuring the following parameters [4]. The 
temperature dependence of the threshold voltage was determined from: 

VT = VT 0 + VT T (T - T 0 ) (3) 

where VTj is the threshold temperature coefficient and T is the absolute 
temperature. The temperature dependence of p was determined from 

P = Po(T/T 0 )- n (4) 

where n is the channel mobility temperature coefficient. The p temperature 
coefficient is determined from 


PTo = - n Po /T o ( 5 ) 

The temperature-independent drain current was calculated from 

IDq = 2Po 3 <-VT T /p To ) 2 (6) 

The temperature parameters and ID 0 values for the three p-FETs in 
package SA10 are listed in Table 7. The results show that a choice of ID 0 = 

100 pA is a good one. This will reduce the temperature effects and allow 
the dosimetry effects to dominate the response. To obtain ID 0 = 100 pA, 

R1 = 100 kQ given that Vref = 10 V. 

DIFFERENTIAL FET OPERATION: The differential characteristics for the 

FET, shown in Fig. 12, reveal the nonlinear behavior of the FET leg of the 
bridge circuit and the linear behavior of the resistor leg. In developing the 
bridge circuit, R2 was chosen so that V+ is slightly more positive than V-. 
As seen in Fig. 12, V+ - V- is about 74 mV. THE EXACT VALUES NEED TO 
BE DETERMINED FROM PRE-FLIGHT GROUND TESTS. 

PRT MEASUREMENTS: For the resistance thermometer the current was 

determined from: 


l(pA) = 50 [Vref(V) - V+(V)] (7) 

where R5 = 20 kQ was used. The resistance was determined from a least 
square fit to the slope of the IV data. A typical example is shown in Fig. 13. 
In the future the value of this resistor should be reduced to 5 kQ to 
improve the accuracy of the current measurement. 
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CAP MEASUREMENTSi Th© capacitance circuit W3s measured by varying 
Vref as shown in Fig. 14. The voltage, Vref, was allowed a delay of about 5 
sec before each transition. The voltage was measured using a National 
instruments DAQPad-1200. The ADC has an input impedance of 1 MQ but 
the sampling time is less than 10 ps which has a negligible effect on the 
measurement. The response for the capacitors Cl (X7R) and for C3 (NPO) 
are shown in Fig. 15. As expected the capacitors have approximately the 
same time response of about 1 second. The capacitor response is 
measured at V- as seen in Fig. 8. The rising response, V r , across the 
capacitor is 


V r = V roo (1-exp(t/x r ) (8) 

where x r is rising time constant and V roo is the voltage at the end of the 

delay cycle. The rising time constant was calculated from a least squares 
fit to: 


t = -x r ln(1 -V r /V roo ) (9) 

where the log time response is shown in Fig. 16. For the falling response, 
Vf, is given by 


v f ■ Vfo'expH^f) (10) 

where xf is the falling time constant and Vf Q is the initial voltage. The 
falling time constant was calculated from a least squares fit to: 

t = xfln(V fo /Vf) ( 11 ) 

where the log time response is shown in Fig. 16. The rise and fall time 
constants are listed in Table 6. 

The temperature dependence of the ceramic surface mount capacitors is 
shown in Fig. 17. As expected [5] the X7R can have a +/- 15% change 
between -55 to 125°C; whereas, the NPO has a +/- 15 ppm/°C temperature 
coefficient. Thus, the X7R has a much larger temperature dependence than 
the NPO capacitor. The failure of the NPO to pull up to the rail is not 
understood. It indicates that the capacitor has started to leak. This same 
behavior has been observed in one other NPO capacitor so the effect 
appears to be real. 
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ELECTRON RADIATION CALIBRATION: The electron radiation test results 
are shown in Fig. 18. The tests were conducted at the JPL Dynamatron 
facility by Dr. Bruce Anspaugh. The electron irradiation was conducted 
using a 1-MeV electron beam that had an 8-inch (20-cm) diameter. The 
dose rate was 10 rads/sec. The devices were irradiated at room 
temperature in a vacuum of IE-5 torr. 

The results show the effect of the shields in reducing the dose received by 
the p-FETs. This data was fitted with an algorithm developed for the ICM1 
dosimeters [6]. The algorithm is based on radiation-induced filling of p- 
FET gate-oxide traps and the generation of oxide-silicon interface states. 
The expression for the p-FET output voltage at room temperature is: 

VO = |VT 0 | + AVT[1 - exp(-D/D v )] + (VO G - |VT 0 |)V(1 + D/D b ) (12) 

where VT 0 is the initial threshold voltage, AVT is the trap filling voltage, D v 
is VT dose coefficient, p Q is the initial transconductance, and D b is the p 
dose coefficient. When operated at the temperature-independent current, 
ID 0 : 


V [2 ID 0 /P 0 ] = V0 0 - |VT 0 | (13) 

where VO Q is the preradiation value for VO. The fitting parameters are 
listed in Table 8. 

This data was further analyzed after realizing that the FET response behind 
the 2-mm aluminum should have been zero for this shield is thicker than 
the 1.46 mm range of a 1-MeV electron beam in aluminum [7]. The range 
and linear energy transfer (LET) for various energies of electrons in 
aluminum are shown in Fig. 19. Thus, the background dose was 
subtracted from the p-FET shielded by the Al-mylar and the 0.25-mm Al- 
shield. The results are listed in Table 9. The results indicate that 
interface-state generation is a weak effect in that D b is large being near 3 
Mrad. The DVT is approximately the same which indicates that each FET 
has the same number of traps. The D v is larger for pFETI since it is more 

heavily shielded than PFET2. 

The initial dose rate follows from Eq. 12 for D b = qo and is: 

dVO/dD| D=0 = AVT/D v (14) 
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The dose rate is listed in Table 9. The value for PFET2 is 4.7 mV/krad 
which agrees with the initial slope of 4.7 mV/krad for the Co-60 results 
shown in Fig. 8 of Ref. 2. The dose rate for PFET1 at 2.3 mV/krad, is too 
low for a 0.25-mm At shield. It should be closer to 4.0 mV/krad according 
to the dose-depth curve shown in Fig. 11.3 of Ref. 8. 

The capacitors did not respond to the electron irradiation. Thus, they will 
serve as reference circuits which should not change in space. However, 
the capacitors have not been irradiated with protons which may cause 
some change to the charging and discharging characteristics. 

DISCUSSION: This experiment will develop the data to aid in developing 
design curves for spacecraft shielding against total dose. However, more 
ground tests are needed using electron irradiation to confirm the dose 
depth curves calculated previously [8]. The dose seen by the 0.25-mm 
shielded pFET is too low by almost a factor of two. 
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Table 1. SAMMES Dosimeter Parts List. 


ITEM 

DIMENSIONS 

MANUFACTURER 

P/N 

PRT: 0.39Q/°C 

2mm x 2.3mm 

OMEGA 

F3105 

Package: 16-pin Au-Plated 
Kovar Platform 


AIRPAC 

Manuf: HES Intern. 

ADI-1 6D-0006 

Lid: Ni-Plated Cold Rolled 
Steel 


AIRPAC 

Manuf: A. T. Wall 

ADC-1 6D-0001 

SMR: 16.9 kQ 1% Metal Film 

120 mil x 60 mil 

Panasonic 

ERJ 6ENF 1692 

SMR: 100 kQ 1% Metal Film 

80 mil x 50 mil 
(2.0 x 1.3 nm) 

Panasonic 

ERJ 6ENF 1003 

SMR: 1 MQ 1% Metal Film 

120 mil x 60 mil 
(3.0 x 1.5 mm) 

Panasonic 

ERJ 6ENF 1004 

SMR: 250 MQ 1% Metal Film 

12 orggi 


MC101822506F 

SMC: 0.1 pF XR7 

I3BEBH 

Panasonic 

ECU-V1H102JCX 

SMC: 1 nF NPO 


Panasonic 

ECU-V1H104KBW 

p-FET Chip: CMOS 

2.2 mm x 2.2 mm 

Hewlett Packard 

1.2-pm CMOS 

Ceramic Substrate: AI203 

mmmu 

Halcyon 

custom 


Table 2. R and C Values 


COMPONENT 

VALUE 

UNITS 

R1 


kQ 

R2 

16.9 

kQ 

R3 

1 

MQ 

R4 


MQ 

Cl 


mf 

Cl Total 

1 

MF 

C3 

1 

nF 

C3Total 

4 

nF 
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ble 3. SAMMES Dosimeters Radiation Shields. 


Material 

rgjgm 



mil 

H39 

Al-Myiar 






Al 

2.7 

35,685-9 

0.05 

1.97 

1.97 

Al 

2.7 

32,686-0 

0.13 

5.12 

5.12 

Al 

2.7 

22,685-2 

0.25 

9.84 

9.84 

Al 

2.7 

26,657-4 

0.5 

19.69 

19.69 

Al 

2.7 

26,695-7 

1 

39.37 

39.37 

Al 

2.7 

26,659-0 

2 

78.74 

78.74 

Ta 

16.6 

35,724-3 

0.05 

1.97 

12.10 

Ta 

16.6 

26,289-7 

0.25 

9.84 

60.51 

Ta 

16.6 

35,725-1 

0.5 

19.69 

121.03 

Ta ; 

16.6 

26,288-9 

1 

39.37 

242.05 


Tal 

t>le 4. 

SAMMES Radiation Shield Matrix 

UNIT 

LOC 

v/wam 

SHIELD | 

1 

2 

FET 



1 

FET 

E-T.TWK\rwi>vrmrrxmmmmmm 


3 

FET 


2 

2 

FET 



1 

FET 



3 

FET 


3 

2 

FET 

1.97-mil Ta (0.050 mm) 1 


1 

FET 

WMSSESnSEMTSES^t^m 


3 

FET 

39.37-mil Ta (1.000 mm) I 

4 

2 

FET 



1 

FET 

39.37-mil Al (1.000 mm) 


3 

FET 

19.68-mil Ta (0.5000 mm) 

5 

2 

FET 

MESESSMSSSEMMSSSmmm 


1 


1.00-mil Al (0.0254 mm) I 


3 

E32EH31 


6 

2 




1 

IsLUU&M 



3 

EMM 

U It R.WBHM 
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Table 5. FET Variant 


VT(F2) 


COMMENT 


E-RAD 

TEMP 


El 

■II 


TEMP, FLT SPARE 
BAD 

TEMP 

BAD 


TEMP 


E-RAD 

E-RAD 

HOLE 

FLT 

FLT SPARE 
FLT SPARE 
FLT 
HOLE 
FLT 

FLT SPARE 
E-RAD 


V 


883 

875 


0.874 

0.866 

0.880 

0.884 

0.873 

0.873 

0.873 

0.881 

0.879 

0.868 

0.879 

0.873 

0.881 

0.883 

0.875 

0.871 

0.870 

0.876 

0.873 

0.876 

0.882 

0.872 

1.549 

0.866 

0.873 

0.867 

0.878 

0.877 

0.877 

0.877 

0.877 

0.877 


Values 


VT(F1) 

V 


.874 

.872 

0.880 

0.874 

0.878 

0.878 

0.876 

0.877 

0.879 

0.877 

0.870 

0.873 

0.867 

0.877 

0.878 

0.880 

0.881 

0.870 

0.878 

0.878 

0.877 

0.876 

0.874 

0.875 

0.873 

0.872 

0.869 

0.876 

0.873 

0.873 

0.873 

0.872 

0.876 

0.873 



VT(F3) 

V 


0.871 

0.874 

0.878 

0.877 

0.876 

0.870 

0.832 

0.867 

0.873 

#NUM! 

0.873 

0.876 

0.871 

0.870 

0.876 

0.869 

0.874 

0.874 

0.873 

0.881 

0.879 

0.869 

0.873 

0.871 

0.884 

0.869 

0.871 

0.866 

0.871 

0.875 

0.863 

0.879 


R(PT) 

Q 


112.8 

111.2 

108.8 

112.8 

293.9 

106.5 

111.2 

108.3 

106.4 
111.2 
111.2 

106.5 

110.7 

112.7 
106.5 
111.2 
112.1 

108.3 

112.7 

133.3 
106.5 

106.4 

112.2 

112.8 
112.8 
112.8 
112.8 
112.8 
111.2 
112.8 
112.8 

108.4 


DAT6 

DAT6 

DAT6 

DAT6 

DAT6 

DAT6 

DAT6 

DAT6 


SA43 

SA44 

SA45 

SA46 

SA47 

SA48 

SA49 


Table 6. CAP Variant Values 


VT(F2) R(PT) 


COMMENTS 


.T SPARE, HOL 
E-RAD 
FLT 
FLT 

FLT SPARE 
FET=?, HOLE 



111.2 
111.2 
111.2 
.872 112.8 

.867 112.8 


0.866 

0.868 

0.842 

0.802 


1.032 

0.939 

0.906 


1.024 

0.944 

0.906 

1 V 

DAT6 

1.040 

0.936 

0.904 

DAT6 

0.987 

0.953 

0.922 

DAT6 

— 

— 

— 

DAT6 

0.992 

0.934 

0.904 

DAT6 
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Table 7. Te mperature Parameters and Operating Curre nt for SA10 


| SA10 

UNITS 

PFET2 

PFET1 

PFET3 


pm 

193.5 

193.5 

193.5 

Wu 

|jm 

9.6 

9.6 

9.6 

■■ 

K 

300 


300 

o 

h- 

> 

V 

-0.875 


-0.870 

VTj 

mV/°C 

1.69 

1.72 

1.72 

n 

unitless 

1.468 

1.572 

1.541 

Po 

HA/V 2 

494.6 

491.9 

489.7 

kp 0 

ma/v 2 

24.5 

24.4 

24.3 

PTo 

ma/v 2o c 

-2.42 

-2.58 

-2.52 

'Do 

ma 

117.9 

106.2 



Table 8. Electron 


Parameters for SA28 (with backgrou nd) 


PARAMETER 

UNITS 

PFET2 

PFET1 

PFET3 

Al-shield 


0.025 

0.25 

2.00 

vo 0 

n 

1.465 

1.455 

1.462 

|VT 0 | 


0.882 

0.874 

0.879 

V[2ID 0 /Po] 


0.583 

0.581 

0.583 

AVT 


336 

258 

40 

D v 

1 

60 

81 

67 

D b 

1 

314 

258 

236 

dV0 0 /dD 

mV/krad 

6.6 

43 

1.8 


Table 9. Electron Calibration Parameters for SA28 (Background removed) 


PARAMETER 

UNITS 

PFET2 

PFET1 

Al-shield 

mm 

0.025 

0.25 

V0 0 

V 

1.465 

1.455 

|VT 0 | 

V 

0.882 

0.874 

V[2-ID 0 /p 0 ] 

V 

0.583 

0.581 

AVT 

mV 

238 

176 

D v 

krad 

51 

76 

D b 

krad 

2177 

2840 

dV0 0 /dD 

mV/krad 

4.7 

2.3 
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SAM26623JWG 

Figure 1. FET variant CAD layout showing the 16-pin platform package, two-layer 
co-fired ceramic substrate and components including the FETs, resistors, and 
PRT. 



Figure 2. FET variant lid consisting of 10-mil (0.25-mm) stainless steel with 280- 
mil (7.1 -mm) holes punched over the three FET locations. 
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Figure 3. FET variant schematic showing three FETs, F2, FI, and F3 under test. 
R5 = 2.46kohm I 71 MEASURE 



lohm-IOkohm •— 4mV-8V 5mV-10V-* *- 2mV-3V 


Figure 4. FET connected to the spacecraft measurement electronics showing the 
Whetstone bridge nature of the p-FET dosimeter. 
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SAM26624.DVG 


Figure 5. CAP variant CAD layout showing the platform package, the two-layer 
co-fired ceramic substrate, and the components including FETs, resistors, 
capacitors, and PRT. 







Figure 7. CAP variant schematic showing the FET and the two capacitors, Cl and 
C3. 



Figure 8. CAP connected to the spacecraft measurement electronics showing that 
the analog MUX is connected across resistor R4 which is in series with the test 
capacitor. 
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Figure 9. p-FET chip fabricated in 1.2-pm n-well CMOS. The chip contains two 
closed-geometry p-FETs plus a resistor and two n-FETs. Contacts to the p-FET 
are SW = p-Source/n-Well, DG = p-Drain/Gate, and B = Body contact to the p- 
substrate. The internal DG connection forces the device to operate in its 
saturation region and provides electrostatic discharge protection to the gate 
oxide. 
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Chal PFET1 
SAM16501.XLS 


V+(V) 


Figure 10. Typical IV characteristics for the p-FETs from FET variant showing 
that all three FET characteristics overlap. The key lists VT in Volts and p in pAA/2 
for three FETs. 
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Cha2FETtemp 

SAM16503.XLS 


1.1 1.2 1.3 1.4 

OUTPUT VOLTAGE, V+ (V) 


1.5 


Figure 11. Temperature dependence of p-FET F2 in SA10. The key lists VT in 
Volts and (3 in pAA/2 at three temperatures. 



SAM16501.XLS VO (V) 


Figure 12. Typical IV characteristics for the PRT. The key lists the resistance in 
ohms and the offset current in pA. 
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Figure 13. Operational mode for the p-FETs is at Vref = 10 V. On the PFET line 
VT has units of Volts and 0 pA/V? and on the R2 line V+ - V- has units of mV. 


CAP RESPONSE, Cl (X7R) 


CAP RESPONSE, C3 (NPO) 



Figure 14. Capacitor experiment voltage profile for Vref and V-. The delay times 
are about 5 seconds and data points are acquired about every second. 
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ChalCAP 

SAM16501.XLS 


TIME, t (aprox sec) 


Figure 15. Typical capacitor transient characteristics for the CAP variant. Note 
that this is an Excel line chart where the 5-second delay times are compressed to 
one second. 



1 3 5 7 9 11 13 15 17 19 21 23 

chaicapLog TIME, t (min) 

SAM16501.XLS 

Figure 16. Logarithmic presentation of the data shown in Fig. 15 
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SA10 



Cha 2 CARemp 

SAM16503XLS TIME, t (approx, sec) 

Figure 17. Temperature dependence of the two capacitors in the CAP variant. As 
expected the X7R showed a larger temperature dependence than the NPO. Note 
that the NPO failed to pull-up at the highest temperature of 45.2°C. 



Figure 18. Radiation response from the three FETs in package SA28 with a Unit 1 
type shield to 1.0 MeV electron beam at a dose rate of 10 rads/sec. 
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Figure 19. Electron range and LET in aluminum. 
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Abstract 

A stacked p-FET dosimeter consisting of a 
RADMON mother chip with three p-FETs and 
multiplexer and an attached RADFET has been 
developed for the STRV-2/MWIR detector. 
Calibration of the dosimeter, using an Am-241 
source, indicates that the RADFET is about 20 times 
more sensitive to radiation than the RADMON. 
This dosimeter is expected to measure radiation 
dose from rads to Megarads. 


I. INTRODUCTION 

This dosimeter has four p-FETs for monitoring the 
radiation dose inside the Space Technology 
Research Vehicle (STRV-2) Medium-Wavelength 
Infrared (MWIR) detector electronics box to be 
launched in 1998. The three year mission dose for a 
450 km by 1600 km orbit behind 2.5-mm A1 shield is 
estimated to be 10 krads. The dosimeter, seen in 
Fig. 1, consists of a thin-oxide RADMON mother 
chip [1] with three p-FETs, multiplexer, and other 
devices. A thick-oxide RADFET [2, 3, 4] is 
attached to the mother chip. The device operates on 
the principle that radiation passing through the gate 

© IEEE 1996 


oxide region of the p-FET creates hole-electron 
pairs. Since the oxides have mainly hole traps, 
radiation-induced positive charge is trapped in the 
oxide. This build-up of positive charge changes the 
threshold voltage of the device which is interpreted 
in terms of radiation dose. 



Figure 1. RADFET mounted on the RADMON 
mother chip in a 16-pin flat package. Three 
radiation sensitivities are achieved: PK6 has highest 
sensitivity, PG0 and PG7 have moderate sensitivity, 
and PG4 has lowest sensitivity being shielded by the 
RADFET. 
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Since the RADFET is about 20 times more sensitive 
to radiation than the RADMON, this dosimeter can 
measure both low and high doses. The objective of 
this effort is to develop a stacked p-FET dosimeter 
for the STRV-2/MWIR detector and to calibrate the 
dosimeter using a low-cost bench-top radiation 
241 Am source. 


II CHIP FABRICATION 

The mother chip, shown in Fig. 2, was fabricated 
using the Hewlett-Packard 1.2-pm n-well CMOS 
process. The chip consists of an 8-bit multiplexer 
(MUX) that addresses four p-FETs (PGO, PG4, 
PK6, PG7), two n-FETs (NG1 and NF3), a resistor 
(RM2), and a user selected device (RX5) which can 
be the on-chip diode (DP5). 



Figure 2. RADMON mother chip 1.8 mm x 2.2 mm 
fabricated in 1.2-pm n-well CMOS. 


The RADMON p-FETs are closed geometry devices 
where the drain-gate is surrounded by the source- 
body. The FETs are covered with a metal-2 light 
shield. The RADFET p-FET has a comb structure. 


The critical dimensions of these devices are listed in 
Table 1. 



Table 1. 

Critical Device Dimensions 


NO. 



1! 

am 

X 0 

nm 

PGO 

p-FET 

Poly Gox 

193.5 

9.6 

20 

NG1 

n-FET 

Poly Gox 

193.5 

9.6 

20 

RM2 

Resistor 

Metal-1 

2.4 

24214 

NA 

NF3 

n-FET 

Metal Fox 

193.5 

9.6 


PG4 

p-FET 

Poly Gox 

193.5 

9.5 

20 

RX5 

Resistor 

External 

NA 

NA 

NA 

PK6 

p-RADFET 

Metal Gox 

500 

10 


PG7 

p-FET 

Poly Gox 

193.5 

9.6 



The untested RADMON chips were mounted in the 
16-pin flat pack using Ag epoxy as shown in Fig. 3. 
The three p-FETs on the mother chip are indicated 
in the figure. The RADFET was attached to the 
mother chip as shown in Fig. 1 with non-conductive 
epoxy in such a manner that PGO and PG7 are not 
covered and PG4 is covered by the RADFET. 



Figure 3. RADMON mother chip mounted in a 16- 
pin flat package. Three p-FETs are located on the 
RADMON 


The RADFET was fabricated using a thick oxide 
process at Southampton University Microelectronics 
Centre so it has a high radiation sensitivity [2, 3], 
The RADFET is 0.5-mm thick and provides some 
radiation shielding for devices MUX, NG1, NF3, 
and PG4 found on the mother chip. The RADFET 
consists of four p-FETs: two are very radiation 

sensitive Type R and two are the less radiation 
sensitive Type K. For the purposes here only one 
Type K is used. 
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Because the RADFET height exceeded the package 
height, an AI2O3 spacer was epoxied to the 
package. Grounding of the mylar lid, the Au ring, 
and the Au platform is provided by the Au wire 
bonds as shown in Fig. 4. A cross section view, 
shown in Fig. 5, provides additional detail. The 
grounding process is extremely important for any 
floating electrode can charge in a radiation field and 
discharge damaging or destroying the dosimeter. 

The package has a 25-pm thick aluminized mylar lid 
as seen in Fig. 6. This lid was chosen to minimize 
the radiation shielding since this dosimeter will be 
located in an electronics box and will receive about 
3.3 krads/year. The packages were burned-in at 168 
hrs at 125°C and centrifuged at 5000 g's. 



Figure 4. Au-wire bonds ground all portions of the 
package to prevent charge build-up 



Figure 5. Cross section of the MWIR Dosimeter 
package showing the wire bonds that ground the Au 
ring and aluminum-coated mylar lid. 



Figure 6. Cut-away view of the Al-mylar lid used to 
protect the chip and wire bonds from mechanical 
damage during board assembly. The A1 side is 
turned toward the package so that it is grounded. 

The electrical equivalent circuit for the chip is shown 
in Fig. 7, where the six on-chip devices are accessed 
through a 3-8 decoder multiplexer (MUX). The chip 
is connected to surround circuitry where the devices 
are connected between the input, OPI, and output, 
OPO, of an Op Amp. During measurement, current 
is forced through each device. 



Figure 7. Chip can address eight devices which are 
measured by an Op Amp connected between OPI 
and OPO. 
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The switches in the MUX are operated in pairs. One 
switch passes the measurement current and the other 
switch is connected to an Op Amp that senses the 
potential, V sen , at the device. 

The bias voltage for these devices is only applied 
during measurement. During irradiation the device 
is not powered and all leads are grounded. In 
operation as seen in Fig. 7, the voltage on the gate in 
the off mode, is the voltage across the n-well at the 
source of the p-FET. Leakage currents through the 
n-well junction keep this node at zero bias. 


III. THRESHOLD VOLTAGE DISTRIBUTIONS 

Integral to the dosimetry process is batch 
calibration. That is, a few devices must be selected 
from a batch and destructively tested with radiation 
with the assurance that the resulting calibration 
curve is applicable to the rest of the batch. This 
requires fabrication of p-FETs with tight electrical 
parameters. 

Room temperature current-voltage measurements 
are shown in Fig. 8 for the five FET’s in this 
experiment. These curves were modeled using the 
saturation drain current expression: 

ID = (p/2)(VO - |VT|) 2 (1) 

where VO is the FET output voltage, p = KP W/L is 
the transconductance, KP = pC 0 , W is the channel 
width, L is the channel length, p is the channel 
mobility, and C Q is the gate oxide capacitance per 
unit area. 

The threshold voltage, VT, is determined from the 
extrapolation of the curves shown in Fig. 8 to zero 
drain current using a least square fitting technique. 
The p values were determined from the slope of the 
characteristics shown in Fig. 8. 

VT distributions for the RADMON FET’s are shown 
in Fig. 9 and VT distributions for the RADFET’s are 
shown in Fig. 10. The VT for the RADMONs is 
0.846 ± 0.003 V and for the RADFET’s is 2.045 ± 
0.010 V. The standard deviation for all FET’s is 
considered to be acceptably small; in addition, no 
outliers were identified. 



ChaPFETI 
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OUTPUT VOLTAGE, VO (V) 


Figure 8. Typical current-voltage characteristics for 
the five FET’s located on the MWIR. Note that the 
FET characteristics for PG0, PG4, and PG7 overlap. 
The numbers in the key are VT in volts and p in 
pAA/ 2 


> 

0) 

*o 

<A 


> 

A 


CD 

O 

K 

CL 


D 

O 



chaicurrfVob THRESHOLD VOLTAGE, VTp (V) 

MW26524.XLS 

Figure 9. The cumulative probability (CUM. PROB) 
threshold voltage distributions for RADMONs; 39 
p-FET values are shown. 
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Table 2. MWIR Dosimeter FET Temperature Parameters at T n = 300 K. 


PARAMETER 

UNITS 

PGO 

NG1 

PG4 

PK6 

PG7 

DEVICE NO. 

— 

25 

25 

25 

25 

25 

VT 0 

V 

-0.8695 

0.7122 

-0.8681 

-2.0500 

-0.8710 

VT? 

mV/°C 

1.1325 

-0.9154 

1.2060 

1.7272 

1.1686 

n 

unitless 

1.0660 

1.1672 

1.0863 

1.3579 

1.0794 

Pa 

pAA/J 

466.3 

1343.1 

465.8 

249.5 

464.7 

k ¥> 0 

pAA/2 

23.13 

66.62 

23.10 

5.0 

23.05 

Pto 

pAA/ 2 °C 

-1.657 

-5.226 

-1.687 

-1.129 

-1.672 

IDn 

UA 

94.73 

148.71 

103.3 

72.65 

98.03 

DEVICE NO. 


13 

13 

13 

13 

13 

VT 0 

V 

-0.8674 

0.7108 

-0.8689 

-2.0368 

-0.8680 

VTj 

mvrc 

1.1094 

-0.8676 

1.1324 

1.6448 

1.1275 

n 

unitless 

1.0266 

1.1086 

1.0363 

1.2950 

1.0321 

Pa 

pAA/2 

465.1 

1330.5 

463.6 

253.7 

462.3 

K^ 0 

pAA/2 

23.07 

65.99 

22.99 

5.07 

22.93 

PTo 

pAA/ 2 °C 

-1.591 

-4.917 

-1.601 

-1.095 

-1.591 

JDq 

LtiA 

97.77 

146.66 

99.63 

73.66 

99.30 



cha2Cunftob THRESHOLD VOLTAGE, VTp (V) 

MWI26524.XLS 

Figure 10. The cumulative probability (CUM. 
PROB) threshold voltage distributions for 
RADFETs; 13 p-FET values are shown. 


IV. OPERATING CONDITIONS 

The operating current for these devices is 
determined from the temperature dependence of the 
IV characteristics. Data from devices PK6 and PG7 
are shown in Fig. 11. A close examination of this 
figure indicates that the IV characteristics cross at 
98 pA for PG7 and at 73 pA for PK6. The cross 
over point defines the temperature-independent 
operating point. 



mwr 5 A 27 .xls OUTPUT VOLTAGE, VO (V) 

Figure 11. Temperature dependence of the p-FETs 
in package MWIR25 showing the temperature- 
independent current of 98 pA for PG7 and 73 pA 
for PK6. The numbers in the key are VT in volts, p 
in pA/V 2 , and T in °C. 

The temperature dependence of the threshold 
voltage was determined from: 

VT = VT 0 + VT T (T-T 0 ) (2) 

where VTy is the threshold temperature coefficient. 
The temperature dependence of p was determined 
from 

P = Po(T/T 0 )- n (3) 


131 






where n is the channel mobility temperature 
coefficient and T is the absolute temperature. The p 
temperature coefficient is determined from 


PTo = ’ n 'Po /T o ( 4 ) 

The temperature-independent drain current was 
calculated from 

ID 0 = 2p 0 3(-VT T /p To ) 2 (5) 

The ID 0 values for devices 13 and 25 are listed in 
Table 2. In this experiment the operating drain 
current, ID 0 is 100 pA. This will reduce the 
temperature effects and allow the dosimetry effects 
to dominate the response. 


V. RADIATION CALIBRATION 

A low-cost bench-top approach was used to 
calibrate the devices using a 1-pCi 241 Am source. 
The Am source emits alpha particles in the 5-MeV 
range and a number of X-rays with a maximum 
energy of 59 keV. The dosing mechanism is 
currently under study. 

The p-FET voltage was measured at 100-pA. The 
RADFET (PK6) was determined directly from this 
voltage and the RADMON (PGO, PG4, and PG7) 
voltage was amplified 10 times. Typical results in 
Fig. 12 show the linear time dependence of the 
voltage for PK6 and PGO measured every 30-min. 

To compare RADFET and RADMON results, the 
differences in the device-to-radiation source distance 
must be determined. These two distances are 
illustrated in Fig. 13 as H-| and H 2 . 

A FET-to-source distance model was developed 
using the dose rate, Df in rad/sec expression [5]: 

D t =K 0 cpL (6) 

where K 0 is in rad-mg/MeV, (p is the flux in 
#/cm 2 •sec, and L is the linear energy transfer (LET) 
in MeV-cm 2 /mg. The flux for a point source is [5]: 

ip = S/4itR 2 (7) 

where S is the source strength in #/sec and R is the 
distance in mm from the source to the p-FET. 



Figure 12. RADFET (PK6) and RADMON (PGO) 
output voltage change due to a 1-pCi 24 1 Am source 
where the RADMON voltage is amplified 10 times 
and the source is 1.27-mm above the reference plate 
seen in Fig. 13. 



Figure 13. Cross section of the apparatus used to 
position the Am source above the dosimeters. 

The distances, H-| and H 2 , seen in Fig. 13, cannot 
be determined by direct measurement. However, 
they can be determined by making dose 
measurements at different heights, h, above the 
reference plate. Thus, R is defined as: 

R = H + h (8) 


where H is the FET-source distance and h is 
measured from the reference plate. The source is 
raised and lowered a distance h above the reference 
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plate seen in Fig. 13. Results from the RADMON 
are shown in Fig. 14. 
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Figure 14. RADMON (PG0) voltage change due to 
a 1-pCi 241 Am source. The RADMON voltage is 
amplified 10 times. The source is located at h with 
respect to the reference plate seen in Fig. 13. 



The results in Fig. 15 show the voltage shift 
relationship between the p-FETs as a function of R. 
The parameters determined from the least squares 
fit to the data are listed in Table 3. Details of the 
parameter extraction are given in the Appendix. 

The ratio of the dose response for the p-FETs is 
Kratio = K(PK6)/K(PG0&PG7) = 21.0. From 
6^Co tests, the RADMON Vq is 4.7 mV/krad [ 1 J. 
Thus using the K ra tj 0 = 21 .0 leads to a RADFET 
V D = 98.7 mV/krad. Previous response [3] of the 
Type K RADFET using 60c o gamma rays was 200 
mV/krad at zero bias and low dose. Recombination 
theory predicts lower response for lower photon 
energy. This is reflected in the LET factor initially 
presented in Eq. 6 and found in Eq. 11. At this 
time, the Am FET dosing is being compared to Co- 
60 dosing. 

The data for each curve in Fig. 14 was acquired 
over a two day period. Using Vq = 4.7 mV/krad, 
the dose rate for the h = 1.27 mm curve is 62 
millirads/sec and the dose rate for the h = -2.54 mm 
curve is 235 millirads/sec. These dose rates are 
higher than the expected space mission dose rate of 
0.1 millirads/sec but smaller than the common 
experimental dose rate of one rad/sec. 


Finally, the p-FET output voltage rate, 
related to the dose rate by: 

VO t , is 

Q 

Q 

> 

il 

■*—> 

o 

> 

(9) 

where Vp is the voltage-dose sensitivity in mV/rad. 
Combining the above equations leads to the fitting 
algorithm: 

VO t = K/(H + h) 2 

(10) 

where 


K = V D K 0 SL/4it 

(11) 

The ratio of the K’s for the FETs leads to: 


Vp(RADFET) _ K(RADFET) 
V D (RADMON) K(RADMON) 

(12) 


This relation allows the calculation of the dose 
sensitivity ratio for the two FET types found in the 
RADMON and RADFET where the difference in 
the FET-to-source distance has been removed. 



mwir613i.xls DISTANCE, R (mm) 

Figure 15. RADFET (PK6) and RADMON 
(PG0&PG7) output voltage shift rate as a function 
of the distance to a 1-pCi Am-241 source. 
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Table 3. RADFET and RADMON Parameters. 


SYM 

UNITS 

RADMON 

RADFET 

H 

mm 

5.536 

5.011 

K 

mVmm^/min 

0.091 

1.994 


VI. DISCUSSION 


The analysis for K and H follows from a least 
squares fit to the following set of equations derived 
from Eq. 1 0: 

hi = -Hi + VKi/VV0 t i (13) 

h 2 + AH = -Hi + VK 2 /VV0 t2 (14) 


The FET response from 1-pCi 241 Am source 
radiation is easily measured. A FET-to-source 
distance model allows determination of relative 
radiation response of RADFET and RADMON 
even though the FETs are located at a different 
distance from the source. The RADFET is 21.0 
times more sensitive to 241 Am than RADMON. 
Radiation doses from rads to Mrads are achievable 
with this detector. 


VII. APPENDIX 


The fitting procedure used to extract the parameter 
listed in Table 3 is illustrated by the graph shown in 
Fig. 16. The analysis is based on the 0.5-mm 
RADFET thickness and the 0.025-mm non- 
conductive epoxy thickness. Thus 

H2 - HI = AH = 0.525 mm (12) 
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Figure 16. Graphical representation of the 
parameter extraction process used to extract the K 
and H parameter listed in Table 3. 


The result of the fitting procedure is shown in Fig. 
16. 
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